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I.  INTRODUCTION 


This  report  contains  theoretical  investigations  on  the  electrical  conduc¬ 
tivity  and  thermodynamic  properties  of  nonideal  plasmas,  which  were  carried 
through  in  the  period  from  1  November  1979  to  31  December  1980  under  ONR  Contract 
N00014-79-C-0073.  In  addition,  the  theoretical  results  were  compared  with  experi¬ 
mental  data  for  nonideal  plasmas.  These  comparisons  are  of  a  preliminary  nature, 
since  the  experimental  conductivities  for  nonideal  plasmas  differ  not  only  quant¬ 
itatively  but  also  qualitatively  in  the  literature. 

CHAPTER  II .  The  dependence  of  the  electrical  conductivity  c  of  dense  (non- 

2  1/3 

degenerate)  plasmas  on  the  nonideality  parameters  y  *  Ze  n  /KT  was  evaluated 
by  summing  the  probabilities  for  v-body  interactions  (v  «  2,3,4,...)  of  the  con¬ 
duction  electrons.  It  is  shown  that  o  is  noticeably  smaller  than  the  binary  con¬ 
ductivity  o 2  for  y  >  10  1.  The  theoretical  decrease  of  o  with  increasing  y  is 
confirmed,  however,  only  by  some  experimental  data,  while  other  experimental  data 
indicate  an  increase  of  a  with  increasing  y  for  the  same  pressure. 

CHAPTER  III.  Based  on  the  classical  and  quantum  Boltzmann  equations,  the 
electrical  conductivities  of  classical  and  degenerate  nonideal  plasmas  were  evalu¬ 
ated.  Although  in  this  kinetic  approach  many-body  interactions  are  taken  into 
account  only  through  an  exponentially  shielded  Coulomb  potential,  in  which  the 
electron-ion  scattering  occurs,  the  results  give,  in  agreement  with  the  experimen¬ 
tal  data,  conductivities  which  are  by  about  one  order  of  magnitude  smaller  than  the 
Spitzer  conductivity  for  ideal  plasmas.  The  increase  of  the  dimensionless  conduc- 
tivity  a  =  m  e  o/(KT)  with  increasing  y  is  confirmed  by  some  experimental 
data  but  not  by  all  of  them.  The  new  Coulomb  logarithm  does  no  longer  go  to  zero 
for  large  y  values  (as  in  the  Spitzer  theory)  but  is  well  behaved  for  large  elect¬ 
ron  densities,  and  even  for  solid  state  densities  due  to  the  consideration  of 
electron  degeneracy. 


1 


1 


CHAPTER  IV«  With  Che  help  of  quantum-field  theoretical  methods  from  the 
theory  of  metals,  the  electrical  conductivity  of  nonideal  plasmas  was  calculated 
under  consideration  of  electron  scattering  by  low-frequency  plasflons  (ion  waves) 
and  high-frequency  plasmons  (electron  waves)  for  classical  and  degenerate  conditions. 
The  resulting  conductivity  formulas  agree  with  the  Spitzer  theory  for  y  -►  0  and  ex¬ 
hibit  numerical  values  which  are  considerably  smaller  than  the  Spitzer  values  but 

are  still  larger  than  the  theoretical  conductivities  obtained  In  III  for  increasing 

* 

y.  The  numerical  values  o  agree  with  the  experimental  data  qualitatively  but  are 
somewhat  too  high. 

CHAPTER  V.  The  possibility  of  anomalous  diffusion  and  conduction  transverse 

to  magnetic  fields  Bq  was  studied  since  large  charged  particle  transport  across 

magnetic  fields  is  of  interest  for  MHD  generators.  For  weakly  nonideal  plasmas, 

the  anomalous  transverse  conductivity  was  shown  to  be  a,  *  a>£/4ir/2w  where  u>  * 

“  r  o  F 

2  1  /2  , 

(ne  /e  m)  is  the  plasma  frequency  and  u)  *  | e | B  /m  is  the  gyration  frequency 

O  DO 

of  the  electrons  (e,m).  This  formula  agrees  with  experimental  data  for  weakly  non¬ 
ideal  plasmas,  but  should  be  also  correct  qualitatively  for  nonideal  plasmas.  There 
are,  however,  no  experimental  data  available  on  anomalous  diffusion  and  conduction 
in  magnetic  fields  for  nonideal  plasmas. 

CHAPTER  VI.  In  connection  with  the  electric  current  transport  in  the  elect¬ 
ric  field  fluctuations  produced  collectively  by  the  electrons  and  ions  in  random 
thermal  motion,  the  electric  microfield  distribution  of  thermal  plasmas  was  deri¬ 
ved  by  equilibrium  statistical  mechanics.  Comparison  with  the  resulting  tempe¬ 
rature  dependent  microfield  distribution  with  the  classical  (T- independent)  Holts- 
mark  distribution  and  its  later  extensions,  indicates  that  the  latter  theories  are 
approximately  applicable  to  strongly  nonideal  plasmas  but  are  invalid  for  ideal 
plasmas  (to  which  they  are  usually  applied  in  literature). 
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CHAPTERS  VII  -  VIII.  By  means  of  Bose  statistics,  the  contribution  of  the 


thermally  excited  (longitudinal)  electron  and  ion  waves  to  the  free  energy  of 
nonideal  classical  and  quantum  plasmas  was  calculated.  It  is  shown  that  the  ran¬ 
dom  low-frequency  ion  oscillations  contribute  more  to  the  free  energy  than  the 
high-frequency  electron  oscillations.  The  free  energy  of  the  random  ion  waves 

is  quantitatively  comparable  to  the  free  energy  of  the  thermal  (non-  collective) 

23  -3 

ion  motions  for  high  densities  (n  <10  cm  )  and  standard  plasma  temperatures 
(T  <  10^°K).  Similar  calculations  were  performed  for  dense  gases,  in  which  the 
random  sound  oscillations  lead,  however,  only  to  a  small  correction  of  the  free 
energy . 

The  theoretical  research  on  nonideal  plasmas  needs  further  clarifications 
by  experiments.  In  particular  more  reliable  conductivity  data  for  nonideal  al¬ 
kali  and  noble  gas  plasmas  are  needed.  This  is  a  preliminary  report  of  research 
results,  which  will  be  communicated  later  in  form  of  publications. 
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II.  CONDUCTIVITY  OF  NON I DEAL  PLASMAS  WITH  MANY -PARTICLE  INTERACTIONS 


By 

H.  E.  Wilhelm 


ABSTRACT 


The  dependence  of  the  electrical  conductivity  of  nondegenerate, 
dense  plasmas  on  the  nonideality  parameter,  y  ■  Ze  n  /KT  (ratio  of 
Coulomb  interaction  and  thermal  energies) ,  is  derived  by  summing  the 
probabilities  for  v-body  interactions  (v  *  2,3,4,...)  of  the  elec¬ 
trons.  As  an  application,  the  dimensionless  probability  coefficients 
for  binary  and  triple  Coulomb  interactions  are  calculated  by  means 
of  simple  physical  models,  and  a  conductivity  formula  for  moderately 
nonideal  plasmas  (0  <  y  <  1)  is  derived  in  which  all  parameters  are 
known.  The  theory  is  shown  to  agree  with  recent  experimental  data. 
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INTRODUCTION 


High  pressure  plasmas  (10*  bar  i  P  i  10^  bar)  produced  by  shock  wave 

compression  are  now  of  considerable  technical  interest.  A  large  number  of 

publications*  *^  are  concerned  with  the  measurement  of  the  anomalous  electrical 

conductivity  of  proper  nonideal  plasmas  (10  *  ~  y  &  1).  Theoretically,  however, 

only  the  conductivity  of  ideal  (y  0)  and  weakly  nonideal  (y  <<  1)  plasmas  is 

11  12 

adequately  understood.  ’  The  degree  of  nonideality  of  a  fully  ionized 

plasma  is  defined  by  the  interaction  parameter  y,  which  represents  the  ratio  of 

2  1/3 

average  Coulomb  interaction  (Ze  n  )  and  thermal  (KT)  energies  (n  *  electron 
density,  Z  =  ion  charge  number,  e  =  elementary  charge), 

y  =  Ze2n*/3/KT  =  1.670  x  10-3Zn*/3T-*  [e.s.u.]  .  (1) 

The  conductivity  theories  of  ideal  and  weakly  nonideal  (0  <  y  «  1)  plasmas 
break  down  for  y  >  10  *,  since  the  Debye  radius, 

D  -  [Z/4„a  +  Z^'V1'2-'1'3  -  ,-1/2n-1/3  ,  (2) 

loses  its  physical  meaning  as  an  electric  shielding  and  Coulomb  Interaction 
length.  This  is  seen  from  the  number  of  electrons  NQ  in  the  Debye  sphere  of  a 
scattering  ion,  which  is  no  longer  large  compared  with  one  for  y  >  10  *, 

Nd  =  ( Att / 3 )  [Z/4ir(l  +  Z)]3/V3/2  ~  y'3/2  .  (3) 

n  a  <i  / 

For  strongly  nonideal  conditions,  n  >  lO^cm  and  T  =  104°K,  we  have 
y  >  0.775,  D  <  4.881  x  10  cm,  and  <  4.87  x  10  2!  Another  reason  for  the 

inapplicability  of  the  conductivity  theory  of  ideal  and  weakly  nonideal  plasmas 
to  proper  nonideal  plasmas  is  the  standard  assumption  of  (shielded)  binary 
Coulomb  collisions  (v  =  2) ,  whereas,  in  reality,  the  conductivity  is  determined 
by  many-particle  interactions  (v  =  2,3,4,...)  for  y  >  10~*. 
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The  many-body  interaction  is  one  of  the  classical,  unsolved  problems  of 

physics.  For  this  reason,  we  calculate  the  conductivity  of  nonideal  plasmas 

and  the  probabilities  for  many-particle  interactions  by  means  of  dimensional 
13  14 

theory.  ’  This  approach  gives  the  exact  dependence  on  the  relevant  dimen- 

13  14 

sional  plasma  parameters  *  and  numerically  correct  results  up  to  a  dimen¬ 
sionless  coefficient,33’^  which  is  in  general  of  the  order  10^.  The  plasma 
is  assumed  to  be  fully  ionized  and  nondegenerate,  i.e., 

n  =  Zn±  <  n  ,  n  =  2(2irmKT/h2)3/2  =  4.828  x  1015T3/2  .  (4) 


ELECTRICAL  CONDUCTIVITY 


In  a  system  of  reference  in  which  magnetic  fields  are  absent,  a  linear 
electric  current  response  j  «  oE  exists,  provided  that  the  generating  electric 
field  is  weaker  than  the  critical  plasma  field  for  electron  heating.  For  any 
gaseous,  liquid,  or  solid  plasma,  the  electrical  conductivity  o  *  |1|/|e|  is 
given  by 


a  = 


(ne  /m)x 


(5) 


since  the  electrons  of  mass  m  <<  M  dominate  the  electric  current  transport  in 
plasmas.  The  interaction  frequency  x  ^  of  the  electrons  is  the  sum  of  the 
interaction  frequencies  for  the  v-particle  interactions, 


-1 

T 


(6) 


since  the  probabilities  (frequencies)  x^  for  many-particle  interactions  of  the 
order  v  are  additive  (v  =  2  for  binary,  v  =  3  for  ternary,  etc.).  A I  is  related 
to  the  total  number  N*  of  (charged)  particles  of  the  system  by  W  =  N*  -  1  »  1. 
For  physical  reasons,  the  conductivity  o[sec  of  a  fully  ionized, 

classical  plasma  can  depend  only  on  the  dimensional  plasma  parameters 

3/2  1/2  -1  -3  2  -2 

e[cm  gr  sec  ],  m[gr],  n[cm  ],  KT(gr  cm  sec  ],  and  the  characteristic 

dimensionless  constant  Z  =  n^/n(KT  =  thermal  energy).  The  conductivity  a  and 

the  parameters  e,  m,  n,  and  KT  have  the  dimensions  V (L  =  dimension  of  length, 

T  =  dimension  of  time,  M  =  dimension  of  mass: 


V[o]  =  r1  ,  V[e]  =  L3/2M1/2r1  ,  V[m]  =  M  ,  Z?[n]  =  L'3  ,  P[KT]  -  ML2!"2  . 


(7) 


Dimensional  theory  is  based  on  the  axioms  of  Duprd.3^’3'^  Accordingly,  the 
secondary  quantity  o  is  given  in  terms  of  the  primary  quantities  e,  m,  n,  and 
KT  by15,16 


(8) 


a  =  c/^n^KT)14** 

is  a  dimensionless  coefficient  which  depends  on  dimensionless  parameters 

such  as  Z  =  n/n^,  and  can  only  be  determined  by  means  of  a  detailed  physical 

model.  C  is  either  a  true  constant  of  order  of  magnitude  one,  C„  ~  1,  or  it 

u  z 

is  a  slowly  varying  function  (quasi-constant).  Comparison  of  the  powers  of  the 
independent  dimensions  L,  M,  and  T  [Eq.  (7)]  in  Eq.  (8)  gives  the  compatibility 
equations. 


~N1  -  3N3  +  2N4  =  0  ,  ^  +  n2  +  N4  =  0  ,  -Nx  -  2N4  =  -1  .  (9) 

These  are  three  independent  equations  (since  only  three  independent  dimensions 
L,  M,  and  T  exist)  which  determine  three  of  the  four  powers  in  terms  of  the 
fourth, 

Nx  =  1  -  2N  ,  N2  =  “  |  >  N3  =  \  ~  3**  *  n4  =  N  .  (10) 

Combining  Eqs.  (8)  and  (10)  yields  a  conductivity  expression  o  *  a„,  which 

N 

contains  a  still  undetermined  power  N, 

CTN  =  cZN<ne2/'n)1/2<e2nl/3/KT)"N  .  (11) 


In  order  to  understand  the  physical  meaning  of  Eq.  (11)  ,  it  is  rewritten  in 
the  form  of  Eq.  (5), 


2  -1 

=  ne  /mx^  ,  v  =  N/3  +  3/2  , 


(12) 


where 


"v1  =  Czi(ne2/m)1/2(e2nl/3/KT)3(V  "  3/2)  .  V  -  2,3,4, ...N  .  (13) 


It  is  now  seen  that  Eq.  (11)  or  the  equivalent  Eqs.  (12)  and  (13)  represent  the 
conductivity  of  a  hypothetical  plasma,  in  which  each  electron  experiences  only 
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many-body  interactions  of  a  fixed  order  v,  since  the  probability  for  a  v-body 
interaction  has  the  n-dependence 


c;V 

Zv 


1 


2 , 3, 4, ... N 


(14) 


For  example,  for  a  hypothetical  plasma  with  2-body  interactions  only  (so-called 
ideal  plasma),  Eqs.  (12)  and  (13)  give 


2/  "I 

o2  =  ne  / mt2 


-1  „-l  4  -1/2  ,  . -3/2 

r 2  =  Cz2e  m  n(KT) 


(15) 


where  the  dimensionless  coefficient  is  known  from  the  kinetic  theory  of  binary 

collisions, 33  C  ^  =  (3/4)  (2/tt) 3/^2/Z  iln  A,  i.e.,  is  a  quasi-constant  which 

varies  only  slightly  with  n  via  the  Coulomb  logarithm  £n  A. 

Since  the  probabilities  for  the  individual  v-body  interactions  are 

additive,  Eqs.  (12)  and  (13)  result  in  the  following  formulas  for  the  inter- 
-1  -1  2 

action  frequency  t  =  and  the  conductivity  o  =  ne  /mT  of  actual  plasmas, 

in  which  v  =  2 , 3, 4, . . . -body  interactions  take  place: 


-1 

T 


N 


l 

v=2 


c,1^ 

Zv  p 


y3(v  -  3/2) 


(16) 


2 , ,  v  3(v  -  3/2), 

o  =  ne  /[m  )  C„  w  y  ]  , 

v=2  ZV  P 


where 


~  ,  1/ 2  3(v  -  3/2)  _  ,  , 

CZv  =  (47T)  Z  CZv  *  v  =  2,3,4,., 


.W 


=  (47rne2/m)1/2  ,  y  =  Ze2n1/3/KT 


(17) 

(18) 
(19) 


are  dimensionless  coefficients,  the  plasma  frequency,  and  the  nonideality 
parameter,  respectively. 

In  order  to  expose  the  many-body  effects  (v  >  3)  of  the  nonideal  plasma 
for  comparison  with  the  corresponding  formulas  of  binary  kinetic  theory  (v  »  2) , 
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Eqs.  (16)  and  (17)  are  rewritten  as 


x-l  =  m-l/2„4_^^-3/2r,-l 


N 


e  n(KT)_‘,'Z[C~2  +  I  ^(y/Z) 

v=3 


-l^,„3(v  -  2) 


(20) 


a  =  m-1/2e-2(KT)3/2/[C~2  +  l  C~J(y/Z)3(V  '  2)] 

v=3 


(21) 


Equations  (20)  and  (21)  give  the  interaction  frequency  and  conductivity  of 
nonideal  plasmas  in  terms  of  a  series  in  y,  which  converges  rapidly  for 
0  <  y  <  1  and  converges  for  any  y  >  1,  since  it  is  finite  (1  «  N  <  ®) , 


l  C~1(v/Z) 3(v-2)  =  C~^(y/Z)  3  +  C‘J(y/Z)6  +  ^(y/Z)9  +  ., 
v=3 


CS(y/Z)3(N~2)*  (22) 


For  ideal  plasmas,  y  -*■  0,  Eqs.  (20)  and  (21)  reduce  to  x  3  =  and  o  *  , 

in  accordance  with  the  kinetic  theory  of  binary  interactions.  For  weakly 
nonideal  plasmas,  y  <<  1,  Eqs.  (20)  and  (21)  show  that  x  3  £  x23  and  a  ~ 

For  moderately  nonideal,  10  1  ~  y  £  10°,  and  strongly  nonideal,  10®  <  y  <  y, 
plasmas,  the  electron  interaction  frequency  x  3  increases,  and  the  conductivity 
decreases  considerably  —  and  by  orders  of  magnitude  —  respectively.  These 
theoretical  results  are  in  agreement  with  measurements  on  nonideal  plasmas,3  3® 
which  exhibit  considerably  smaller  conductivities  than  expected  from  binary 
collision  theory. 

It  should  be  noted  that  Eqs.  (16)  and  (17)  or  (20)  and  (21)  are  applicable 
to  nondegenerate  plasmas  only,  i.e.,  to  densities  n  <  n  or  interaction  parameters 


0  <  y  <  y  ,  y  =  21/3  - —  /h.,„  =  2.823  x  102ZT-3/2  .  (23) 

(KT/2irm)i/ 


For  a  physical  interpretation  of  the  above  results,  the  partial  collision 


frequency  in  Eq.  (13)  for  the  v-body  interaction  of  a  conduction  electron  with 
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v-1  other  charged  particles  (e,i)  is  rewritten  in  the  form 


Ty1  ~  (nQr) v  1  v/r  ,  v  =  2,3,4 ,...N  ,  (24) 

where 

n  ~  n±  ,  Z  ~  1  , 

v  ~  (KT/m)1/2 

2  (25) 

r  ~  e  /KT 

Q  ~  (e2/KT)2 

are,  as  to  order  of  magnitude  (~) ,  the  electron  or  ion  density,  the  speed  of 

the  conduction  electron  relative  to  the  interaction  partners  (e,i),  the  Coulomb 

interaction  radius,  and  the  Coulomb  scattering  cross  section,  respectively 
-2 

(Q  *  irr  ).  The  bar  designates  the  thermal  average. 

_  3 

In  Eq.  (24),  Qr  [cm  ]  is  the  interaction  volume  of  one  scattering  partner, 
w  =  nQR  [1]  is  the  spatial  probability  for  a  binary  interaction,  wv-^  = 

(nQr)V  1  [1]  is  the  spatial  probability  for  v-1  simultaneous  binary  interactions 
(in  a  At  ~  r/v)  of  the  conduction  electron  with  v-1  other  charged  particles 
(v-body  interaction),  and  v/r  [sec  is  the  frequency  of  interactions  of  the 
conduction  electron,  which  occur  at  distances  r.  Accordingly,  the  frequency  for 
a  v-body  interaction  of  a  conduction  electron  with  v-1  other  charged  particles 
is 

tv1  =  w^1  v/r  ,  w  =  nQr  ,  v  =  2,3,4,...N  .  (26) 
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APPLICATION 


For  the  practical  use  of  the  conductivity  formula  {Eq.  (21)],  the  dimen¬ 
sionless  coefficients  have  to  be  determined  either  experimentally  or  by 
physical  arguments,  since  a  complete  kinetic  equation  for  many-particle  Coulomb 
interactions  is  not  available.  For  moderately  nonideal  plasmas,  the  conductivity 
is  in  good  approximation  given  by  Eq.  (21)  as 


o  — 


(KT) 


3/2 


1/2  2,.-l  , 
m  e  [CZ2  + 


cz3(Y/z)3] 


0  <  y  <  1 


(27) 


since  the  next  higher  term  in  the  denominator  is  of  the  order  y  •  Thus,  for 

0  <  y  <  1,  it  is  sufficient  to  calculate  the  conductivity  of  nonideal  plasmas 

from  a  physical  model  of  two-  and  three- particle  Coulomb  interactions.  In  this 

case,  only  two  dimensionless  constants,  Cz2  and  C^,  have  to  be  determined. 

For  ideal  plasmas,  the  coefficient  Cz2  has  been  evaluated  by  means  of  the 

Boltzmann  equation  for  an  unshielded  Coulomb  potential  $  =  Ze/r  (Rutherford 

-1  2  1/2 

scattering  cross  section)  as  C„„  <*  Z  /In  A,  where  A  *  [1  +  b  / b  )  ]  and 

Z2  max  o 

2 

bQ  =  Ze  /3KT  (see,  e.g.,  Ref.  11).  Different  authors  prefer  either  the  mean 
ion  radius,^  b  =  (3Z/4irn)^^,  or  the  Debye  radius, ^  b  *  D,  as  upper 
impact  parameter  in  order  to  avoid  the  Coulomb  divergence  of  the  binary  colli¬ 
sion  integral.  Accordingly,  either^"*"  A  =  [1  +  9(3Z/4ir) ^y”^]^^  or^ 

-3  1/2 

A  =  [1  +  9(Z/4tt(1  +  Z))y  ]  .  For  both  choices,  in  A  +  0  for  y  >  1  and 

In  A  <<  1  for  y  <  1. 

For  a  classical  nonideal  plasma,  0  <  y  <  y,  a  physically  meaningful  Coulomb 
logarithm  is  obtained  by  evaluating  the  binary  collision  integral  for  a  shielded 
Coulomb  potential  <t>  =  Ze  exp(-r/6)/r  (Wentzel  scattering  cross  section).  This 
approach  does  not  require  an  artificial  cutoff  of  the  impact  parameter  and 
gives 


Cz2  =  (3/4)  (2/tt)1/2/Z  In  A 
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(28) 


1 
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where 

A  =  8KT/(h2/mS2)  »  1  (29) 

for  nondegenerate  plasmas.  A  is  proportional  to  the  ratio  of  thermal  (KT)  and 

2  2 

quantum  potential  (f»  /m6  )  energies  of  the  electrons,  since  the  scattering  in 
the  shielded  Coulomb  potential  is  a  wave-mechanical  process  (independent  of 
n  $  n).33  Equation  (29)  contains  not  only  the  effects  of  binary  interactions 
at  distances  r  ~  6  but  also  the  collective  many-body  interactions  at  distances 
r  >  6. 

The  electric  shielding  length  6  of  the  Coulomb  potential  of  the  classical, 
nonideal  plasma  can  depend  only  on  Z  and  the  dimensional  parameters  e,  m,  n, 
and  KT.  Dimensional  analysis  shows  that  a  6^  with  a  still  undetermined  power  N 
exists  which  is  independent  of  m, 

6N  =  CZNn“1/3(e2n1/3/KT)“N  .  (30) 


The  interaction  length  is  the  linear  superposition  6 
to  satisfy  the  limiting  conditions. 


6„  +  <5™  .  since  6  has 

N- 


6  =  [Z/ 4tt (1  +  Z)]1/2y  1/2n  1/3  ,  y  «  1  (N  =  1/2)  ,  (31) 

-1/3 

6  =  (47rn/3Z)  '  ,  y  »  1  (N  =  0)  ,  (32) 


corresponding  to  the  Debye  and  mean  ion  radii,  respectively.  Elimination  of 
the  dimensionless  constants  by  means  of  Eqs.  (31)  and  (32)  leads  to  the  shield¬ 
ing  length 


6  =  (4W3Z)"1/3[1  +  (4tt/3Z)1/3(Z/4tt(1  +  Z) ) 


(33) 


By  Eq.  (29),  A  »  1  for  0  <  y  <  l,  and  A  >  A(n)  =«  4tt ( 3Z/tt ) 2 / 3  for  1  <  y  <  y 


13 


E 


It  is  seen  that  Eqs.  (29)  -  (33)  provide  a  satisfactory  theory  of  the  Coulomb 
logarithm  of  nondegenerate  nonideal  plasmas. 

In  order  to  calculate  the  coefficient  for  triple  e-i-i  interactions 
(e-i-e  collisions  are  disregarded  consistent  with  the  disregard  of  e-e  colli¬ 
sions  in  the  evaluation  of  C^) ,  application  is  made  of  the  model  Eq.  (24) 
which  yields 


\  ~  (ni%i?ei)V"1  *ei/?ei  »  0  <  Y  <  1  ,  (34) 

where 

n1  =  n/Z 

vei  =  (4/3)(8KT/irm)1/2 

2  (35) 

r  =  Ze  /3KT 

Qei  =  (Ti/4)(Ze2/KT)2  In  A 

are  the  exact  thermal  averages  known  from  kinetic  theory. Comparison  of  Eq. 
(34)  with  Eq.  (13)  shows  that 


CZv  s  4(8/tt)1/2(jj  Z2£n  A)V_1Z_1  ,  0  <  y  <  .  (36) 

-1  1/2 

Accordingly,  C  —  (4/3)(ir/2)  Zin  A  for  v  =  2  in  agreement  with  Eq.  (28),  and 
for  v  =  3 

-1  i  /?  y  9 

CZ3  =  (»/9)(ir/2)A,V(ln  A)Z  ,  0  <  y  <  1  .  (37) 


Substitution  of  Eqs.  (28)  and  (37)  into  Eq.  (21)  yields  for  the  electrical 
conductivity  of  moderately  nonideal  plasmas: 


a  — 


(3/4)  (2/tt)1/2(KT)3/2 


m1/2Ze2[in  A  + 


(7T/12)Z~1(X.n  A) 2y3 ] 


0  <  y  <  1 


(38) 
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The  Coulomb  logarithm,  in  A,  is  evaluated  in  Eq.  (29)  in  dependence  of  n,  T,  and 
Y- 

Equations  (34),  (36),  and  (37)  are  based  on  binary  e-i  and  triple  e-i-i 
collisions  and  collective  many-particle  interactions  which  are  considered 
through  the  shielded  Coulomb  potential  with  the  interaction  length  6  =  6(n,y)» 

Eq.  (33).  The  approximately  equal  signs  in  these  equations  are  a  reference  to 
the  disregard  of  binary  (e-e)  and  triple  (e-i-e)  collisions.  Furthermore ,  we 
have  restricted  the  applicability  of  the  results  to  moderately  nonideal 
(0  <  y  <  1)  plasmas,  since  the  Coulomb  logarithm  has  been  calculated  without 
considering  the  influence  of  triple  collisions.  The  latter  effect  is  in  all 
probability  not  quantitatively  significant  since  £n  A  is  a  slowly  varying 
function  of  A.  A  more  accurate  determination  of  the  dimensionless  coefficients 
C^  and  Cz^  has  to  be  postponed  until  a  kinetic  equation  for  nonideal  plasmas 
is  available,  which  takes  into  account  not  only  binary  but  at  least  also  triple 
interactions  and  correlations. 

The  experimental  data  for  nonideal  alkali  and  noble  gas  plasmas^ 

(0.1  <  y  <  1)  indicate  that  the  electrical  conductivity  is  roughly  an  order 

of  magnitude  smaller  than  predicted  by  the  theories  of  ideal^  and  weakly 
,  ,  12 

nonideal  plasmas.  In  Fig.  1,  isobars  of  the  dimensionless  Coulomb  conductivity 
1/2  2  3/2 

a*  =  m  e  a/(KT)  are  reproduced  versus  the  number  of  electrons  in  the  Debye 
3  -3/2 

sphere,  =  4irD  n/3  ~  y  [Eq.  (3)],  showing  (1)  conductivities  according  to 

the  ideal  plasma  theory^  (b  =  D) ,  (2)  computer  conductivities  from  molecular 

max 

dynamics  and  Monte  Carlo  methods  (with  error  estimates)  by  Valuev  and  Norman, ^ 

and  (3)  experimental  conductivities  for  a  cesium  plasma  at  a  pressure 
4  10 

p  =  5  x  10  Pa  by  Dikhter  et^jil.  Conductivity  curves  based  on  the  present 

analytical  theory  [Eq.  (38)]  are  shown  for  (4)  lithium  and  (5)  cesium  plasmas 

4 

at  a  pressure  p  =  5  x  10  Pa  for  comparison. 
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Figure  1  demonstrates  that  the  theoretical  conductivity  values  from 
Eq.  (38)  are  correct  as  to  order  of  magnitude  and  lie  well  within  the  errors 
of  the  experimental  data^’^  and  the  computer  experiments.^  Equation  (38) 
predicts  a  slight  decrease  of  a*  for  NQ  <  0.1,  which  is  due  to  the  contribu¬ 
tions  from  the  triple  interactions  [Fig.  1,  curves  (4)  and  (5)].  This  effect 
was  observed  by  Kulik  et^  al. ^  in  measurements  on  cesium  plasmas. 
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III.  CONDUCTIVITY  OF  NONIDEAL  CLASSICAL  AND  QUANTUM  PLASMAS 


By 

H.  E.  Wilhelm 


ABSTRACT 


The  electrical  conductivity  of  fully  ionized,  moderately  nonideal 

2  1/3 

plasma  with  Coulomb  interaction  parameters  0.1  <  y  <  1,  where  y  =  Ze  n  /KT 
is  the  ratio  of  Coulomb  and  thermal  energies,  is  calculated  for  displaced 
Maxwell  and  Fermi  electron  distributions,  respectively.  The  electrons  are 
scattered  by  an  effective  Coulomb  potential  <J>  (r)  =  Zer  ^exp(-r/6),  which 
considers  binary  (0  <  r  <  6)  and  many-body  (6  <  r  <  “)  interactions.  The 
shielding  distance  is  given  by  <5  *  a(4irn/3Z)  with  a  -  aQy  N-1  for  classi¬ 
cal  plasmas  and  6  =  8(4im/3Z)  with  8  =  8QY  NT  M~1  for  quantum  plasmas, 

2  1/3  2  -1  2/3 

where  r  -  Ze  n  /fi  m  n  is  the  ratio  of  Coulomb  interaction  and  quantum 
potential  energies  of  the  electrons.  It  is  shown  that  the  resulting  conductivity 
formulas  are  applicable  to  densities  up  to  four  orders  of  magnitude  higher 
than  those  of  the  ideal  conductivity  theory,  which  breaks  down  at  higher 
densities  because  the  Debye  radius  loses  its  physical  meaning  as  a  shielding 
length  and  upper  impact  parameter. 
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INTRODUCTION 


1-3) 

The  theory  of  the  electrical  conductivity  of  fully  ionized  plasmas  based 

on  the  Boltzmann  equation,  the  Fokker-Planck  equation  (derived  by  expanding  the 

binary  collision  integral  for  the  small,  successive  velocity  changes  of 

Coulomb  scattering) ,  or  the  Lenard-Balescu  equation  (taking  into  account  the 

dielectric  properties  of  the  medium)  is  in  agreement  with  the  experimental  data 

for  rarefied  high-temperature  plasmas,  y<<:1*  The  interaction  parameter  is 

2  1/3 

defined  as  the  ratio  of  (average)  Coulomb  interaction  (Ze  n  )  and  thermal  (KT) 
energies  (n  is  the  electron  density  and  Z  the  ion  charge  number) , 
y  «  Ze2n1/3/KT  =  1.670  x  10~3Zn1/3/T 

in  cgs-units  which  will  be  used  throughout.  The  conventional  transport  calcu- 

1-3)  3/2  1/2  2 

lations  give  an  electrical  conductivity  of  the  form  a  ~  (KT)  /m  e 

2  1/2 

for  classical  ideal  plasmas,  where  [1  +  (D/p0)  ]  s  D/pQ  ^or  D>>po‘  D  is 

the  maximum  impact  parameter  for  which  the  Deby*  length  is  used,  and  pq 

2 

is  the  average  impact  parameter  for  90°  deflections  (Landau  length),  pQ  =  Ze  /2KT. 

The  condition or  is  satisfied  only  for  not  too  low  temperatures 

4) 

T  and  not  too  high  densities  n.  Conductivity  formulas  with  this  Coulomb 
logarithm  break  down  for  large  interaction  parameters  y  and  densities  n,  since 
the  Debye  radius 


D  =  [ Z/4tt (1  +  Z)  ]1/2Y~1/2n_1/3 


-8 


becomes  smaller  than  the  atomic  dimension  10  cm  and,  thus,  completely  loses 
its  physical  meaning  as  an  electric  shielding  length  and  maximum  impact  parameter. 
E.g.,  for  T  =  10^  °K,  Y>10*3and  D<10  ^cm  if  n>102*3cm  3.  Moderately  nonideal 
plasmas  with  y~1  are  readily  generated  through  shock  wave  compression  and  exhibit 
conductivities  of  the  order  a  ~  10^  -  102mho/cm3  which  are  much  smaller  than 
those  which  would  be  obtained  by  applying  the  conductivity  formula  for  ideal 
plasmas  in  the  nonideal  regime. 
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Although  there  are  some  bulk  measurements  of  the  electrical  conductivity 

5-8) 

of  nonideal  cesium  and  noble  gas  plasmas  available  ' ,  theoretical  explanations 

of  these  results  are  still  missing.  The  momentum  and  energy  transport  in 

9) 

weakly  nonideal  plasmas, y  <<1.  was  treated  by  Wilhelm  by  means  of  an 
exponentially  shielded  Coulomb  potential,  which  permits  to  consider  not  only 
short-range  binary  (r  5  D)  but  also  long-range  many-body  (r  >D)  interactions.  This 
interaction  model  was  used  shortly  afterwards  by  Rogov^^  for  the  calculation 
of  the  conductivity  of  weakly,  nonideal  argon  and  xenon  plasmas  with  Debye 
shielding. 

For  moderately  nonideal  plasmas ,  0. 1  <  y  S  1,  various  phenomenological  approaches 

have  been  used  to  extend  the  conductivity  formula  of  ideal  plasmas,  e.g.,  Goldbach 

et  al  multiply  the  Debye  length  D  with  a  free  parameter  x(p)which  is  chosen 

to  match  the  experimental  data,  i.e.  to  compensate  for  the  too  rapid  decrease  of 

D  with  pressure.  A  kinetic  equation  has  been  proposed  for  nonideal  plasmas 

12) 

by  Klimontovich  ,  which  considers  spatial  correlations  and  temporal  retardation 
in  the  collision  integrals.  This  equation  appears  to  have  not  yet  lead  to  transport 
coefficients  because  of  the  mathematical  difficulties  associated  with  its  solution. 

In  the  following,  the  momentum  relaxation  time  and  the  electrical  conductivity 
of  (i)  classical  and  (ii)  quantum  plasmas  is  calculated  for  intermediate  non¬ 


ideal  conditions,  0.1  <  y  <,  1.  For  this  region  of  interaction,  the  concept  of 

Debye  shielding  already  breaks  down  since  the  number  of  particles  in  the  Debye 
3 

sphere  4ttD  / 3  is  no  longer  large  compared  with  one  for  y  >  0.1.  This  difficulty 

13) 

can  not  be  remidied  by  replacing  D  with  the  quantum  mechanical  shielding  length 


separates  the  region  in  which  an  electron 
21 


r 


experiences  few-body  encounters  (r  <  6)  from  the  region  in  which  an  electron 

experiences  many-body  interactions  (r  >6)  in  a  nonideal  plasma,  as  long  as 
—8  24  —3 

6  >  10  cm  (n^  <10  cm  ) .  Thus,  the  mean  ion  distance  evolves  naturally 
as  the  characteristic  interaction  distance  for  nonideal  plasmas,  for  which 
Debye  and  Fermi  shielding  fail.  We  demonstrate  mathematically  that  6  *  y  Nn^ 
with  0  <  N  <  \  for  classical  plasmas  and  6  =  y  Nr  Mn^  with  0  <  H  £  h  for 
quantum  plasmas  (T  =  Ze^n^^/if^m  i.e.  y  ^  ~  1  and  r  ^'1  are  correction 

factors  which  are  insignificant  since  the  plasma  conductivity  depends  logarith¬ 
mically  on  6. 

We  calculate  the  electrical  conductivity  of  plasmas  with  (i)  Maxwell  and 

(ii)  Fermi  distributions  of  the  electrons,  when  all  ions  have  the  same  charge 

number  Z.  The  electrons  are  assumed  to  be  scattered  by  the  exponentially 

shielded  Coulomb  potential  <fi  *  Zer  1exp(-r/6)  which  takes  many-body  interactions 

at  distances  r  >  6  into  account.  The  considerations  are  applicable  only  to 

24  -3 

moderately  nonideal  conditions,  0.1  <  y  <  1,  up  to  densities  n  <<  10  cm 
Thus,  the  following  theory  is  limited  to  densities  n  well  below  the  electron 
density  in  (solid)  metals. 
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PHYSICAL  FOUNDATIONS 


The  electrical  conductivity  a  of  any  gaseous,  liquid,  or  solid  medium, 
in  which  the  electrical  current  transport  is  due  to  electrons,  is  proportional 
to  the  electron  density  n  and  the  relaxation  time  t  of  the  average  momentum 
<mv  >  of  the  electrons  (m  is  the  electron  mass  and  e  >  0  is  the  elementary 

P 

charge) 

o  =  (ne2/m)r. 

The  relaxation  time  t  is  determined  by  the  scattering  potential  and  the 
(classical  or  quantum  statistical)  kinetics  of  the  electron  gas  in  the  electric 
field. 


In  nonideal  plasmas,  the  region  0  <  r  s  8  of  binary  and  few  -  body 

collisions  and  the  region  6  <  r  <  °°  of  many  -  body  interactions  are  separated 

by  the  electric  shielding  radius  6.  Dimensional  theory  gives  for  classical 

(n  <<  n)  and  quantum  (n  £  n)  plasmas  (see  Appendix)  : 

1/1 

6  =  a(3Z/47rn)  ,  n  <<  n  , 

-N 

a  =  aQY  ,  0  <  N  <  4  , 


and 


5  =  8(3Z/4irn) 

■  •  e0r-Vs 


1/3 


n  >  n  , 

0  <  M,N  <  h  , 


aQ  =  ZN(4tt/3Z)  1/3/ [4tt(H-Z)  =  ZMaQ, 

n  =  2(2TrmKT/h2)3/2 


The  nonideality  parameters  of  the  classical  and  completely  degenerate  plasmas 
are  defined  by 


■1— 


i 


■  *+ 


The  classical  formulas  (2)-(3)  are  the  special  case  M  ■  0  (fj  ■  o.Q)  of 

4) 

the  general  quantum-mechanical  Eqs.  (4)-(5).  Eq.(2)  becomes  the  Debye  radius 

of  the  weakly  nonideal  (y<<l)  classical  plasma  for  N  =  6(N“4)“D,  whereas 

13) 

Eq. (4)  becomes  the  Fermi  radius  of  the  completely  degenerate  (n>>n)  plasma 

for  M  =  \  and  N  =  0,  6^=4$,  N=0)'D  .  For  M=N=0,  Eqs.  (2)  and  (4)  reduce  to  the 

F 

shielding  radius  of  the  strongly  nonideal  plasma,  in  which  the  kinetic  energy 

2  -1  2/3  l /3 

(KT  or  Ti  m  n  )  of  the  electrons  is  negligible,  6 (M=N=0) ~ 0Z/4tm )  ,  which  is 

the  mean  ion  distance  up  to  a  factor  aQ=  B0  ~  1.  For  these  reasons,  the  powers 

M  and  N  in  Eqs. (2) -(6)  are  limited  to  the  interval  0  <  M, N  <  %.  For  non-solid 
24  -3 

(n  <<  10  cm  )  plasmas  of  intermediate  nonideality,  0.1  <  y  <  1,  which  implies 
2  2-1  2/3 

1  <  F  <  10  since  r  =  (KT/1i  m  n  )y,  extremely  simple  relations  hold  as  to  order- 
of-magnitude: 

ct ( y)  '  1,  6(y,r)  -  1,  6  -  (3Z/4™)1/3.  (9) 

Based  on  the  above  considerations,  the  scattering  of  electrons  by  Z-times 
charged  ions  in  plasmas  of  intermediate  nonideality  is  described  by  the  shielded 
Coulomb  potential 

4> (r)  =  Zer  1  exp(-r/6),  0  <  r  <  °°  ,  0.1  <  y  <  1.  (10) 

which  contains  the  binary  and  few  -  body  collisions  at  distances  0  <  r  <  6  and  the 

many  -  body  interactions  at  distances  6  <  r  <  00  .  A  similar  Coulomb  potential  is 

used  in  the  conductivity  theory  of  metals,  although  the  use  of  such  a  "binary 

22  -3 

quasi-potential"  is  questionable  for  densities  n  >  10  cm 

^  JU 

The  differential  cross  section  o(9,g)  for  the  scattering  (g  ■+•  g  )  of 

14) 

electrons  by  the  potential  (10)  is  in  the  center  of  mass  system 

a(6, g)  =  (Ze2/2m)2/[g2sin2(0/2)  +  u2]2  ,  u  =  1i/2m6  (11) 

where  0  =  Mg,g  ),  g  =  vg  -  vi#  g  -  v*  -  v^*  ,  and 


24 


u 


the  electron  and  ion  velocities  before  and  after  the  interaction  are  designated 

by  v  .  and  v*  ,  respectively.  The  speed  u  corresponds  to  a  de  Broglie 
0*1  0*1 

wave  length  of  the  order  X  -  6.  For  u  +  0  or  {  +  «,  Eq.  (11)  reduces  to  the 

14) 

Rutherford  cross  section. 

The  scattering  cross  section  a(0,g)  is  strictly  valid  only  in  the  Bom 

14  ) 

approximation  .  Contrary  to  what  one  might  expect  in  general  for  the  latter, 

Eq.(ll)  describes  in  good  approximation  the  scattering  in  the  exponentially 

decaying  potential  (10)  because  of  the  peculiarity  of  the  Coulomb  interaction. 

The  Coulomb  interaction  <j>  -  1/r  has  the  unique  property  that  the  Born  approximation 

and  the  exact  wave  mechanical  approach  give  the  same  scattering  cross  section'*^ 

(identical  with  the  Rutherford  formula).  In  the  region  0  <  r  <  6,  the  interaction 

potential(lO)  is  practically  Coulombic,  and  thus  the  Born  approximation  gives 

the  correct  solution.  In  the  region  6  <  r  <  °°,  the  interaction  potential  (10) 

is  effectively  screened,  i.e.,  the  Born  approximation  gives  the  correct  solution 

because  <J>(r)  is  small.  In  the  transition  zone  r  a  6,  the  Born  approximation 

holds  fairly  well  for  reasons  of  continuity. 

The  relaxation  time  t  is  obtained  by  evaluation  of  the  collision  integrals 

for  the  electron  momentum  mvg  for  the  (i)  classical  and  (ii)  degenerate  plasma, 

respectively.  Both  in  the  cases  of  classical  and  Fermi  statistics,  the  particle 

velocities  v  and  v*  before  and  after  the  interaction  are  interrelated  by  the 
e  y  i  e  9 1 

classical  conservation  equations  for  momentum  and  energy. 
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CONDUCTIVITY  OF  CLASSICAL  PLASMA 

According  to  kinetic  theory,  the  average  momentum  density  nm(<vg>  -  <v^>) 
exchanged  per  unit  time  between  electrons  and  ions,  interacting  with  the  Coulomb 
potential  (10),  is  given  by  the  collision  integral  for  mv^,  which  determines 
the  momentum  relaxation  time  t  , 


-nm(<v  >  -  <v.>)/x  = 

e  1 

m  /•••/  ve[fe(v*)  f^v*)  -  fe(vg)  fi(vi)]g  0(6, g)  dO  dvfi  d^  .  (12) 

The  scattering  cross  section  a(6,g)  is  given  in  Eqj(ll)  and  the  solid  angle  element 
is  dO  =  sin  0  d6  d4>.  In  response  to  an  applied  electric  field  £,  the  electrons 


and  ions  drift  with  velocities  <v  >  and  <v.>  so  that  their  distribution  functions 

e  i 

are  displaced  Maxwellians, 

fs(vs)  =  ns(mg/27rKTs)3/2  exp[--|mg(vg  -  <vg>)2/KTg],  s*e,i.  (13) 

Eq^l3)  represents  a  5-moment-approximation  to  the  nonequilibrium  solution  of  the 

Boltzmann  equation.  The  perturbations  of  f  (v  )  due  to  viscous  stresses  and 

s  s 

heat  flows  are  neglected  in  Eq.(13),  since  they  yield  only  corrections  of 


higher  order  to  the  conductivity. 

9) 

The  collision  integral  (12)  is  integrated  by  standard  methods  for  subsonic 

drift  velocities,  |  <vg>  -  <v^>  |  <  (2KT/m)3^2,  with  the  usual  approximations 

(m  =  m  m  /(m  +  m  )  s  m  =  m,  T  =  m  [ (T  /m  )  +  (T. /m, )]  s  T  =  T) . 
ei  e  l  e  i  e  es  es  e  e  i  i  e 

(For  supersonic  drift  velocities,  a  linear  response  1  =  oE  between  current  density 

-±  *  Q) 

j  and  electric  field  E  does  no  longer  exist.  )  The  resulting  relaxation  time  is 


given  by: 


•T1  =  |(2KT/rm)1/2  n  Q 


Q  =  f(Ze2/  KT)2  L 
L  =  eA  E1(A_1) 

where 

A  =  2KT/mu2=a+2(8m/^)(4W3Z)"2/3KT 


(15) 

(16) 

(17) 
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by  Eqs.(2)  and  (11)  for  the  classical  plasma.  Furthermore 


El(x) 


T*  -  In  x  -  Z  (-l)m  xm/m(ml) 
m=l 


(18) 


is  the  exponential  integral  of  order  one  (r*  =  0.477..  =  Euler's  constant). 

1  ON 

The  latter  satisfies  the  inequalities,  for  x  >  0,  ' 


18) 


*sZn(l  +  2/x)  <  eXE1(x)  <  Zn( 1  +  l/x),(l  +  x)  ^  <  exE^(x)  <  x 
Accordingly,  Eq.  (17)  gives  formally  for  small  and  large  arguments  x  =  A  \ 

L  -  In  A,  A  >>  1  ;  L  =  A  ,  A  <<  1 

Rewriting  A  in  terms  of  the  thermal  and  quantum  potential  energies  shows  that 
for  classical  plasmas 

A  =  8Et/Eq  »  1,  ET  =  KT,  Eq  =-h2/m62. 


-1 


(19) 


(20) 


(21) 


Combining  of  Eqs.  (14)-(17)  with  Eq.(l)  yields  the  desired  electric 
conductivity  of  the  classical  plasma  of  intermediate  nonideality,  0.1  <  y  $  1: 

a  =  3(KT)3/2/2(2irm)}4  e2  ZL  (22) 

where 

-2  +2  -2/1 

L  =  In  [8mtT  aT  (4irn/3Z)  ;  KT],  A  »  1,  (23) 

by  Eqs. (17)  and  (20), 


The  conductivity  formula  (22)  differs  from  the  conductivity  of  the  ideal 
1-3) 

plasma  mainly  through  the  term  L.  The  latter  has  the  form  of  a  Coulomb 
logarithm,  L  =  tn  A  for  A  >>  1,  i.e.  for  all  densities  n  and  temperatures  T  for 
which  the  plasma  is  nondegenerate,  >  E^,  Eq.(21).  Numerically, 

A  =3.482  x  1011a+2(n/Z)"2/3T  .  (24) 

2  i  o\ 

The  corresponding  argument  A =  2KT  D/Ze  of  the  ideal  Coulomb  logarithm  ' 

In  Ap,  is 


Ap  =  1.464  x  10^  Z_1(l  +  Z)"5  n”^2  T3/^2 


(25) 
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Table  I  compares  A  of  the  nonideal  plasma  and  A^  of  the  ideal  plasma  for  large 
densities  n  and  the  typical  temperature  T  *  104  °K.  It  is  seen  that  in  A^ 


18  —3 

the  ideal  plasma  is  unacceptably  small  for  densities  n  >  10  cm  ,  whereas 

22  —3 

in  A  of  the  nonideal  plasma  has  reasonable  values  up  to  densities  n  <  10  cm  , 
4 

if  T  =  10  °K.  The  conductivity  formula  (22)  holds,  therefore,  for  densities  n 

up  to  4  orders  of  magnitude  higher  than  the  conductivity  formula  of  the  ideal 
plasma.  Eq. (22)  is  not  applicable  to  n  -  T  regions  for  which  A  <<  1,  i.e. 

Et  «  Eq,  which  imply  degenerate  electron*. 


TABLE  I:  A  and 


Ap  versus  n  for  T  =  10^  °K ,  t  =  1,  and  a  -  1. 


n[ctn  3] 


10 


18 


10 


20 


10 


22 


10 


24 


A 

"d 


3.482  x  10' 


1.035  x  10 


1.616  x  10 


1.035  x  10 


0.750  x  10 


1.035  x  10 


-1 


0.348  x  lO1, 


1.035  x  10 


The  conductivity  formula  (22)  becomes  in  cgs  -  units  or  practical  units 
(9  x  1011  sec  1  =  1  mho  cm-1) , 

o  =  1.394  x  108  T3/2  /Z  In  A  [sec'1]  =  1.549  x  lO-4  T3/2/Z  In  A  [mho  cm"1]  (26) 

where  A  is  given  in  Eq.(23).  Accordingly,  if  T  =  104  °K  and  Z  =  1,  o  =*  1.899  x  101 

—1  t  g  _3  I  _ ■»  o 0  “3 

mho  cm  for  n  =  10 _L  cm  and  o  =  3.046  x  10'L  mho  cm  for  n  =  10  cm 

1/2  2  3/2 

In  Fig.  1,  isobars  of  the  dimensionless  conductivity  a*  =  m  e  o/(KT) 
are  reporduced  versus  the  number  of  electrons  in  the  Debye  sphere, 

3  -3/2 

Np  =  4ttD  n/3  ~Y  ,  showing  (1)  conductivities  according  to  the  ideal  plasma 
4 

theory  ,  (2)  computer  conductivities  from  molecular  dynamics  methods  (with  er¬ 
ror  estimates)  15  and  experimental  conductivities  16  for  (3)  Cs  and  (4)  Li  plasmas 
4 

at  p  =  5  x  1C  Pa.  Conductivity  curves  based  on  Eq.  (22)  are  shown  for  (5)  Li 

4 

and  (6)  Cs  plasmas  at  p  =  5  x  10  Pa  for  comparison.  It  is  seen  that  Eq.  (22) 

is  in  good  agreement  with  the  machine  calculations  13  and  in  reasonable  accord 

16 

with  the  experimental  data 
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CONDUCTIVITY  OF  QUANTUM  PLASMA 


The  electrons  in  a  plasma  become  degenerate  if  their  thermal  DeBroglie 
wave  length  is  larger  than  the  mean  electron  distance,  i.e.  at  densities 
n  >  4.828  x  1015  T 3/2  . 

lx  9 1  —3 

E.g.,  for  T  =  10q  °K,  degeneracy  requires  n  >  5  x  10^  cm-  .  In  view  of  their 

large  mass  m^  >>  m,  the  ions  can  be  treated  as  classical.  The  momentum  relaxation 

time  t  of  the  degenerate  electron  gas  is  determined  by  the  quantum  statistical 

19) 

collision  integral  for  mve. 


-  nm(<v  >  -  <v  >)/t  = 
e  i 

m  /•••/  v  {f  (v*)f  (v*)[l 
6  6  e  i  l 


f  (V  ) 
ft  e 


f  i  (v±)  [i 


x  g  o(0, g)  dQ  dvfi  dv^  (27) 

where  the  scattering  cross  section  o(6,g)  between  electrons  and  ions  is  given  by 
Eq.  (11).  The  solutions  to  the  velocity  distributions  are  the  displaced 
Maxwellian  (13) for  the  ions  (s  =  i)  and  the  5-moment  Fermi  approximation  for  the 
electrons. 


f  (v  )  =  2(m/h)3  {1  +  exp[4n(v  -  <v  >)2  -  y]/KT  }_1  . 

“  “  t  6  G 

The  chemical  potential  y  =  y(n,  T)  is  determined  by  the  integral  functional 

n  =  Jfe(V  u)  d\  ' 

Again,  a  linear  response  j  =  aE  exists  for  small  drift  velocities  <v  > 

e,i 

or  weak  electric  fields  E.  Integration  of  Eq.(27)  yields,  after  standard 
approximations ,  for  the  relaxation  time  of  the  degenerate  electron  gas: 

r1  - 1 A 1/2  AV«i 

3  (KT)3/2R(n,T) 


(28) 


(29) 


where 


L 


Q 


), 


cao) 


KT  <2(n,T) 


(31) 
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and 


£(n,T)  -  |(1  +  2n  "  +  ...), 
z  n 

R(n.T)  =  (1  +  2'2  "  +  .. 


3,3ir  '  ns  2/3, ,  ,  5tt  ,3ir  ' 
Q(n,T)  =  5  (—4—  ")  U  +  "i2  (— 4- 

,  1/2  n  .  o_2  »  1/2 

R<n.T)  =  f/tl  +  ^(^Y- 


n  <  n(T) , 

(32) 

) ,  n  <  n(T) , 

(33) 

2  -4/3 

-  ")  +...], 
n 

n  >  n(T) , 

(34) 

n  -4/3 

n>  n(T). 

(35) 

Equations  (32)  -  (33)  and  Eqs.  (34)  -  (35)  result  from  expansions  of  the  Fermi 
distribution  (28)  in  the  collision  integral  (27)  for  densities  n  <  n(T)  and 
n  >  n(T),  respectively.  Eqs. (32)  and  (34)  indicate  that  A^  »  1  for  n  <<  n 

and  «  1  for  n  »  h  since  A^  =  (2/3)(g/ot)2  A(£(n,T)  by  Eq.  (31). 


The  series  are  based  on  expansions  of  the  normalization  integral  id  Eq.(28), 
which  gives  the  chemical  potential  p  explicitly  as  a  function  of  n  and  T, 

=  ln{~  (1  +  2~3/2£)1  +  (|  -  3~3/2)(")2  +  ...]},  n  <  n(T),  (36) 

KT  n  n  n 

v  /3ff 1/2  n.  2/3  u2  3tt1/2  n  -4/3  n4  3tt1/2  n.-8/3  .  ~m  frn 

KT  ~  (  4  ^  [1_12(— 80  ^  4  n  +”<]*  n>n(T)(37) 

Combining  of  the  conductivity  formula  in  Eq. (1)  with  the  relaxation  time  of 

Eq.(29)  yields  for  the  electrical  conductivity  of  the  degenerate  electron  plasma 

of  intermediate  nonideality,  0.1  <  T  S  1: 

3(KT) 3/2  R(n,1) 

0 - 1/2  2  (38) 

8(2m)  7  e  Z 

where  is  given  by  Eq.(30).  In  the  limiting  cases  of  large  and  small  values  of 


Lq  "  *«tj(f)2A£(n,T)], 

LQ  -  -3  (|)  V(n,T) 


Aq>>  1, 

Aq<<  1. 
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since  / l  =  (2/3)  (8/a)2/lO  by  comparison  of  Eqs.(31)  and  (17). 

For  n /n  ->  0,  Eqs.(38)  and  (39)  reduce  to  the  classical  conductivity, 
Eqs .  (22)  and  (23),  since  ^(n.T)  -*■  4//F  and  Q(n,T)  -*■  3/2  for  n/fl  -*■  0  by 
Eqs .  ( 32)  -  (33).  On  the  other  hand,  Eqs.(38)  and  (40)  give  in  the  limit  of 
complete  degeneracy,  n/fl  -*• 


3 

9h  n 


„9  2  2 

2  M  e  Z  L. 


(41) 


where 


.  12  ,  -  2  2  ,.,1/3  2 

A()  =  -r-  (3tt  t  /2)  8 


(42) 


by  Eqs.  (34)  and  (35).  Since  'i  =  h  and  N  =  0  for  n/n  -*•  “>, 

2  2/3  -1 

8  =  (Z/4it(1+Z)  ](4tt/3Z)  '  T 

by  Eq. (4) ,  i.e.,  8^1  denending  on  the  magnitude  of  T  =  Ec/E^,  Eq.(8). 

Equation  (43)  agrees  with  the  expression  for  the  conductivity  of  a  low 
,  20) 

temnerature  metal.— 


(43) 
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VWPpIl 


IP1 


GENERALIZATION 


Nonideal  plasmas  may  exhibit  not  only  a  high  degree  of  single 

ionization  but  also  multiple  ionization,  due  to  lowering  of  the  ionization 

energies  by  the  internal  Coulomb  fields,  and  overlapping  of  the  atomic  wave 

21) 

functions  at  sufficiently  high  pressures.  In  an  electrically  neutral  plasma 
with  N  species  of  ions  (i)  of  charge  Z^e  and  density  n^,  the  electron  density 

n  and  entire  ion  density  n(i)  are  related  by 

N  N  .... 

n  =  E  Z  n  ,  n(i)  =  En  .  (44) 

i=l  i=l 

Since  the  probabilities  for  interaction  between  the  electrons  and  ions 
of  type  i  =  i,  2,...N  are  additive,  the  momentum  relaxation  time  of  the 
electrons  is  in  presence  of  N  ion  components  (i)  given  by 


-1 

T 


N 
=  I 
i=l 


-1 


(45) 


From  Eqs. (44)  to  (45)  follows  that  the  derived  conductivity  formulas  are 


generalized  to  many- ion-component  plasmas  by  means  of  the  substitutions: 
-  N 

Z  -*  1  =  E  n  Z./n(i)  , 
i=l 

— N 


Z2  Z2  =  l  n  Z2/n(i)  , 
i=l  1  1 


Z"1  =  z/z2 


(46) 

(47) 

(48) 


Since  it  is  extremely  difficult  to  calculate  accurately  the  ion  densities  n^ 

22) 

in  many-component  nonideal  plasmas,  it  is  advisable  to  make  use  of  equivalent 
approximations  [see,  e.g.,  Eq . ( 9) ]  to  avoid  too  cumbersome  conductivity  expressions 
in  practical  applications. 


APPENDIX 


Since  adequate  methods  for  the  solution  of  the  many-body  problem  are 

not  available,  we  derive  the  shielding  radius  <$  of  nonideal  plasmas  from 

23) 

first  principles  of  dimensional  theory  .  This  method  gives  6  in  dependence 
of  the  relevant  dimensional  plasma  parameters  (e,  m,  n,  KT,  fi) ,  but  leaves  a 
dimensionless  proportionality  constant  of  order  one  undetermined,  Cq  -  1. 

The  latter  can,  however  be  found  from  the  physical  argument 

2  U 

6  =  D  =  [KT/4ir(l  +  Z)e  n]^  for  y  <<  1*  n  <<  n*  (Al) 

since  5  approaches  the  Debye  radius  D  in  the  limit  of  the  weakly  nonideal 
classical  plasma.  The  dimensions  of  the  characteristic  plasma  quantities, 
which  determine  6  of  dimension  L  are  given  in  terms  of  the  fundamental  dimen¬ 
sions  of  length  (L) ,  time  (T) ,  and  mass  (M) : 

D[e]  =  i3/2  M"2  T~l,  D  [m]  =  M,D[  n]  =  l“3,  D[KT]  =  ML2T~]  D[K]  ■=  WlV1.  (A2). 
A.  Classical  Plasma.  Since  in  a  classical  plasma  5  depends  on  the 

dimensional  parameters  e,  m,  n,  and  KT,  Dupre's  fundamental  theorem  of 

23) 

dimensional  analysis  demonstrates  that 

6  =  C0eNl  mNz  nN3  (KT)Nl*  (A3) 

where 


*1^  -  3N^  +  2N^  =  1,  -|n^  +  N2  +  =  0,  -N^  -  2N^  =  0,  (A4) 

by  comparison  of  the  powers  N  of  L,  M ,  and  T  in  Eq.  (A3).  Elimination  of 

,  N2  in  terms  of  N^  =  N  reduces  Eq.  (A3)  to 

6  =  Co(e2n1/3/KT)"N  n_1/3  .  (A5) 

-1/3 

It  is  seen  that  5  =  CQn  for  N  =  0  (strongly  nonideal  plasma,  y  »  1,  with 
negligible  thermal  energy)  and  6  =  D  for  N  =  1/2  (weakly  nonideal  plasma,  y  «  1) . 
Accordingly, 

Co  =  [477(1  +  z)]~1/2  (A6) 
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by  Ec .  (Al),  Eqs.(A5)  and  (A6)  are  combined  in  an  illustrative  form  for  classical 


plasmas : 

5  ‘  • 0 : » : 1/2  •  " 
This  is  Eq.  (2)  where  aQ(Z)  is  defined  in  Eq.  (6). 

B.  Quantum  Plasma.  Since  in  a  quantum  plasma  6  depends  on  the  dimensional 

23) 

parameters  e,  m  ,  n,  KT,  and  IT,  Dupre's  theorem  gives 

6  =  C  eNl  mN2  nN3(KT)N4  1^5  (/ 

o 

where 


|nx  -  3N3  +  2N4  +  2N5  =  1,  ^  +  N2  +  N4  +  N5  =  0,  -NL  -  2N4  -  N5  =  0  (A9) 

by  comparison  of  the  powers  of  L,  M,  and  T,  in  Eq.  (A8) .  Elimination  of  N  , 

N2>  in  terms  of  N4  =  N  and  =  2M  reduces  Eq.  (A8)  to 
.1/3  21/3  „ 


,  1/3  „  2  1/3  M 

6  ‘  Co (£“  KT  )  (“2“-"l  2/30  n 

"fi  m  n 


(A10) 


It  should  be  noted  that  Eqs.  (A3)  and  (A7)  are  the  special  case  =  0  of  Eq. 

(A8)  and  M  =  0  of  Eq.  (A10) ,  respectively.  For  M  =  1/2,  Eq.  (AlO)  reduces  to 

13) 

the  Fermi  shielding  length  of  the  completely  degenerate  plasma  (n  >>  n) . 


Eqs.  (A6)  and  (AlO)  are  rewritten  in  an  illustrative  form  for  quantum  plasmas: 

.  ,Ze2n1/3  -N  Ze2n1/3  -M  3Z  1/3  ......  1 

6  ~  6o(~lS - }  („2  -1  2/3}  W  ,  0  <  N,  M  <  2  .  I 

■n  in  n 


(All) 


This  is  Eq.  (4)  where  Bq  is  defined  in  Eq.  (6). 

In  the  above  formulas,  the  power  0  <  N  <  1/2  characterizes  the  nonideality 
(N  =  1/2  for  y  <<  1,  N  =  0  for  y  »  1),  whereas  the  power  0  <  M  <  1/2  characteri¬ 
zes  the  degeneracy  (M  =  0  for  n  <<  n  ;  M  =  1/2  for  n  >>  n) .  Although  dimensional 
theory  alone  does  not  provide  expressions  for  M  and  N,  it  is  recognized  that  a(y) 
[Eq.  (3)]  and  8(y>  O  [Eq.  (5)]  are  of  magnitude-of-order  one,  i.e.  6  -  (3Z/4im)^3 
for  nonideal  classical  and  quantum  plasmas  of  intermediate  nonideality,  0.1  <  y  <  !• 
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L- 


20  24  2 5) 

In  the  conductivity  tl  ory  of  metals  ’  ’  ,  the  mean  ion  radius 

1/3 

r  =  (3Z/4im)  is  used  widely  as  shielding  length  of  the  ion  potential, 
based  on  phenomenological  arguments.  The  presented  dimensional  analysis 
provides  the  first  mathematical  justification  not  only  for  nonideal  plasmas 
but  also  for  metals. 
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ELECTRICAL  CONDUCTIVITY  OF  NONIDEAL  QUASI-METALLIC  PLASMAS 
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Abstract 


Electrical  conductivity  formulas  are  derived  from  first  principles  for  fully 

ionized  nonideal  plasmas.  The  theory  is  applicable  to  an  electron-ion  system  with 

a  i)  Maxwell  electron  distribution  with  an  arbitrary  interaction  parameter 
2  1/3 

Y  =  Ze  n  /rT  (ratio  of  the  mean  Coulomb  interaction  and  thermal  energies)  and 

2  1/3  2  -1  2/3 

ii)  Fermi  electron  distribution  with  an  interaction  parameter  F  =  Ze  n  /h  m  n 
(ratio  of  the  Coulomb  interaction  and  Fermi  energies).  The  momentum  relaxation 
time  of  the  electrons  in  the  plasma  is  calculated  based  on  plane  electron  wave 
functions  ’nteracting  with  the  continuum  oscillations  (plasma  waves)  through  a 
shielded  Coulomb  potential  Ug(r)  =  e^^exp (-r/5g) /r ,  which  takes  into  account 
both  electron-ion  interactions  (s=i)  and  electron-electron  interactions  (s=e) . 

It  is  shown  that  the  resulting  conductivity  formulas  are  applicable  to  higher 
densities,  for  which  the  ideal  plasma  conductivity  theory  breaks  down  because 
the  Debye  radius  loses  its  physical  meaning  as  a  shielding  length  and  upper 


impact  parameter.  The  conductivity  obtained  for  classical  plasma  is  of  the  form 
*  3/2  1/2  2 

oc  =  a (KT)  /m  e  and  agrees  with  the  ideal  plasma  conductivity  formula 
with  respect  to  the  temperature  and  density  dependence  for  y/Z  -*■  0,  but  its 


magnitude  is  significantly  reduced  as  y/Z  increases.  For  quantum  plasmas,  the 

*  3  2  2 

conductivity  obtained  is  of  the  form  =  o^fi  n/m  Ze  ,  which  shows  that  the 


degenerate  plasma  behaves  like  a  low  temperature  metal. 
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1.  INTRODUCTION 


The  electrical  conductivity  of  nonideal  plasmas  has  been  subject  to 
many  experimental  and  theoretical  investigations- — The  theory  of  the  electrical 

Q\ 

conductivity  of  fully  ionized  plasmas,  based  on  the  Boltzmann  equation,— 

9) 

the  Fokker-Planck  eouation—  (derived  by  expanding  the  binary  collison  integral 
for  the  small,  successive  velocity  changes  of  coulomb  scattering),  or  on  the 
Lenard-Balescu  equation^— ^ (taking  into  account  the  dielectric  properties  of 
the  medium)  is  in  agreement  with  the  experimental  data  only  for  rarefied  high 
temperature  plasmas,  y/Z<<l.The  interaction  parameter  y/Z  is  defined  as  the 
ratio  of  (average)  Coulomb  interaction  (Ze^n1^3)  and  thermal  (KT)  energies 
(n  is  the  electron  density  and  Z  the  ion  charge  number). 


y  =  Ze2n 1 / 3/rT  =  1.670  x  10_3Zn]/3/T 

(cgs-units  are  used  throughout).  The  conventional  transport  calculations 

8— 10)  2—2 ) 

for  classical  ideal  plasmas—  and  weakly  nonideal  plasmas-g—give  an  electrical 
conductivity  of  the  form  o  -  (kT)  3/2/m1/2e?Z£nA,  where  A  =  [  l+(y^)  2  ) 1  -  dj)/P0 

for  JD>>P0*  The  impact  parameter  for  90°  deflections  (Landau  length)  is  pq= 
Ze2/3KT.  The  condition,  AQ>>1  or  ^HA^IO1  is  satisfied  only  for  not  too  low 
temperatures  T  and  not  too  high  densities  n.  Conductivity  formulas  with  this 
Coulomb  logarithm  break  down  for  large  interactions  parameters  y /z  and 
densities  n,  since  the  Debye  radius 

=  [Z/4ir(l+Z)  ]  3 /2y  -1  /2n-^/2 

becomes  of  the  order  and  smaller  than  the  mean  particle  distance  n"3^3  for  y/Z  ~ 

10  and  y/Z  ">  10  ^  ,  respectively.  Thus  d^  loses  its  physical  meaning 

as  an  electrical  shielding  length  and  maximum  impact  parameter.  For  this 

reason,  a  new  shielding  length  A  is  introduced  through  dimensional  analysis—^-. 

Moderately  nonideal  plasmas  with  Y/Z  1  are  readily  generated  through  shock 

1 2* 13) 

wave  compression  and  exhibit  conductivities  of  the  order  o~ 101 -102mho/cm — - , 

which  are  much  smaller  than  those  which  would  be  obtained  by  applying  the 
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conductivity  formula  for  ideal  plasmas  in  the  nonideal  regime.  Although 

there  are  bulk  measurements  of  the  electrical  conductivity  of  nonideal 

Cesium  and  noble  gases  plasmas  available*—  the  formulas  are  valid  only 

for  small  interaction  parameters — —  y/Z.  The  momentum  and  energy  transport 

in  weakly  nonideal  plasmas  (y/Z  <<  1)  was  treated  by  Wilhelm-^  by  means  of 

an  exponentially  shielded  Coulomb  potentia1 ,  which  permits  to  consider  not 

only  short-range  binary  (r^d^)  but  also  long-range  many-body  (r>d^)  inter- 

21 

actions.  This  interaction  model  was  used  shortly  afterwards  by  Rogov— 
for  the  calculation  of  the  conductivity  of  weakly  nonideal  Argon  and  Xenon 
plasmas  with  Debye  shielding.  Later,  Wilhelm-^--'*  applied  his  theory  to 
nonideal  plasmas  by  deriving  a  shielding  length  and  Coulomb  logarithm  which 
are  valid  for  0  <  y/Z  <  1. 

For  moderately  nonideal  plasmas  (0.1  <  y/Z  <  1)  various  phenomenological 

approaches  have  been  used  to  extend  the  conductivity  formulas  of  ideal 

4) 

plasmas,  e.g. ,  Goldbach  et  al—  multiply  the  Debye  length  d^  with  a  free 

parameter  x(p)  which  is  chosen  to  match  the  experimental  data,  i.e.  to 

compensate  for  the  too  rapid  decrease  of  dQ  with  pressure.  A  kinetic  equation 

19) 

has  been  proposed  for  nonideal  plasmas  by  Klimontovicb-—  ,  which  considers 

spatial  correlations  and  temporal  retardation  in  the  collision  integrals,  but 

does  not  take  into  account  many-particle  collisions.  Ebeling  et  al— ^used 

20)  21) 

recently  kinetic - -  and  correlaton  function —  methods  to  derive  a  resistance 

formula  for  nonideal  plasmas,  which  is  applicable  only  for  y/Z  <<  1. 

28) 

Herein,  we  extend  and  apply  the  Bloch — '  transport  theory  to  nonideal 

23) 

plasmas,  based  on  concepts  similar  to  those  used  for  solids — -  and  liquid 

2  A  23) 

metals - - - .  The  application  of  this  model  to  nonideal  plasmas  is  justified 

since  a  plasma  exhibits  a  quasi-crystalline  structure  for  y/Z  >  10  * ,  becomes 
a  liquid  for  y/Z  1,  and  undergoes  a  diffuse  transition  into  a  solid,  metallic 
state  at  a  critical  value  y  .  The  role  of  the  longitudinal  phonons  in  the 
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theory  of  metals  is  assumed  by  the  quanta  of  the  plasma  oscillations  (plasmons). 

The  theory  to  be  presented  provides  a  momentum  relaxation  time  for 

i)  classical  plasmas  (n<h,  fl  =  2(2mnK,I)3/2/'fi3)  with  an  interaction  parameter  y/Z 

2  1  /  3 

and  ii)  quantum  plasma  (n>n)  with  an  interaction  parameter  r/Z(r=Ze  n  / 
fi2m-1n?/3).  The  results  are  compared  with  previous  theories  and 
experiments . 
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II.  ELRCTRON-PLASMON  INTERACTION 


We  consider  the  electrical  conductivity  of  nonideal  plasmas  due  to  the  many- 

bodv  interactions  of  the  electrons  with  longitudinal  plasma  waves  (similar 

to  the  interaction  of  electrons  with  phonons  in  liquids  or  solids).  The 

plasma  under  consideration  is  a  continuum  of  volume  fi  containing  N  electrons 

and  N/Z  ions,  which  exhibits  3N  (high  freouency  branch)  and  3N/Z  (low 

frequency  branch)  characteristic  frequencies  w  (q)  of  longitudinal  oscillations 

s 

(s=e,i).  The  high  frequency  branch  corresponds  to  electron  plasma  oscillations 
and  the  low  frequency  branch  are  the  ion  sound  waves. 

The  motion  of  electrons  in  a  continuum  is  affected  by  the  continuum 
oscillations  (many-body  interactions).  In  ideal  plasmas,  the  change  in 
motion  is  caused  by  binary  collisions  of  the  electrons  with  the  plasma 
particles.  In  nonideal  plasmas,  however,  the  electrons  interact  with  the 
fluctuating  Coulomb  ield  of  all  charged  particles.  Therefore,  this 
interaction  can  be  treated  as  a  scattering  of  the  electrons  by  the 
random  longitudinal  waves  of  the  plasma  continuum,  which  are  thermally 
excited. 

As  in  the  theory  of  metals2—’-  --^we  are  considering  a  free  electron 
model,  which  is  applicable  to  nonideal  plasmas.  For  dense  plasmas  with 
Z  electrons  per  ion,  the  electron  wave  functions  are  approximated  by  plane 
waves  -  exn(iic*r).  The  electron  energy  E  is  given  in  terms  of  the  wave 
vector  k  by  E=fi2  k2  /  2m ,  so  that  the  Fermi  surface  is  spherical. 

Let  oie(q)  be  the  e  eigenoscillation  with  wave  vector  q  of  an  electron  wave 
(e-plasmon)  and,  ui^(q)  the  i*"^  eigenoscillation  with  wave  vector  q  of  an  ion 


sound  wave  (i-plasmon).  Taking  into  consideration  conservation  of  energy  and 
momentum,  i.e.  -RojOc)  =  E'  -  E  and  -ft q  =  p’  -  p,  where  E',  p'  and  E,  p  are 
respectively  the  energy  and  momentum  of  an  electron  before  and  after  a  collision 
with  a  plasmon  of  energy  -fiw(q)  and  momentum  "hq,  we  see  that  an  electron  or  ion 
interacting  with  the  plasma  as  a  whole  can  emit  and  absorb  plasmons  which  are 

•>  -V 

quasi-particles  with  energy  "fiu^q)  an d  momentum  "fiq. 

The  quasi-particles  or  plasmons  obey  Bose-Einstein  statistics,  and  their 
distribution  function  is 

Nq  =  - -  — r -  .  (1) 

<1  r-fitoCq),  , 

6XP  KT~  '  1 

Let  P(k,  k')  be  the  transition  probability  per  unit  time  that  upon  a  collision 
of  the  electron  with  a  nlasmon,  an  electron  in  a  state  k  moves  to  another  state  k'  which 
is  not  occupied  by  any  other  electron.  If  f(k)  is  the  Jistribution  function  of  the 
electron  occupying  the  state  k  and  f(k*>  the  distribution  function  of  the  electron 
in  the  state  k' ,  the  number  of  electrons  which  move  from  the  state  k  to  the  state  k' 
is  (Pauli  principle) 

P(k,  k')f(k)[l  -  f(k')] 

Since  there  exists  always  an  inverse  transition  to  the  above  foreward  interaction, 
the  total  rate  of  change  in  time  of  f(k)  due  .  electron-wave  interactions  is 
obtained  by  summing  over  all  k'  , 


6f(k)  _ 

6t 


{P(k',k)f(k’)[l  -  f (k) ]  -  P(k,  k')f(k)[l  -  f (k ’ ) ] }  .. 


The  interaction  processes  are  calculated  by  perturbation  theory.  According  to 
Akhiezer-^  or  Schiff— ^  the  probability  of  transition  from  an  initial  state  lc  to 
a  final  state  it'  is 

P(t  {■) -l^iv.l2  «<v  v 


E^,  and  E^  are  the  energies  of  the  electron  in  the  states  k'  and  k, respectively. 

(m^ 1 1  is  the  matrix  element  of  the  transition  k"^'.  For  the  absorption  of  a  plasmon, 
|Mkk.|is  proportional  to  fi  ,  and  for  the  emission  of  a  plasmon  it  is  proportional  to 


(N  +  1).  The  more  remarkable  dependence  of  |m  .  , j  is  on  the  Fourier  transform 

M  K.K. 

of  the  potential  U  (r) ,  by  means  of  which  the  particles  in  the  plasma  interact. 
Similar  to  the  matrix  element  of  a  metal  given  by  Sham  and  Zimanil—/ ,  |Mkk’l  f0r 
electron  (s=e)  and  ion  (s=i)  oscillations  with  different  frequency  is 
l^.l2  =  l«  I2  (<i  ’  eq )  2  |Us(q)|2,  s=e,  i 


where 


TIN3 


2m  n  a)  ((J) 
s  s  s 


tsi  2  "fi(Nq  +  1) 

^  2m  n  u  (?) 

s  s  s 


for  the  absorption  of  a  plasmon 


for  the  emission  of  a  plasmon 


(A) 


(5) 


For  plane  waves  normalized  in  a  unit  volume,  the  Fourier  transform  of  U  (r) 

s 


(e  =  -e,  e.  =  Ze)  is  given  by 
s  ‘  4irn  e  |e  I 

|us(q)|  q 

6s  +q 


k'-k 


si  2 


(6) 


is  the  mean  square  amplitude 


m  is  the  mass  of  the  particle  s  ,  n  is  its  density,  |  a 

S  84 

of  the  mode  of  an  oscillation  of  frequency  w  (q) ,  e  is  the  unit  vector  direction  of 

s  q 

the  propagation  vector  q,  and  |Ug(q)|2  is  the  square  of  the  Fourier  transform  of 
the  potential  (through  which  the  electrons  and  ions  are  interacting  in  the  plasma). 
Instead  of  using  a  phenomenological  pseudo-potential— ^  U  (r)  with  adjustable  parameters 


we  describe  the  nonideal  plasma  with  classical  (y/Z)  and  quantum  mechanical  (r/Z) 

interaction  parameters  by  means  of  a  Yukawa  potential  with  a  shielding  radius^-^ 
-1/3 

6  -n  '  (n=n  =N/f2.  n.=N/Zfi): 
s  s  e  1 

e  e  exp(-r/6  ) 

Vr>  =  -T  S  »  s=e,i  .  (7) 

where 


5^  =  (4Trn/3Z)“1/3{l  +  (4it/3Z)1/3  [4*(1  +  Z)  ]-1/2  (y/Z)“1/2} 
with  =  dD  for  y/Z  «  1  and  6  =  (4wn/3Z)-1/,3  for  y/Z  »  1. 


(8) 
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The  plasma  system  consists  of  two  components,  electrons  and  ions  where 


every  charge  (-e  for  s=e  and  Ze  for  s=l)  is  assumed  to  be  uniformally  spread 

over  a  spherical  cell  Qf  Wigner-Seitz  type  of  radius  ^g(s*e,i),  from  which 

the  electrons  are  scattered  in  accordance  with  the  exponentially  shielded 

potential  being  defined  above,  a  radius  is  to  be  defined  in  order 

to  take  into  account  the  electron-electron  interaction.  6  is  assumed  to 

e 

coincide  with  the  shortest  wavelength  ?  2v/q^  for  n^n  and  T /Z>1,  and  q 
is  defined  by  the  conservation  of  the  total  number  of  degrees  of  freedom  of 
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Kor  a] 1  wave  numbers,  0  <  q  <  qg "  211/6^,  che  electron  oscillations  propagate 
with  frequency  u)  =  u>  (q)  >  u)  in  nonideal  plasmas. 


2.  Ion  Oscillations.  The  low  frequency  branch  of  the  space  charge  waves  is 

due  to  ions  sound  waves,  which  are  coupled  with  the  electrons.  Since  the 

311 

ions  are  nondegenerate,  the  frequency  of  the  ion  oscillations  is  — 


where 


and 


W.  (q)  =  v(q)C  q  ,  C 
*  s  s 


(k . kt\i 

V) 


/2 


v(q)  =  [1  + 


Z(k  /k.) 


e  l 


l+Zy^tqS  )2/(36it) 
e  e 


,1/2 

2^3  J  .  n«n  , 


(18) 


(19) 


v  (q)  ?  1,  n>>n  .  (20) 

v(q)  is  a  correction  factor  of  magnit ude-of-order  1,  which  shows  the  influence 
of  the  electrons  on  the  ion  oscillation  (M=ion  mass,  =  5/3  for  the 

ions)  where  0  <  q  <  q^  ~  2 tt /6 ^  . 
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III.  relaxation  time 


The  distribution  function  f(k)  in  Eq.(2)  is  not  symmetric  with  respect  to  the 

+  23)  +  ± 

origin  in  the  k-space  ,  since  it  is  'polarized"  by  the  electric  field  E.  If  E 


is  in  the  x-direction,  f (k)  has  the  form 


f (E)  =  f0(E)  +  *  , 


f  (1-f  ) 
o  o 

KT 


The  Fermi  distribution  function  describing  the  thermal  equilibrium  of  the 
electrons  is 


fo(E)  _  [  1  +  e 


(E-O/KT, 


(22) 


where  S  is  the  Fermi  energy, and  $  is  proportional  to  a  function  CtE)  of  the  energy 

22) 

E  of  the  electrons — 

4>  =  e  |e|v  C(E)  .  (23) 

x 

If  equation  (2)  is  changed  from  the  discrete  summation  over  to  an  integral 
[where  the  volume  fi  is  set  to  unity  because  the  electron  plane  waves  have  been 
normalized  for  a  unit  volume  in  Eq.(6)]  Eq.  (2)  becomes 

The  full  expressions  for  ?<’£,£')  and  P(lcjlc)  are  given  by  Eqs.  (3)  -  (6).  For 
absorption  of  a  plasmon, 

tiN8  q2  „ 

P(£,k')  =  nV  7~T  |us(q)  |  5(E’-E-ftm)  ,  (25) 

msnsws(q)  °  s 

where 

E'  =  E  +  1uus  ,  lc'  =  k  +  q  .  (26) 

For  emission  of  a  plasmon. 
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1 


P(k|£)  =  ~ ?+1)?a  Jus(q)l2  6(E,-E+tW) 
msns  “sU)  ^ 


where 


E’  =  E  -  riws,  k'  =  k  -  q 


Due  to  the  fact  that  in  equilibrium  —  -  0  for  a  Fermi  distribution  function, 
it  is  easy  to  show  that  (detailed  microscopic  balance  for  direct  and  inverse 
interactions) 

w(kjk)  =  p(k,£')  fo0c') u-f0(£)]  =  p(£;£)f0(£)[i-f0(£')]«  w(£,£')  ,  (29> 

Accordingly,  the  linearized  interaction  integral  of  the  Boltzman  equation  is 
obtained  by  substituting  Eq.(21)  in  Eq .  (24).  under  cr nsideration  of  the 
relation  (29)-  Limiting  ourselves  to  first-order  terms  in  <j>,  Eq.  (24)  becomes 


=  i_ -  f  W(k\k)  [$(£')-  $(£)  ]d 3£' 

(2tt)  KT  J 


According  to  Haugl^a  momentum  relaxation  time  exists  in  a  closed  system  of 
particles  in  presence  of  an  electric  field,  and  the  interaction  integral  can  be 


written  as: 


6f  =  f-fo  _ 

<5t  T  ~  T  ’ 

where  the  relaxation  time  t  is  in  general  a  function  of  energy  E.  Hence 


(E)  =  - 


Eq • ( 32)  indicates  how  the  relaxation  time  t(e)  is  related  to  the  collision  term. 

We  distinguish  the  cases:  i)  A  classical  plasma  of  low  density  n<fl,  at  any 

degree  of  nonideality  y/Z  <  l,  for  which  we  expect  the  thermal  energy  of  the  electron 

to  be  much  greater  than  fiu^ ,  i.e.,  KT.  In  this  case,  only  elastic  scattering 


(T 


1 


of  the  electrons  by  the  plasma  waves  is  considered,  E’  =  £.  ii)  A  nonideal  quantum 
plasma  with  T/Z  <  ]  at  high  densities,  n>ft,  and  for  which  in  general, 
we  can  no  longer  neglect  compared  to  the  electron  energy  E,  i.e.  for  which 
liu-s2  KT.  It  is  recognized  also  by  inspecting  the  relations  of  the 
frequencies  to  the  wave  vectors  [Eqs.  (12)  and  (18)]  that  the  effects  of  the 
electron-electron  interaction  (electron-electron  waves)  and  electron-ion 
interaction  (electron-ion  waves)  on  the  relaxation  time  are  to  be  studied  as 
two  distinct  cases. 
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IV.  CLASSICAL  PLASMAS,  n  <  n 


The  two  frequency  branches  of  the  longitudinal  oscillations  are  given  by 
Eqs.  (12)  and  (18)  for  the  electron-electron  interaction  (s*e)  and  electron-ion 
interaction  (s=i).  The  relaxation  times  Teg(E)  (s=e,i)  are  by  Eqs.  (21),  (30) » 


(32 j,  (25),  and  (27) 

(3f  / 3E) <J> 

Tes^^  (6f /6t) 
s 


(33) 


where 


(«). 


8tt  m  n  KT 
s  s 


/ 


q2lus(q)|2 

ws(q) 


-s 

N 

q 


[f0(k’)(l  -  f Q (k) ) ] 6 (E ' 


E  +  -fiWg) 


+  fQ(k)  f  1  -  f  Q (k ' )  ]  6 (E '  -  E  -  HWs)J  [1  -  —  ]d3k'  .  (34) 

For  the  evaluation  of  Eq.  (34)  we  assume  that  a)  E  (classical  plasma), 

S)  E'  =  E,  |  k'  |  =  | k I  (elastic  scattering)  and  y)  isotropic  scattering  (no 

angular  dependence  before  interaction).  With  these  assumptions,  the  transformations 

in  the  Appendix  A  yield  with  Eqs.  (33)  and  (34), 

8/2Tim1/2E3/2tn  n 
s  s 


T  (E)  = 
es 


qs  q5|Ug(q) I 2dq 


-R 


I 


(35) 


0  oj  (q)(exp( 


Riu  ' 

TCT- 


1) 


Before  evaluating  the  integral  in  Eq.(35)  for  the  two  frequency  branches  under 

2 

consideration, we  first  observe  the  behavior  of  |  U  (q)  |  in  Eq.(6).  In  the 
2 

q-domaine,  |ll^q)|  is  bound  between  the  limits 
2  2  2  2.4 

16tt  neeOg  „  ????/. 

- Vt"  1  lu  (q)l  £  16ir  n  e  eV,  s=  e,i  ,  (36) 

(1  +  4u  r  s  s  5 


1 .  Electron-Electron  Wave  Interaction 

For  the  interaction  of  the  electrons  with  the  high  frequency  plasma  oscillation 
of  frequency  wg(q)  TF.q.  (12)  ],  the  relaxation  time  is  evaluated  by  means  of  Eq.(35) 

-Q 

and  the  plasmon  distribution  function  of  Eq.(l)  as 
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(37) 


/2  m3/2E3/2 


T  (E)  =  -Z - r - 

ee  „  K  4  T 

2tt h  ne  I 


re  3  . 

_  f  _ q  dq _ 

\  .  .-2  2.2. ,  L  2  2.5s.  /&(»-  2  2.5$ 

(5  +q)(l+aq)  (exp[^t(l  +  a  q  ) 


!]  ~  1) 


=  h(qe) 


J  2  2  4  hcjp  2  2  5$  ’  0  <  %  * 

o  (1  +  aqV(exP[^(l  +  a  qVl  ~D 


where 


h(^  ^—^272  "  (5e/lT)  • 


(1  +*%r 


a  =  2ir/5 

e  e 


in  accordance  with  the  mean  value  theorem  for  integrals^— '  since  U  (q)  is  bounded 

s 

[Eq.  36]  in  the  interval  (0,q  )  and  q  is  approximated  by  the  mean  value  of  the  two  limits 

e  e 

of  the  integral.  With  a  proper  change  of  variables  the  integral  in  Eq. (38)  is 
evaluated  (Appendix  B)  as 

4mw 

1  =  %  (e  ,  aq  )  , 

e  daic  e  p  e 
e 


n  »  2v  n  2v 

00  B  r  00  Be 

V  ,  '  v  ,  2  '  .  ,  v  2v  e  1  ,  2  .  ,  v  2v  P 

X  (e,aq)  =  l-  -z-e  +  4)  -r= — r  ~,\n  1 - ,  1  -  -re  +  4  )  -77 — , 

e  p  e  5  e  (2v  +  4)2v!  ^4  5  p  (2v  +  4)2v! 


/  00  B.  e 

+  2  y  - 

,  *  v2  3  1  .  (2v  +  2)2v! 

(aq  )  v=l 


/  e  00  B„  e 

+  — 1  _  _£  +  2  y  2_v.  p 

,  '  .4  3  +  4  (2v  + 

<a%)  v-l 


(2v  +  2)2v! 


(aqe) 


-  1'  ' 


B2V  are  Bernoulli  numbers—^,  and 
A  ^  d 

£P  "(l+Z)dD*aqe  "  2lr,/^  T~  ’  ^e  (mKT)1 


e  (UaV)1/2 
P  e 
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where  d^  is  the  Debye  length.  By  Eqs.  (37)  and  (38),  the  relaxation  time  due 
to  the  electron-electron  wave  interaction  is 


t  (E) 
ee 


/2  mV2*  E^2 

e 


8ffne  R  (e  ,aq  ) 
e  P  e 


(42) 


where  R  (e  ,aq  )  -  1  for  n  <<  n. 
e  p  e 


2.  Electron-Ion  Wave  Interaction 


For  the  interaction  of  the  electrons  with  the  low  frequency  oscillations , of  frequency 
^(q)  =  C^q  [Eq.  (18)].  the  relaxation  time  t  .(E)  is  evaluated  by  means  of  Eq.(35) 


and  Eq.(i8)  with  v(q)~l,  and  |U1(q)|  from  Eq.  (5). 
theorem  for  integrals  as  in  Eq.  (38),  one  obtains 

C  m  ^/2  M  E 


By  applying  the  mean  value 


Tei(E)  = 


(43) 


2TfRzne  I. 


where 


0  4 


/ 

PfiC  ql 

\ 

/exp 

s 

.1 

H 

_ 1 

-1) 

0  -  qi  -  qi 


(44) 


’1  ■  ^  ;  <y»>4 


-RC 


(45) 


The  integral  I  is  of  Debye-type^-^ ,  which  becomes  with  b  =*  — —  and  bq  «  x 

i.  l\i  * 


bq 


h  -  <y4  7  J 


i  4 


x  dx 


„  4 
,4qi 


X  1 
e  -  1 


(6i/7T)  Ri^i> 


(46) 


.  2v 


2bqi  "  B?u(bqi) 

^i(bqi)  1  ~  ~5  +  4  l  (2v  +  2)2v! 

v=l 


bq^  <  2it 


(47) 


32) 


Here  again  B ^  are  Bernoulli  numbers — .  The  magnitudes  of  b  and  bqi  follows  from 

the  relations, 

•fiC 


4  x 


•fi 


b 


ITT 


,  bq.  =  2 n /ITT"  X .  /6 ^ 


(48) 


Combination  of  Eqs.  (42)  and  (443  yields,  for  the  relaxation  time  due  to  electron- 
ion  waves  interactions: 

/2  m^K.E*2 

T  .(E)  =  - 7 - .  (A9) 

61  arnZe  Ri(bqi) 

where  R.  (bq.)-l  for  n<<n. 

1  i 

3.  Electrical  Conductivity. 

The  electrical  conductivity  0,1s  related  to  the  energy  dependent  relaxation  times  by 


ne 

a  =  -  <t>  , 

c  m 


1  V  1  4 

-  =  )  -  ,  s=e,i  , 

<T>  L  <T  > 


u  <r 
s  es 


where 


<Tes>  =  ^es(E)f(E)dE 


on 

•  f 


f (E)dE  =  1  , 


2  /E  — E/KT 

For  classical  plasmas,  f(E)  is  the  Maxwell  distribution,  f (E)  ■ - .  e  '  , 

/if  (K l )  * 

which  gives 

mV2K  K.(kT)^ 

<T>  =  - - - - - -  ,  (52) 

mre  [k.R  (e  .aq  )  +  Z  tc  R  (bq,)] 
i  e  p  ^e  e  i  i 

From  Eqs. (50)  and  (52),  we  obtain  the  conductivity  for  a  nonideal  classical  plasma: 
K  K  (kT)V2 

a  =  - LS. - - -  ,  (53) 

C  Ti(2inn)V2e  L 


where 


L  =  k  R  (e  ,  aq  )  +  Zk  R. (bq  )  , 

l  e  p  e  ell 

The  above  results  will  be  discussed  in  section  [VI]. 
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4.  Thermal  Conductivity  of  Nonideal  Plasmas. 


Charge  and  energy  in  a  Coulomb  system  are  transferred  simultaneously 

during  the  motion  of  the  electrons,  no  matter  whether  this  motion  is  caused 

by  an  electrical  field  or  a  temperature  gradient.  In  each  case  both  an 

electrical  and  thermal  current  appear.  The  same  methods  used  for  the 

evaluation  of  the  electrical  conductivity  may  also  be  applied  to  the  thermal 

2  2 

conductivity.  They  are  related  by  the  Lorentz  number  (it  /3)  (K/e)  through 
the  Wiedemann-Franz  law  which  assumes,  however,  that  the  collision  processes 
are  such  that  a  common  relaxation  time  exists  for  an  electric  and  thermal  field, 
which  is  always  satisfied  at  high  temperature  i.e.  'Rco<<KT,  and  hence 


A 


-  -  —  <V 

ct  T 
c 


and  the  thermal  conductivity  is  simply  by  Eqs.  (53)  and  (55) 


A  = 


Vi(71 


/2^C 


(KT) 


5/2 


3eV'2L 


(55) 


(56) 


where  L  is  given  by  Eq. (54). 
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V.  QUANTUM  PLASMA,  n>n. 


For  plasma  densities  n<n,  it  was  proper  to  ignore  -fiwg  compared  with  the 

energy  E  of  the  electrons  which  is  of  the  order  of  magnitude  of  KT.  At  high 

pressures  where  n>>n,  -Korean  no  longer  be  ignored  compared  to  E. 

Our  approach  to  the  quantum  Boltzman  equation  follows  the  Kohler  variational 
33) 

method —  ,  which  is  frequently  used  in  connection  with  the  resistivity  of 
metals.  Combining  the  collision  term  (30)  with  the  electric  field  1  yields 


(| £)-  - ,£)[*(£')  -  *(£>]d3jt'  =  -  v—  °-eE  , 

'0t/s  (2rr)JKTJ  S  x  aE 

Only  v^  appears  since  E  is  assumed  in  the  x-direction,  E  =  (E,0,0). 

According  to  Eqs.(32)  and  (55),  (21)  and  (23),  the  interaction  term 

(<5f/6t)  and  the  function  d>  are 
s 

...  3f 

(«)■  ‘  \~i eE  • 


3f 

*  =  "  eE_3E  VXC(E)  * 

In  order  to  determine  the  function  C(E),  we  follow  the  idea  of  Bloch 
and  expand  C(E)  in  a  power  series  of  (E  -  £)• 


(57) 


(5  8) 

(5  9) 


C(E)  -  l  Cy(E  -0V  -  Co  +  CX(E  -  O  + .  (60) 

where  this  series  is  treated  as  a  trial  function  in  the  Kohler  variational  method 

to  determine  the  coefficients  of  the  series  C  . 

V 

By  Eq. (C-ll) (Appendix  C) ,  the  are  determined  from  the  system  of  equations 

l  CvDuv  =  Nu  (61) 

v 


with 


D 

uv 


(2tt)  KT 


H 


v  (E  -  W(k’,k)(v’(E'  -  o 

x  X 


-  vx(E  -  0V]  d3kd3k, 


(62) 
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f  3f 

-  J  -g|  (E  -  0MF(E)dE 


- .//  I 


E=const 


fdI7dk| 


where  ds  is  a  surface  element. 

Eq. (14)  indicates  that  for  the  high  frequency  branch  (electron  oscillations)  and 
n  >>  n,  the  frequency  ws(q)  depends  on  v^,  the  Fermi  speed  (high  pressure 
quantum  plasma).  Contrary  to  the  classical  case  where  n«n,  in  the  quantum 

34) 

case  the  thermal  energy  of  the  electrons  is  small  compared  with  the  Fermi  energy —  . 
Accordingly,  the  series  exoansion  (59)  is  approximated  by  the  first  term  .  This  assumption 
is  widely  used  in  the  transport  theory  of  metals  ("Bloch  approximation")  and  gives 
good  results  especially  for  the  electrical  conductivity.  Eqs.  (32)  ,  and  (56)- 
(61;  give  as  relaxation  time  in  the  approximation  C(E)  £  Cq  , 


F(E)dE 


>tt)3KtJJ  X 


ch  m  vx.  ,3:*-  ,3£,  D 
.k  ,k)  [  1 - Jd  k  d  k  c 


where  N  and  D  are  to  be  evaluated.  The  integral  N  in  the  numerator  of 
o  00  o 

Eq.  (62)  is  of  Fermi-type  and  since  the  contribution  to  the  intergral  arises  in 

.  the  vicinity  of  E  =  (,  the  limits  have  been  extended  to  0  and  00  .  Following  the 

22  ) 

approximation  made  by  Haug — ,  we  have 

Nq  =  F(E)  +0((KT)2)  ,  (65 

E^C 

2  2  2  2  2 

where  F(E)  is  given  in  Eq.(61),  with  v^  =  k  fcosa/n)  ,  ds  -  k  sina  da  d<f>  and 
2 

dE/dk  = -fi  k/m  (the  angles  a,  and  <j>  are  defined  in  Appendix  A.).  Accordingly, 


-  < 


2  ,  ,  ,  47r/2m 

cos  a  sina  da  a <p  =  — £ 


a»o  ij>=0 
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i.e. 


_  4TL/2m\3/2  3/2  ,  2  2/3  IT 

N0  "  350)  (V  *  5o  (3t  n)  2m  ' 


(67) 


The  denominator  D  ofEq.(S2)  contains  W(k,k’)  defined  in  Eq.(29),  which 

contains  the  transition  probability  P(k',k).  This  latter  function  has 

different  forms  for  the  absorption  f Eq . <25 ^ ]  and  for  the  emission  [Eq. (27) ]  of 

a  plasmon.  First  consider  the  case  of  the  absorption  of  a  plasmon  with 

E'  =  E  +  Hai.  >and  D  =  D+  . 

s  oo  oo 

Using  Eq. (25)  and  (29)  in  the  denominator  of  Eq. (62)  we  get  after 
integrating  over  E 

+  *  rrr  v2;ua(’>i\s 

. 

s  s 


1  - 


44  £  dE  , 
d5v  hv 


(63) 


where  E'  =  E  + -fiojg,  ds'  and  ds  are  elements  of  the  energy  surface.  With 
a  =  angle  between  k  and  the  x-direction,  y  =  angle  between  k'  and  the 
^-direction,  0  =  angle  between  k  and  k*  ,  and  <p,  ip  =  azimuthal  angles 
around  k  and  S'  respectively,  ds  and  ds'  are  expressed  as 
2  E 

ds  =  k  sina  da  d<p  *  — j  2m  sina  da  dip  , 

-fi 


(6?) 


2  k' 

ds'  *  k'  siny  dy  dip  =  —  q  dq  dip  , 


where 


(70  > 


V  •  til 
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k 
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,  .  -fiat  -  -fi  q 

r  cos  Y  -  1  +  2E 


k  cos  a 


2-2/2m 


,2  ,  k'  J2 
ft  v  v  -h 


(71) 
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After  evaluating  the  angular  integrals  and  expressing  N8 ,  |U(q)|  by  their 

q  s 

respective  expressions  in  Eqs. (1)  and  (6),  we  obtain 


7  2  2  +°°  q  n 

32ir  m  n  e  e  C  fS  q  f  (E)[l  -  f  (E  +  TiwB)] 
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35) 


The  double  integral  of  Eq.  (72)  is  evaluated  using  the  mean  value  theorem — '  as 
*r  ?s  q3f  (E)  (1  -  f  (E  +  *uij 


»«,> 
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J  J  u)  iq 
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(73) 


where 


h(q  ) 


(1  +  Ssq  ) 


2.2.2  (5s/7r) 


(74) 


In  Eq.(70),  w(q) (s=e)  is  the  dispersion  of  the  longitudinal  waves  of  the  degenerate 
s 

(n>>n)  electrons  and  f(E)  is  the  Fermi  distribution.  With  the  change  of  variables. 
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which  gives 
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Eq.(72)  becomes 
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(78) 


The  integral  over  n  in  the  H-  integral  is  easily  carried  out  and  is  equal  to 
e/(l  -  e  e) ,  while  the  integral  over  e  is  developed  in  a  series  of  integrals 
of  Einstein  type-12^  , 
e 
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-H 


v  x  , 
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x  .2 
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(.79)  is  evaluated  by  expanding  the  denominator  in  a  series  of  e 
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Substitution  of  (80)  into  (77)  shows  that 

I  a  J 
n 


(80) 


(81) 


where  the  are  to  be  defined  shortly,  and  the  sum  over  n  has  only  a 

few  terms.  For  the  emission  transition,  It'  -  is  -  q  ,  E'  »E  -  ,  D~  is 

s  oo 

obtained  from  Dqo  by  replacing  q  by  -  q  and  u>sby  -  ojg,that  is  e  by  -  e,  and  by 
adding  a  factor  (-l)n+^  to  the  numerator.  Accordingly 

°:0  -  £<-»""*„  j„  •  <«, 

The  total  denominator  D  -  D"  +  D*  of  Eq.(62)  is  by  Eqs.(81)  and  (82) 
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(83) 


The  terms  with  an  even  power  of  e  in  the  integrand  of  Eq.(77)  vanish  according 

to  (&3).  Only  terms  of  odd  power  remain  whose  coefficients  are  the  same  as 

those  of  D+  apart  from  a  factor  2.  Hence, 
oo 


,4  4*4 

2  ne  o 


°°  3ir2mKKT 


(ae  )6  l 


e  ,2  2,2  e, 

(e  -  ep )  e  e  de 

,  e  , ,  2 
(e  -  1) 


(84) 


In  the  range  of  pressures  and  temperatures  in  which  this  equation  will  be  applied. 


the  variables  of  interest  will  be  ep  given  by  Eq(75)  and  aqg  both  of  which  can 


be  expressed  in  terms  of  nondimensional  parameters,  with  q^  -  2ir/ , 
aqe  =  (9/5)1/2(27r)2/3(n/n  )  1/3(y/z)“1/2, 

e  =  2ir25/6(n/n)1/3(Y/Z)1/2  . 

V 


(85) 

(86) 


The  integral  in  Eq.  (83)  is  a  sum  of  integrals  ,  v  =  5,  3,  1,  defined  in 

2  4 

Eq.  (81 ) ,  with  different  coefficients  an>  namely  a^  =  1,  a^  «•  -  2ep  ,  a^  *  ep 


Hence 


D  *  8[JS  -  2e2  J.+  J  ]  =  &F(e  ,  aq  )  , 
oo  j  p  3  pi  pe 


(87) 


where 


„4  4-4 

2  ne  Se 


3 x2  mRKT(aep)6 


(88) 


The  relaxation  time  for  electron-electron  wave  interaction  in  a  degenerate  quantum 
plasma  (Eq.(62))  is  obtained  from  Eq.(65),  (87)  and  (88)  in  the  form 


3Tr3fiKT(aep)^ 

‘  .Vf(c0.  > 


(89) 


At  high  pressures  for  which  ep  >  1,  we  can  easily  terminate  the  series  of  the 

function  F(e  ,aq  )  , 
p  ^e  ’ 
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(90) 


F<ep>aqe)  =  8e3  exp(-ep)[l  +  f^]  +  0(e  */ep) ,  ep  >  1 . 

,  -  -E 

Letting  (1  +  — )  =  a,and  F(e  ,  aq  )  =  8a  e  e  p,  (91) 

ep  pep 

3ir5-fiKTa6e3 

Tee  =  — 4 — 4“^  ’  %  >  1*  for  n>n  and  r/Z?l,  (92) 

8e  a<5  e  p  v 

e 

where  a,  ep,  6^  are  given  by  Eqs.  (14),  (43)  and  (10)  respectively.  For  the 
quantum  plasma,  the  ions  are  still  classical,  and  a  relaxation  time  resulting  from  the 
electron-ion  sound  wave  interaction  can  be  derived  from  Eq.(43)  with  E  5 
(Fermi  energy)  since  •hw  <<  E  for  all  densities  under  consideration.  Thus 

3ir  K'l i3 

Tei  24mZe4  R^bqp  ‘  ^ 

For  quantum  plasmas  bq^  can  be  quite  large,  so  that  the  series  in  Eq.  (47)  does  not 
converge  for  bqi>2Tr.  Henc«  in  the  present  case  where  b^<2x  (bq^  *  x) , 


,  «  B,  X 

R-t(x)  =  1  ~  T  +  4  (2v+2)2ul 


,  x<2tt  , 


R.(x)  = 

■L  4 


164(51-4?  5VX  +  +  + + -iill  x> 

4  *  |v+  2  +  32  +  4  3  +  5  4  I  » 

x  v=l  I  vxvx  vx  vxj 


where  ^(5)  =  ][  k  is  the  Riemann  Zeta-functiorr^—  . 
k-1 

The  effective  relaxation  time  due  to  both  electron-electron  wave  and  electron-ion 

wave  interactions  in  a  nonideal  quantum  plasma  is 

T  .  T  T  , 

ei  ee  ei  .... 

Teff  '  t  +  T  ”  1+G(T)  ’  (96) 

ee  ei 

where  the  relaxat  ion  time  ratio  is 

T  . 

G(T)  =  ~  (97) 

ee 


For  -*■  E  ~  kT  and  'fiu>e<<kT,  Eq.  (92)  and  (93),  averaged  over  a  Maxwell  distribution  give 
the  classical  electrical  conductivity  derived  in  Eq.(53).  On  the  other  hand,  if  we 
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assume  ilto  >kT  and  complete  degeneracy  of  the  electrons,  the  electrical  conductivity  (50) 
of  the  quantum  plasma  with  the  relaxation  time  (96)  is 

°0  =  A  2  2  ~  - -  •  (98) 

‘  2  m  ZeZR,(bq.)(l+G(T)) 

where  for  all  cases  of  interest,  i-e-«  10  cm  snS  10  cm  and  T=10  -10  *X, 

G(T)<<1.  But  for  higher  densities  and  lower  temperatures  G(T)  can  make  a  significant 
contribution.  The  electrical  conductivity  Oq  is  a  linear  function  of  the  density 
n  but  less  sensitive  to  the  temperature  T.  Eq.(98)  agrees  with  the  expression 


for  electrical  conductivity  of  a  low  temperature  metal— Further  discussions 


and  applications  follow  in  section  VI. 
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VI.  APPLICATIONS 


1.  Classical  nonideal  plasmas. 


In  order  to  apply  the  formulas  derived  above,  we  propose  a  study  of  the  electrical 

conductivity  as  a  function  of  the  relevant  parameters  involved.  For  applications 

of  the  theory  to  strongly,  intermediate  and  weakly  nonideal  plasmas,  it  should  be 

noted  that  the  dimensionless  y/2,  £p  =  Hujp/KT,  aq^,  bq^  and  n/n  occuring  in 

Eqs .  (53)  and  (54)  for  the  electrical  conductivity  of  classical  nonideal  plasmas  can  not  be 

varied  independently.  Since  y/z  increase  with  increasing  n  and  decreases  with  increasing  T, 

c  varies  over  a  large  n-T  region  and  hence,  so  does  i  ,  similar  to  y/Z.  Numerically, 

P  ® 

y/z  =  1.67  x  10_3n1/3T_1,  cp  =  A. 328  x  lCfV^T"1 ,  n/n  *  2.07  x  10~16nT*3/2 ,  (99) 


where  f°r  Y/Z  <  ^» 


.,.1/3.  .1/2.  ._t1/2 

aq^  =  (Ait/3)  (tik  (y/Z)  , 

£p  is  given  by  Eq.  (86). 

For  — <1,  aq  >>1  and  hence  from  Eq.(AO)  e  ^  aq  ,  i.e., 
Z  g  e  p  6 


(100) 


=  w11/6  (32k  )1/2  (n/n) 1/3 


(101) 


and  formal,  aqe<l,  Eq.(AO)  reads. 


£  =  2tt(32)1/6(y/Z)1/2[1  +  (Ati/3)2/3ttk  (y/z)_1]1/2(n/h)1/3t  (102> 


In  accordance  with  Eqs.  (102)  and  (103),  it  is  clearly  seen  that  Re(e p,  aqfi)  =  1 
for  all  ^>1,  as  long  as  n<<n. 

So  far  the  parameters  studied  are  related  to  the  electron-electron  interaction. 

For  the  electron-ion  contribution  only  one  characteristic  parameter  bq^  occurs 
in  R^(bq^)  of  Eq.  (A5).  For  this  case  of  relatively  low  density  plasma,  n<n  and  y/ZSl, 
6  given  by  Eq.  (8);  the  expression  of  b^  (3^  =  2tt /6 ^)  as  a  function  of  the  relevant 


dimensionless  parameters  is 


bqA  =  (2ir)  (8ir/3Z)  k^1  2  (m/M)1//2  [1  +  ZQ{y/Z)~1^2  ]'-1(n/n)1^3 

Z  =  (Aw / 3Z)  1/3  [A7t(1+Z)]-1/2 


(103) 


(104) 


and  numerically. 
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(105) 


bq.  =  4.129  x  101 (m/M) 1/2z"1/3  [1+Zo(y/Z)"1/2]_1  (n/n)1/3 
Hence  bq^<<l  for  n<fl  and  R^(b<^)-1  according  to  Eq.  (45).  In  view  of  this  order 

of  magnitude  of  Re(ep,  aqg)  and  R(bq1)  as  n/ft-K)  .the  electrical  conductivity  of  the 
classical  nonideal  plasma(Ea  .  (53)]iS  a  very  weakly  varying  function  of  the  shielding 
parameters  6g(s=e,i)  defined  in  Eqs.  (8),  (10),  and  (11).  The  weak  dependence  is 
attributed  to  the  many-body  interaction  character  of  the  present  theory,  since  the 
electron  is  interacting  with  the  plasma  as  a  whole  rather  than  with  individual 
particles. 

The  electrical,  conductivity  formula  presently  derived  for  classical  nonideal 

plasmas  is  only  weakly  dependent  on  the  density  of  the  electrons  n,  and  goes  as 
3/2 

-  T  .  Hence,  in  this  respect  it  does  not  only  agree  with  the  usual  kinetic 

1-11).. 


theory  results 


(ideal  plasmas)  as  n  gets  very  small  compared  to  n  or  y/Z<<l, 


but  also  with  the  recent  theories  for  nonideal  plasmas  as  well— — 7^. 

The  nonideal  effects  of  the  plasma  on  the  electrical  conductivity  are  then 
governed  by  the  Coulomb  interaction  alone  through  the  nondimens ional  Coulomb 
conductivity  o* ( y/Z)  which  we  define  by 

0  =  o  *(WZ)  .0a>3/2 

°c  °c  1/2  2 

m  e 

and  with  Eq.  (53)  it  is  shown  that, 

*  3/2 
°c  (Y/Z)  =  l(y/Z,  n/n) 


(106) 


L(Y/z,n/fi)  is  defined  in  Appendix  D  as 

L( Y/z,n/h)  =  KiRe(Y/z,n/n)+ZiceR1(Y/Z,  n/fl) 


(107) 


(108) 


where  for  n<n 


Re(y/Z,  n/n)  =  1  +  y^n/n)173;!  +  a^Y/Z)'1] 1/2 (Y/Z) 1/2 


+  a^ (y/Z)  +  . . . , 


(109) 


and 


R^.(y/Z,  n/n)  =  1  -  bo  a+Zo(r/Z)"1/2]'1(n/S)1/3+ b  Il+Z  (Y/Z)“1/2]"2(n/n)2/?.. , 

°  (110) 

The  constants  y^,  bQ  etc.,  are  given  in  Appendix  D.  Comparison  with  the  existing 
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results  such  as  ideal  olasma-^— ^o*(  y/Z)  would  be  o  ( y/Z)*=  [  3  (2  /tt)  1  ^2 /A  ]ln  ta  (Y 

s  o v 

1/2 

with  «  =  3/Z(4tt( 1+Z) )  and  is  valid  only  for  y/Z<<l.  Similar  nondimensional 

o 

Conlomb  conductivities  can  be  expressed  from  the  formulas  which  are  of  special 


applicability  to  only  weakly  and  intermediate  nonideal  classical  plasmas,  i.e.. 


Y  1) 

Y<1.  From  the  first  case  (Ebeling  et  al—  ),  o*(y/Z)  can  be  shown  to  be: 

0* (y/Z)  =  [l-(l/2Za  )(y/Z)3/2]/f£n[a  (y/Z)"3^], 
b  o  o 

17)  * 

where  f  -  1.73.  And  from  the  second  case  (Wilhelm — ),  o  (y/Z)  is: 


(111) 


*,  _  3/ (8tt)  1/2z 

°C(^Z)  M  J2L]  ■ 


r  ,47m  s~l/3 


[41r(l+Z)]‘1/2(y/Z)“1/2 


—  1/3 

where  6^  =d^  for  y/Z<<land  6^=(4nn/3Z)  y/Z>>1-  Table  I  compares  these  different 

nondimensional  Coulomb  conductivities  over  a  wide  range  of  densities  expressed  as  a  ratio  n/fi  at 

4  *  * 

typical  plasma  temperature  T  =  10  K  with  Z=l.  At  high  densities  both  o  and  o  show 

s  E 

19  -3  20  -3 

a  sudden  jump  at  n>10  cm  ,  and  would  give  a  negative  value  at  n>10  cm  and  thereby 

* 

their  applicability  comes  to  an  end  at  these  and  higher  densities,  o  (y/Z)  on  the 

w 

other  hand,  shows  significant  non  ideal  effects  of  the  plasma,  which  can  be  attributed 

”2  /3 

to  the  argument  of  the  logarithm  which  is  an  '  for  y/Z>>l  and  behaves  like  the  ideal 

-1/2 

Coulomb  logarithm  argument  for  y/Z<<l  which  is  "n  .  This  behavior  is  observed  in 

* 

Fig.  1  where  we  draw  a  (y/Z)  as  a  function  of  y/Z  along  with  the  other  Coulomb 

w 

k 

conductivities,  o^y/Z)  shows  evidence  of  its  limitation  to  only  weakly  nonideal 

plasmas  as  it  behaves  (besides  a  very  small  difference  in  magnitude)  identically  with 
*  * 

a  (y/Z)  of  the  ideal  plasma,  a  (y/Z)  of  the  present  theory,  on  the  other  hand, 
s  c 

shows  much  more  evidence  of  the  effects  of  the  nonideality  of  the  electrical 
conductivity.  It  shows  an  important  difference  in  magnitude  for  ^>1  and  yet  converges 

k  k 

to  o  (y/Z)  of  the  ideal  plasma  as  y/Z-K).  a  (y/Z)  does  not  show  this  later  behavior  as 


y/Z^-0  which  is  due  to  the  quantum  mechanical  scattering  involved.  Curves  similar 

29  30) 

to  3  and  4  of  Fig.  1  have  been  reported  in  experiments—  -  -  . 
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Fig.  1.  Dimensionless  Coulomb  Conductivity  Versus  Interaction 
Parameter  y/Z  of  Different  Theories  at  T=10*  °K. 

1:  ideal  plasma,  2:[1]  3:  present  theory,  4:[17], 
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TABLE  I :  The  Nondimensional  Coulomb  Conductivity  o*( y/Z) 

at  T-104  eK,  Z-l. 


n/n 

* 

a 

c 

present  theory 

* 

0 

s 

ideal 

* 

°w 

[17] 

* 

°E 

[1] 

8. 40x10" 7 

0.109 

0.130 

0.045 

0.129 

-6 

8.40x10 

0.102 

0.175 

0.053 

0.173 

8. 40x10" 5 

0.120 

0.268 

0.064 

Q.  261 

-4 

1 

l 

8.40x10 

0.139 

0.571 

0.080 

0.530 

-3 

8.40x10 

0.191 

0.105 

— 

Fig.  2  shows  the  electrical  conductivity  of  classical  nonideal  plasma  presently 
derived  and  given  by  Eq.(53),  versus  the  variable  interaction  parameter  y/Z  and  constant 
density  n.  It  is  seen  that  the  value  of  a  is  slashed  by 

_3 

several  order  of  magnitude  as  y/Z  goes  from  ~10  to  y/Z^10,  through  variation 

of  the  temperature  T  at  a  fixed  density  n.  It  should  be  noted  however  that  a  slower 

increase  is  expected  in  the  electrical  conductivity  with  increasing  y/Z  by 

varying  the  density  n  at 'fixed  T.  This  later  behavior  is  due  to  the  dependence  of  the 

Coulomb  conductivity  Oc  on  y/Z  and  n/ft  through  L(y/Z,n/n).  Such  an  increase  would  not 

—1  1/3  3/2 

have  existed  if  L(y/Z,n/ft)  were  a  constant,  since  y/Z~T  n  and  o  ~a*/( y/Z)  .  A 

c  c 

numerical  illustration  of  this  point  and  of  the  order  of  magnitude  of  is  shown  in 

4 

Table  II  for  a  typical  Temperature  T  =  10  °K  at  different  densities  for  a  H-plasma.  The  relation 
between  the  cg«  units  and  the  practical  units  is  9  x  10^  mho-cm  1  sec 


4 

TABLE  II:  The  Electrical  Conductivity  a at  T  =10  °K  and  Z=l. 


-3 

n  cm 

y/z 

o  mho-cm  ^ 
c 

io17 

0.077 

2.930X101 

1019 

0.36 

3. 6 30x10 1 

1021 

1.67 

8. 590x10 1 

Fig.  2.  Electrical  Conductivity  [mho-cm'^]of  Classical  Nonideal 
Plasma  Versus  the  Interaction  Parameter  y/Z  for  Constant 
Density  n. 


2.  Comparison  with  Experiments 

The  alkali  plasmas  are  being  the  more  commonly  used  in  the  measurements, 
due  to  their  lower  ionization  potential  and  the  material  characteristics  they 

offer  as  wires,  powders  and  vapors.  Although,  a  bulk  of  experimental  results 

.  5,12,15,29) 

has  been  accumulated  over  the  years  on  the  electrical  conductivity  — , 

the  measurements  reported  are  mostly  those  of  weakly  nonideal  plasmas,  i.e. 

0  <  y/Z  <  1.  Here  we  present  in  Fig.  3  the  Coulomb  conductivity  given  by  Eq.(107) 

corresponding  to  a  classical,  weakly  nonideal  Lithium  plasma.  For  Lithium  (first 

ionization  potential  1^  =  5.39  ev)  conductivity  data  are  reported  in  Refs.  36,37. 

By  letting  KT[1+(1/Z)]  =  p/n,  where  p  is  the  plasma  pressure  we  have 


and 


n  =  [  (  z+ifc- 


3/4 

] 


..  3/8,  15/8  ,3/ 2  15/4  n  /7,  -.3/8 

n/n  =  a^p  (y/Z)  ,  /2)(2irm)  e  [l+U/Z)  J 


(113) 

(114) 


Thus  Eq.(107)  becomes 

a*  (Y/Z,p) 


k  tc./(2ir)1/2 

e  x _ 

ttL(y/Z,  p) 


(115) 


where 

L(y/Z,P)  «  (5/3)  {2  +Vf11/3  [l+a0(Y/Z)"1]1/2(Y/Z)19/8p1/8+a1(y/Z)1 

+  U2aL2/3  [l+ao(Y/Z)"1]1(Y/Z)13/8p1/4  -b^173  [1+Zo(y/Z)"1/2  ]_1 
•  (y/Z)5/8  p1/8  +  bla12/3  [l+Zo(Y/Z)'1/2]“2(Y/Z)5/4p1/4  + 

(116) 

with  defined  in  Eq  (114),  aQ,  a^  etc...  are  given  in  Appendix  D,  and 

* 

-  k.  =  5/3.  In  Fig.  3  is  shown  for  two  different  pressures  500  and  125 

atmospheres  at  Z  =  1.  The  theoretical  curves  are  isobars  in  the  low  interaction 

parameter  range  which  agree  with  the  experimental  results  for  a  Lithium  plasma. 

*  * 

Along  with  a  ,  given  by  Eq.  (112)  is  depicted  and  its  agreement  with  the 
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experimental  results  is  also  good.  One  can  observe,  however,  that  shows  a 

slight  increase  of  the  Coulomb  conductivity  with  increasing  the  pressure  p  over 

* 

the  entire  range  0$y/Z<l.  The  pressure  dependence  of  0y  is  given  by: 


A  =(8m/fi2)(W3Z)~2/3  (Ze2)5/4  ~^r  [l+cY_1/2]2 

w  zl/4  pi/ 4 


with  c  =  (4ir/3Z)  [Z/4ir(l+Z)  ]^2  ,  and  hence 


(117) 


o*  =  (3/4z)C2/tt)1/2/^i  .  (118) 

w 

k  k 

It  should  be  noted  that  while  o^,  a w  and  the  experimental  results  of  Fig.  3 

are  in  a  close  agreement  in  the  range  0.3  $  y/Z  $  0.5,  the  Coulomb  conductivity 
* 

og  of  the  ideal  plasma  is  clearly  inadequate  over  the  entire  range  of  the  non¬ 
ideal  regime. 
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3 •  Nonideal  Quantum _ Plasmas 


As  it  was  shown  in  section  V,  in  connection  with  the  relaxation  time  ratio  G(T) 

given  by  Eq.  (97)  ,  the  electron-electron  wave  interaction  contributes  little  to  alter 

20 

the  relaxation  time  of  the  electron  in  high  density  plasmas  in  the  range  n«10 
24  -3 

-10  cm  ,i.e.,  C(T)<<1.  The  electrical  conductivity  is  mainly  due  to  electron-ion 

waves  interaction.  In  order  to  apoly  the  conductivity  formula  (98)  to  nonideal  quantum 

plasmas,  we  express  the  dimensionless  formula  o_/ui  as  a  function  of  the  relevant  para- 

Q  P 

meters,  the  quantum  interaction  parameter  T/Z,  n/n  etc...  From  Eq. (98)  we 
observe  that  the  traditional  logarithmic  term  associated  with  the  ideal  and  weakly 
nonideal  classical  plasmas  is  represented  in  our  formula  by  Q  -  R^(b^i) (1+G(T)) . 

Since  G(T)<<1  in  the  range  of  densities  of  interest,  Eqs.  (94)  and  (95)  permit 
to  express  Q  in  the  form 

Q(n/n)  -  l-CQ(n/n)1/3  +  (^(n/n)273  +  ...,  bqi<2ir,  (119) 

and 

Q(n/n)-  B  (n/n)  ^73~4exp[-A  (n/ft)^"73]  [A..+A(n/A)373+A-(n/A)  73+ 
o  o  j.  2  J 

+A^(n/fl)+Aj.  (n/A)^73]  +  ...  ,  bqi>2u,  (120) 

where  the  constants  A  ,  A1 ,  B  etc  .  .  .  are  defined  in  Appendix  D. 

o  1  o 

In  dimensionless  form.Eq.  (98)  becomes, 

o  /to  =  a*(n/n)  (d„/Xe)  (n/n)l73(T/Z)  1  ,  (121) 

Q  P  Q  * 

where 

o*(n/n)  =  3l73(K./Z)(3/2Tr)576/Q(n/n)  .  (122) 

Xg  is  given  by  Eq.  (41)  d^,  is  the  Thomds-Fermi  screening  length, 

d  =[  (fi2/4me2)  (ir/3  )^73]l/2  given  numerically  with  n/fl  in  Eq.  (99),  and 
F 

d  /X  =  12.30  n~176  T172,  r/Z  =  1. 884x10®  n"173.  (123) 

F  6 

By  Eq.  (115),  0^/(0^  decreases  with  increasing  •  interaction  parameter  T/Z. 

Since  T/Z  is  independent  of  the  temperature  T  and  decreases  with  increasing  n  like 
-1/3 

T/Z'n  ,  the  nonideal  quantum  plasma  (n>n)  becomes  more  ideal  with  increasing 
density  n  in  contrast  to  the  classical  nonideal  plasma  (n<n) ,  which  behaves 
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1/3 

more  "nonideal”  as  n  increases  since  y/Z  -  n  .  Consequently,  the  electrical 
conductivity  of  nonideal  plasmas (n>ft)  has  a  minimum  as  a  function  of  the  pressure 
(or  density  n)  at  some  pressure  p  (or  density  n).  Along  these  lines,  similar 
conclusions  have  been  reported  in  experiments  on  the  electrical  conductivity  of 
alkali  plasmas — ! —  .  This  behavior  should  not  be  attributed  only  to  the  exponential 
electron  density  increase  (with  increasing  pressure)  due  to  (nonideal)  ionization 
potential  lowering,  but  also  to  the  quantum  effects  (n>ft)  resulting  from  electron 
shell  overlapping,  electron  tunneling,  and  electron  transport  in  an  ordered 
liquid-like  ionized  medium.  All  these  effects  contribute  to  the  minimum  observed 
in  electrical  conductivity  data.  The  behavior  of  the  electrical  conductivity 
derived  for  nonideal  quantum  plasmas  in  Eq.(98),  is  shown  in  Fig.  4,  which  gives 
Cq  versus  the  interaction  parameter  T/Z.  A  numerical  illustration  of  the  order 
of  magnitude  of  a^[Eq. (98)]  is  shown  in  Table  III  for  a  typical  temperature  T=104°K 
at  several  densities  of  the  degenerate  electrons  of  a  hydrogen  plasma. 

TABLE  III.  Electrical  Conductivity  of  Quantum 
Hydrogen  Plasma  at  T=1G4°K(Z=1) . 


n[cm  3 ] 

22 

10  * 

CO 

CM 

o 

«— 4 

1024 

Oplmho-cm  3] 

7. 240x10 1 

7. 870x10  2‘ 

9. 480x10 3 
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in  or? 


APPENDIX 


A.  TRANSFORMATIONS. 


The  integration  over  k'  in  equation  (32)  is  transformed  to  an  integration 
over  E'  by  the  following  change  of  variables 


2  2 

,,  -R  k1 
*  =  2m 


<'  2dk'  =  ^  ra3/2  /F  dE’ 


The  square  bracket  of  Eq.(32),  defining  4>  ,  transforms  to 


i-!' 


k'  C(E') 
x 

k  C(E) 
x 


where  k'  and  k  are  related  to  k'  and  k.  In  accordance  with  the  scattering  diagram 
x  x 

we  have 


cos  a  =  ——  ,  k1  =  k'  cos  0.  ,  ,  cos  8,  ,  =  cos  a  cos  0  +  sin  a  cos  (<J) ■  -  4>)  , 

rC  X  iC  X  K  X  X 


With  this  change  of  variables,  the  integrals  in  Eq. (32)  are  over  the  energy  E1 ,  the 
angle  0  and  the  azimulthal  angle  <J>,  where  the  latter  integration  over  <j>  cancels  the 
term  containing  cos((^-  <J>) ,  since  for  each  interaction  we  have  a  constant  q* 


Hence , 


rr  3/2 

d3k’  =  2irsine  d6  ”  -  *F  dE'  , 


—*■  T  f 

Considering  the  triangle  (k,  k  ,  q)  in  the  above  diagram  indicates  that 
=  k'2  +  k2  -  2kk'cos6  , 


where 


|k'|  =  |£|  ,  (1  -  cos0)  «  -^2  *  sin0  de 
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B.  INTEGRAL  I 


The  integral  in  Eq. (35)  is  of  the  form 


L-cy.)4]'  — — C - nr - 

O  (1  +  aZqZ)^[Exp[— P-  (1+  aVn  -  U 


2  2  %  ^wp 

Let  e  =  (1  +  a  q  )  and  ep  =  ,  then  I  becomes 


|e  ,  2  2.2  . 

(e  -  eD)  de 


11  (ae 


(  (£  “  eF 

>6  e  . 

)  J  e  -  1 


P'  e 


Eq.[B2]  contains  integrals  of  Debye^—  type  which  are  given  by 

X  n  T  °°  R  X2v  " 

f  t°dt  n  1  x  r  B2v 

)  ee_  x  X  [n  "  2(nhl)  L  (2v  +  n) (2v) !  J  ’ 

o  v=l 

with  n  =  4,  n=2 ,  and  n=0.  For  n=0  Eq.  [B3]  has  a  logarithmic  solution  and 

-  2~2 
for  n  <<  n  we  have  z  «  1,  Z  «  1,  a  q  >1.  Thus,  I  becomes 


1e  - - 2 -  Re(eP’  ^  ’ 

a  ep 


where 


B  g2V 

Re(V  a<U  =  1  “  5  ee  +  4  JL(2v  +  4)2v!  "  (aqgP 

,  G  “  B„  G2v 

— i  _  _*  +  2  y  2?  e. _  + 

(aq  )2  3  v=l  (2V  +2>2v! 

e  — 


2  v  t 

*  B-> 
y  2v  ? 

(;ml  (2v  +  5)2v! 


i _  [i_!£+2  y  Bzv  cp  1  +  4 

q  )4  L  3  '2v+2)(2v)!J+  > 


— T  tn\- - —  +  e  -  e 

e  e  -  1 


It  is  seen  that  Re(ep,  aqg)  -  1  for  ee«  1,  i.e.  for  n«n. 


mm 


C.  Expansion  of  C(E) 


Introducing  an  integral  operator  L  representing  the  integrals  over  k'  , 

Eq.(55)  can  be  formally  written  as 

3f 

L(vxc<E))  -  -  »,-i|  • 

where  the  function  C(E)  is  to  be  expanded  in  a  power  series  of  (E  -  O •  Define 
*  -  VxC(E)  =  vx  l  Cm(E  - 

In  order  to  determine  the  coefficients  of  the  series  C  ,  we  use  the  Kohler 

u 

33  ) 

variational  principle —  ,  which  makes  use  of  the  expression 
(ip  ,  Lip)  =  fff  4)(Lt|i)d3k  . 


,  L*)  -JJftfWd3! 

■  -  -A  ■  •  -Iff 


3f  ~  + 

v  tj/— rjr  d  k 
x  3E 


where  the  integrations  are  carried  over  a  constant  energy  surface  of  the  electron 
before  the  interaction. The  operator  L  is  defined  by  Eq. (55)  with  the  integration 
over  ,  and  ip  is  the  series  [C2].  Furthermore, 

(ik  Lip)  =  y  C  C  D  I 

vr»  r/  L  y  V  UV 
U,V 


(4)  ,  -  V  -rf)  =  l  C  N 
'  x  3E  L  y  y 


v  (E  -  £)W  W(k',k)[v  (E  -C)V  -  v'(E’  -  C)V]*  d3k  d\' 

X  xx 


f  9  8f 

.T  I  2  o 

Ny  '  "  J  Vx  3E 


(E  -Oyd3k  , 


We  seek  the  maximum  of  (\J>  ,  L^)  with  the  supplementary  condition — , 


0J».  M>)  *  (<l>»  -  vx~gf  ) 


Tccd  =Ycn, 

u  y  v  yv  L  y  y 
y,v  y 


[CIO] 
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by  Eqs.[C5]  and  [ C6 ] .  For  this  purpose  we  add  the  constraint  [C9]  to  Eq.[C6] 

multiplied  by  a  Lagrangian  parameter  X  and  obtain  the  maximum  from  the  condition 


0  =  4—  y  CCD  +  X  J  C  N  =2TCD  +  XN 
dCy  y  v  yv  £  y  y  u  v  yv  y 

Multiplication  by  C  and  a  summation  over  y  yields 


X  =  -2.  Hence,  Eq.[Cll]  reduces  to  the  system  of 


by  comparison  with  Eq. 
equations 


[CIO] 


I  CD 
L  v 


uv 


which  determine  C 


N 

M 

In  particular  Cq 


N 

o 

D 

oo 


[Cll] 


[C12] 
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D.  Evaluation  of  L( y/Z,  n/n) 

In  accordance  with  pq . ( 54 )  L(y/Z, n/n)  =  K^e  1  Z<eR  .  where  R  and  R^  are 
given  by  Eqn.  [B5]  and  (47 )  respectively.  By  Eqs.  (100)  -  (103),  R^(y/Z,n/n) 
becomes  for  n<fi: 

R  (y/Z,n/n)=l  +y. (n/n) 1/3 [1+a  (y/Z)"1 ] 1/2 (y/Z) 1/2 
e  i  o 

+  ja1+y2(n/n)2/3[l+ao(y/Z)_1]|  (y/Z)1 
+  u3(n/n)1/3[l+ao(y/Z)"1]1/2(y/Z)3/2 

+  Ja2+P4(n/n)2[l+ao(y/Z)“]l  }(n/h)2/3[l+ao(y/Z)_1] (Y/Z)2+. 

[Dl] 

where 

a  =  (4tt/3)2/3ttk  ,  a  =-(36)1/3/it5/3k  ,  a9=-A(18n)1/3/K  [D2] 

o  el  e  l  e 

pi=  -4tt(32)1/6/5 

V  28/j  ,2  h 

y3=4/2  /(3ir2)1/3K£  ,  tD3] 

y  =  — ir4/90 
4 

In  accordance  with  Eqs.  (47)  and  (104)  R^(y/Z,  n/ft)  is  given  for  n<fi  by: 

oo  . 

R-(y /Z,n/rt)  =  1  -  b  (n/H)1/3  [1+Z  (y/Z)'1/2]_1  +  I  b  (n/h)  /3  • 

°  v— 1 


•  [l+Zo(y/Z)_1/2]'2v 

where 

b  =  (4/5)  (2tt) 3/2  (TT/3Z)1/3<Y2(m/M)1/2 

O  'v 

b  =  4B  (5b  /2)2v/3/(2v+  2) 2v ! 
v  2v  o 


[D4] 

[D5] 

[D6] 


Based  on  the  definition  in  section  VI,  Q(n/n)  =  R^(1+G(T))  -  ^(n/ft)  for 
G (T) <<  1 .  In  accordance  with  Eqs.  (94),  (95)  and  (8)  with  y/Z>>l,  R^(n/n)  is 
given  by: 

R^(n/n)  =  l-Co(n/n)1/3+4  [  Cv(n/fl)2v/3  ,  bq.<2* 

v=l  1 
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where 


B2v  are 


where 


C  =(8/5)11(2™  )1/2(3Z)'1/3(m/M)1/2 
o  i 

[D8] 

VB2v(5Co/2^/3/(2v+2)2v! 

[D9] 

Bernoulli  numbers  and  v  is  an  index  of  summation.  Furthermore, 

R  .(n/ii)  =  B  (n/n)  -  4  j[  exp [-vAv0 (n/ft)  ^3]  /a  +  A  (n/ft)  ^ 

^  0  v=l  V 

+A  _(n/h)_2/3+A  .  (n/ft)_1+A  (n/il)~^/3i  >  bq,>2»,  [DIO] 

\>3  v4  J  i- 


B  =964 (5) / (5C  /2)4 
o  o 


A 

=  5C  /2 

VO 

o 

A  i 
vl 

=  1/v 

A  0 
v2 

=  4(5Co/2)" 

V 

A  n 

v3 

=  12 (5C 

c 

12) 

1 

-2/v3 

Av4 

=  24 (5C 

0 

/2) 

"V 

Av5 

=  24 (5C 

o 

/2) 

'4/v5 

[Dll] 


[D12] 
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V.  ANOMALOUS  DIFFUSION  ACROSS  MAGNETIC  FIELDS  IN  PLASMAS 


By 

H.  E.  Wilhelm 


Abstract 

The  anomalous  diffusion  transverse  to  a  homogeneous  magnetic 
field  resulting  from  the  interaction  of  the  charged  particles 
with  the  electric  microfields  in  plasmas  with  an  approximate  local 
thermal  equilibrium  is  analyzed  by  means  of  statistical  methods 
based  on  the  Langevin  equation.  The  correlation  functions  of  the 
stochastic  velocity  and  electric  microfields  are  calculated  in  clo¬ 
sed  form,  from  which  an  anomalous  transverse  diffusion  coefficient 
=  (kT/m) /O/*^)  | to |  and  momentum  relaxation  time  t  *  (*^F|w|) 
are  derived  for  particles  of  charge  e  >  0,  mass  m,  and  gyration  fre¬ 
quency  m  «  eB0^m  “  thermal  energy).  Comparison  with  the  pheno¬ 
menological  Bohm  diffusion  coefficient  D®  -  kT/l6|e|Bo  indicates  that 
anomalous  diffusion  in  nonturbulent  plasmas  is  considerably  stronger 
than  in  turbulent  plasmas. 
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INTRODUCTION 


According  to  the  kinetic  theory  (binary  interactions)  of  plasmas-  ,  the  charged 

particles  (electrons  or  ions)  should  diffuse  across  an  external  magnetic  field  with 

2  2 

a  diffusion  coefficient  D,  *  (kT/m)T/(l  +  w  t  ) ,  where  t  is  their  momentum  relaxation 

time  and  u>  =  e Bo/m  is  their  gyration  frequency  (kT  -  thermal  energy,  m  »  mass,  e  >  0 

_2 

*  charge  of  particles).  The  dependence  D,  «  B  on  the  induction  has  been  observed 

o 

2  2 

for  a)  t  >>  1,  e.g.,  in  weakly  ionized  plasmas  in  which  the  charged  particles  inter- 

2) 

act  mainly  with  the  neutrals.—  For  plasmas  with  predominant  Coulomb  interactions, 

the  experiments  frequently  indicate  an  anomalous  diffusion  D,  <*  Bq\  e.g.,  for  low 

3)  4) 

pressure  arcs—  ,  hollow  cathode  discharges  of  ion  thrusters—  ,  magnetically  contained 

fusion  plasmas—^ ,  magnetically  insulated  diodes—^  and  magnetically  Insulated  ion 

beams— ^ .  Bohm  was  the  first  to  investigate  anomalous  diffusion  across  magnetic  fields 

in  low  pressure  mercury  arcs  and  derived  from  the  experimental  data  the  transverse 

8 ) 

diffusion  coefficient— 

D?  =*  kT/16|e|BQ 

Following  the  original  suggestion  by  Bohm  that  the  anomalous  diffusion  is  caused 
by  turbulent  particle  transport  across  the  magnetic  field  lines,  some  not  quite  success¬ 
ful  attempts  have  been  made  at  explaining  "Bohm  diffusion"  within  the  frame  of  plasma 
turbulence  theory^-’— ^ .  Except  for  the  dimensionless  phenomenological  coefficient 
C  ■  1/16,  the  Bohm  formula  can  be  readily  deduced  by  means  of  dimensional  analysis. 

More  recent  experiments  indicate  that  1/20  <  C  <  1/2,  depending  on  the  type  of  plasma 
and  the  level  of  turbulence  present— . 

We  consider  herein  nonturbulent  plasmas  in  a  homogeneous  magnetic  field,  which 
are  fully  ionized  or  contain  neutrals  only  in  small  concentrations  so  that  their 
interactions  with  the  charged  particles  are  negligible.  For  the  diffusion  processes 
the  usual  assumption  is  made  that  they  perturb  the  Maxwell  distributions  of  the 
electrons  and  ions  only  slightly.  Based  on  the  Langevin  equations^-’— ' *for  the 
electrons  and  ions  in  a  homogeneous  $  -field,  the  correlation  function  for  their 
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stochastic  velocity  fields  and  for  the  random  electric  microfields  are  calculated. 
The  thermal  equilibrium  fluctuations  of  the  electric  microfield  are  shown  to  pro¬ 
duce  random  particle  drifts  across  which  result  in  an  anomalous  diffusion  coef¬ 
ficient  D,  «  B-^.  The  anomalous  diffusion  coefficient  is  of  the  magnitude  of  the 
'  o 

g 

maximum  diffusion  coefficient  in  a  magnetic  field,  i.e.  D,  >  D,_ . 

Bohm  diffusion  in  turbulent  plasmas  can  be  treated  in  an  analogeous  manner. 

The  evaluation  of  the  correlation  functions  of  the  macroscopic  velocity  and  electric 
field  fluctuations  in  turbulent  plasmas  requires  ,  however,  different  mathanatical 
methods. 
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THEORETICAL  FOUNDATIONS 


The  stochastic  motion  v(t)  of  a  charged  test  particle  (m,e)  under  the  influence 

■f 

of  the  random  fluctuations  E(t)  of  the  electric  microfield  of  a  quasi-homogeneous  plasma 
is  commonly  described  by  the  Fokker-Planck  equation  or  the  corresponding  Langevin  equa¬ 
tion—’—^.  Within  this  theoretical  approach,  the  diffusion  coefficients  perpendicular 
(D^)  and  parallel  (Du)  to  an  external  homogeneous  magnetic  field  Bq  m  (0,0, Bq)  are 
given  in  terms  of  the  mean  square  particle  displacements  at  time  t  by  — ^ 

-  <Ax^(t)>/2t,  <Ax(t)^>  -  <Ay(t)^>;  D()  =  <Az(t)^>/2t.  (1) 

Eq.  (1)  contains  the  assumption  that  the  plasma  is  isotropic  in  all  planes  perpendi¬ 
cular  to  Bq  (D^  »  D^  =  DL  y  D„).  The  particle  displacements  Ar(t)  -  Jfcv(t)dt  are  de¬ 
termined  by  the  Langevin  equation  for  the  magnetoactive  plasma^-’— 

dv,(t)>Ht  -  %l,(t)  +  v,(t)x|  ]  -  -  v(t)  ,  (2) 

m  o  t 

dv„(t) /dt  *  — Ei.(t)  -  -  v„(t)  .  (3) 

m  t  ii 

In  Eqs.  (2)  and  (3)t  ^(t)  is  the  stochastic  longitudinal  (nonrelativist ic)  mlcrofleld 

produced  collectively  by  the  charged  field  particles  (electrons  and  ions)  of  the  thermal 

equilibrium  plasma  at  the  position  of  the  test  particle  (m,e),  whereas  v(t)*B  is  the 

o 

stochastic  Lorentz  field  generated  by  the  random  motion  v(t)  of  the  test  particle  charge 
across  the  field  lines  Bq.  The  averages  of  the  stochastic  fields  for  an  ensemble  of  . 
test  particles  vanish, 

<v(t)>  -  ft,  <ft(t)>  -  ft  .  (4) 

The  ensemble  averages  over  the  stochastic  fields  are  time-independent,  since  they  are 
identical  with  the  time  averages  over  periods  t»  x,  t„,  in  statictical  equilibrium. 
These  averages  are  calculated  by  means  of  the  velocity  distribution  of  the  test  part¬ 
icles,  which  is  a  Gauss  or  Maxwell  distribution  (Markoff  process)  for  times  t  >>  t ,  V 


Accordingly , 


o 

<v“(t)>  “  kT/m,  i  *  x,y, z .  (5) 

The  test  particle  experiences  a  friction  force  -m(vi(t)/f  +  v„(t)/TM)  as  it  moves 
through  the  random  impulses  from  the  fluctuating  microfield  £(t) ,  i.e.,  x  m  x{^t  } 

and  t  =  are  integral  functionals  of  ^(t).  Since  the  transport  mechanisms  in 

the  directions  perpendicular  and  parallel  to  1  involve  and  are  free  from  induction, 

o 

respectively,  t  and  x  are  necessarily  different. 

In  accordance  with  the  theory  of  the  Langevin  equation,  the  velocity  v (t)  of 
the  test  particle  changes  significantly  during  a  relaxation  time  t,  but  may  vary  only 
by  a  fraction  of  the  change  |v(t  +  t)  -  v(t)|  during  a  single  field  impulse  ^(t)-r^-’— ^ 
In  rarefied  plasmas  this  condition  is  satisfied  if  the  plasma  frequency  w ~uj  is 

t  £• 

large  compared  with  the  relaxation  frequency  (e^  -  dielectric  permittivity  of  vacuum, 

14) 

eQ  *  elementary  charge,  m^  »  electron  mass) — 

-1  -1  ,  2 /  xl/2 

>>  i  ,  ,  o)  -  (ne  /e  m  )  .  (6) 

For  the  Brownian  motion  of  a  macroscopic  test  particle  (colloid)  in  a  viscous 
fluid,  the  relaxation  time  t  is  given  by  the  Einstein-Stokes  formula.  For  the  sto¬ 
chastic  motion  of  a  microscopic  particle,  the  relaxation  frequency  x  ^  has  to  be  cal¬ 
culated  as  a  correlation  integral  of  the  microfield.  The  fascinating  idea  of  the 
interrelation  between  relaxation  time  and  microfield  has  been  first  formulated  by 

Einstein  in  his  investigation  on  the  thermal  equilibrium  between  atoms  and  Planck 

,  15) 

radiation-: — 
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ANOMALOUS  DIFFUSION 


We  consider  a  fully  ionized  plasma  of  n(r)  electrons  (e  -  -e  ,  m  ■  m  )  and 

o  o 

n(r)/Z  ions  (e  =  Ze^,  m  =  M)  per  unit  volume  in  a  homogeneous  magnetic  field  * 
(O.O.B^).  A  weak  density  gradient  Vn(r)  is  assumed  to  be  present  so  that  i)  the 
plasma  can  be  considered  to  be  statistically  isotropic  in  planes  perpendicular  to  Is 


and  ii)  the  associated  diffusion  currents  perturb  the  local  thermal  equilibrium  only 

slightly.  The  Langevin  equations  (2)  and  (3)  represent  two  coupled  differential 

equations  for  the  transverse  components  v  (t)  and  v  (t)  and  one  independent  differ- 

x  y 

ential  equation  for  the  axial  component  v^Ct)  of  the  velocity  fluctuation  of  the  test 
particle,  which  have  the  formal  solutions: 


v  (t)  *  (u  coswt  +  v  sinwt)e 
x  o  o 


-t/T 


+  (e/m)  /  e  ^  S^T[E  (s)cosio  (t-s)  +  E  (s)sinoi(t-s)  ]ds, 
/x  y 

0 


v  (t)  -  (v  coswt  -  u  sinwt)e 
y  o  o 


-t/T 


+  (e/m)/  e  S^T[E  (s)cosu>  (t-s)  -  E  (s)  sinu>  (t-s)  ]  ds , 

y  * 


where 


v  (t)  -  w  e't/T"  +  (e/m)/  E  (s)e_(t_s) /t" ds 
z  o  /  z 


«) j  E 
0 


w  =  eB  /m  >  0,  e  S  0 
o 


(7) 


(8) 

(9) 

(10) 


and  0*  ,v  , w  )  *  v(t»0)  is  the  initial  value  of  the  stochastic  field.  Scalar  raulti- 
o  o  o 

plication  of  the  solution  vector  v(t)  with  the  initial  vector  v(t-0)  and  averaging 
yields,  after  generalizing  the  resulting  correlation  functions  for  "initial"  times  t'  , 
in  the  limit  t  >>  r,  th  : 
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(11) 


r 


<v  (t)v  (t')>  =  <v  (t)v  (t')>  =  (kT/m)costo(t 

xx  y  y 


<v  (t)v  (t‘ )>  =  (kT/m)e 
z  z 


lt-t' |/t„ 


(12) 


2  2  2 

since  <u  >  =  <v  >  =  <w  >  =  kT/m.  The  mean  square  displacements  in  the  directions 
o  o  o 

perpendicular  and  parallel  to  l!  are  proportional  to  t. 


<Ax(t)2>  =  <Ay(t)2> 


kT  7  7 

(t)v  (t1  )>dt  dt*  =  2 — Tt/(1  +io  T  ), 
x  m 


0  0 


t  t 


<Az(t)" 


<v  (t)v  (t1  )>dt  dt*  =  2(KT/m)x  t 
z  z  11 


(13) 


(U) 


by  Eqs.  (11)  and  (12).  Accordingly,  the  transverse  and  parallel  diffusion  coefficients 
defined  in  Eq.  (1)  are 


D,  =  (kT/m)t/(l  +  <o2t2),  D„  =  (kT/m)T„ 


(15) 


D„  and  t„  are  known  from  the  diffusion  theory  for  plasmas  without  magnetic 
fields—^  . 

For  the  evaluation  of  the  transverse  relaxation  time  t  in  Dj^  ,  the  solution  (7) 
is  squared  and  averaged. 


<v  (t)2>  =  <(u  coswt  +  v  simot)2>e  2t^T 
x  o  o 

t  t 


t 


+  (e/m)2e  2 t^T//e  ^r+s^  [<e^  (r)E^ (s)>cos<o(t-r)cos<o(t-s) 


+  <E^(r)E^(s)>sinu(t-r)sinw(t-s)]dr  ds, 


(16) 


since  v(o)  and  ^(t)  are  statistically  independent  and  <E^(r)Ej(s)>  «  0  for  i  ^  j. 
Eq . (16)  becomes  after  integration  over  the  coherence  strip  by  means  of  the  trans¬ 
formation  £  =  r  -  s,  n  =  (r  +  s)/2,  in  the  limit  t  >>  t: 
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where 


1  2  f 

kT/m  =  j  (e/m)  t  /  $(£)cosuj£d£ 


<t(C)  =  <I>(r-s)  =  <E  (r)E  (s)>  =  <E  (r)E  (s)> 

xx  y  y 


(17) 

(18) 


in  view  of  the  isotropy  in  planes  1sq.  Eq.  (17)  gives  the  relation  between  x  and 
Ei(t)  through  an  integral  over  the  microfield  correlation  function  $(£),  which  will 
be  shown  to  decrease  exponentially  with  increasing  |^|. 

The  correlation  function  for  E,(t)  can  be  expressed  in  terms  of  the  correlation 
functions  for  v,(t),  since  these  stochastic  vector  functions  are  interrelated  through 
the  Langevin  equation.  By  Eq.  (2), 


E,  =  -v,  x  Bq  +  (m/e)v,/T  +  (m/e)dv| /dt  (19) 

or 

E,  *  -v,  x  I  +  (m/e)v. /t  ,  t  >>  x,  (20) 

-  -  o  - 

as  known  from  the  general  theory  of  the  Langevin  equation— \  Since  the  correct 
<Ax2>  and  D,  are  obtained  in  the  limit  t  >>  t  (dvf  /dt  -*■  3)— \  the  correlation  func¬ 
tion  of  the  microfield  is  directly  calculated  from  Eq.  (20)  as 

<E  (t)E  (t1  )>  =  B2<v  (t)v  (t1  )>  +  (m/ex)2<v  (t)v  (t1  )> 
xx  o  y  y  y  y 

-  (mBQ/ex)  [  <v^  (t)v^  (t1  )>  +  <v^(t)vx(t'  )>].  (21) 

The  method  used  for  the  determination  of  the  xx  and  yy  velocity  correlations  in 
Eqs.  (11)  and  (12)  yields,  by  means  of  the  solutions  (7)  and  (8),  for  the  asymmetric 
velocity  correlations 


<v  ( t ) v  ( t ' ) >  =*  +(kT/m)sinu(t-t,)e'lt_t' ^/t,  (22) 

x  y 

<vy(t)vx(t'  )>  =  -(kT/m)sinu>(t-t'  ^/t.  (23) 

These  weak  correlations  are  caused  by  the  rotation  of  the  charges  in  the  magnetic 
field,  and  vanish  for  w  =  eBQ/m  ->-  0.  Substitution  of  Eqs.  (11)  -  (12)  and  (22)  -  (23) 
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into  Eq.  (21)  gives  the  desired  correlation  function  for  the  microfield  components 


E  (t)  or  E  (t) , 
x  y 

<E  (t) E  (t')>  =  (1  +  a)2T2)(wT)~2(kT/m)B2cosu)(t-t,)e_'t‘t'  ^  /t  .  (24) 

X  x  o 

Combining  of  Eqs.  (24),  (17),  and  (18)  results  in  an  eigen-value  equation  for  the 
relaxation  time, 

+00 

1  =  J  (1  +  w2T2)T_iy e~^/Tcos2oj£d£  .  (25) 

— oo 

Since  the  integral  is  2t  (1  +  2<j)2t2)  /  (1  +  4uj2t2),  Eq.  (25)  has  the  solution 

(A2  =  1/2;  |  (o  |  t  =  l/v'T  .  (26) 

This  remarkable  result  gives  for  the  microfield  driven  diffusion  of  charged  particles 
(e,m)  transverse  to  a  magnetic  field  Bq  [Eq.  (15)]  the  transport  coefficients: 

D,  =  (kT/m)/[(3//2)|w|]  ,  T_1  =  /2  |u|  .  (27) 


For  electrons  (e  =  -eQ,  m  =  nO  and  ions  (e  =  ZeQ,  m  =  M)  at  the  temperature  T,  the 
anomalous  transverse  diffusion  and  relaxation  frequencies  are: 


D?  *  kT/(3//2)e  B  ,  t_1  =  /2  e  B  /m  , 

o  o  e  o  o  o 

Df  =  kT/(3//2)Ze  B  ,  t"1  =  Jl  Ze  B  /M  . 

•  oo  i  oo 


(28) 

(29) 


This  completes  the  theory  of  the  anomalous,  microfield  driven  diffusion  of  charged 
particles  across  a  homogeneous  magnetic  field  in  nonturbulent  plasmas,  which  are  in 
an  approximate  local  thermodynamic  equilibrium. 
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DISCUSSION 

The  anomalous  diffusion  of  charged  particles  across  magnetic  fields  is  due  to 
the  eigenvalue  character  of  ut  [Eq.  (26)],  which  reflects  the  interrelation  of  elec¬ 
tric  and  velocity  field  correlation  [Eq.  (24)].  The  mathematical  maximum  of  the 
transverse  diffusion  coefficient  is  obtained  from  Eq.  (15), 

dD, /dt  =  (1  -  w2t2)(1  +  u2T2)_1kT/m  =  0  ,  (30) 

as 

D,  =  kT/2m|  oj  |  ,  f"1  =  |u|  .  (31) 

The  actual  diffusion  coefficient  in  Eq.  (27)  is  somewhat  smaller,  but  nearly  equal 

to  the  mathematical  maximum  in  Eq.  (31),  D,  s  D,  since  3//F  =  2.121  >  2.  Thus,  in 

approximate  thermal  equilibrium,  transverse  diffusion  is  practically  optimum. 

8)  8 

Comparison  of  Eq.  (27)  with  Bohm  diffusion—  indicates  that  D,  £  8D7.  Although 
the  phenomenological  coefficient  "8"  is  subject  to  experimental  errors,  it  appears 
safe  to  conclude  that  transverse  diffusion  in  plasmas  is  considerably  weaker  in 
presence  of  turbulence  than  in  approximate  thermal  equilibrium. 

The  condition  (6)  for  the  applicability  of  the  Langevin  equation  to  electrons 
(e)  and  ions  (i)  becomes 


(e)  up  >>  /2|uJ  ,  (i)  up  >>  ,  (32) 

by  Eq  .  (26).  These  inequalities  provide  upper  limits  for  the  magnetic  field  intensity, 

(e)  B  «  (mn/2e  )1/2,  (i)  B  «  (M/Z  )(mn/2e  )1/2.  (33) 

o  o  o  in  o 


(32)  or  (33)  corresponds  for  Z~1  to  the  (illustrative) 

induced  field  <(v,  x  1  )  > ^2  of  the  random  electron  (s  = 
-s  o 

is  small  compared  with  the  r.m.s.  electric  microfield  >f‘- 


cond i t i on 

e^  and  Im 

1 


t  hat 


t  he  r . a  .  s  . 


i* 


1  •  -V. 


plasmas  without,  magnetic  fields 


These  inequalities,  which  can  probably  be  relaxed  from  small  (<<)  to  smaller  (<)  in 
applications,  are  in  general  realized  in  plasmas  in  which  anomalous  diffusion  is 
observed. 

In  experiments,  the  plasma  is  not  always  in  an  approximate  local  thermal  equi¬ 
librium.  Since  the  characteristic  times  for  thermal  relaxation  T  -  T.  ■+•  0  between 

e  i 

electrons  and  ions  and  thermal  anisotropy  relaxation  T±  -  T«  0  in  the  magnetic 
field  are  large  compared  with  the  momentum  relaxation  time  (t),  transient  plasmas 

g  i 

may  be  encountered  with  4  Ti#  In  this  nonequilibrium  situation  the  diffusion 

0  1 

formulas  are  still  applicable  if  one  sets  T  -  TA  in  Eq.  (28)  and  T  *  Tji  in  Eq.  (29). 

Plasmas  with  diffusion  in  weak  density  gradients  are  stable  and  remain  so  even 

for  larger  density  gradients  due  to  the  stabilizing  effects  of  the  homogeneous  magn- 
19) 

etic  fielcF—  .  In  more  complicated  plasma  systems,  e.g.,  low  pressure  arcs  with 

19) 

current  flow  due  to  external  electric  fields,  various  convective  instabilities — ■ 

19) 

and  electron-ion  streaming  instabilities — ■  may  arise.  Anomalous  diffusion  in  un¬ 
stable  and  turbulent  plasmas  will  be  treated  in  a  seperate  investigation. 
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VI.  COLLECTIVE  MICROFIELD  DISTRIBUTION  IN  THERMAL  PLASMAS 


By 

H.  E.  Wilhelm 


ABSTRACT 


The  temperature  and  density  dependent  probability 

-f  ■+ 

distribution  W  *  W(E;T,n)  of  the  collective  microfield  E 

in  fully  ionized,  ideal  plasmas  is  calculated  from  first 

principles  of  statistical  mechanics.  For  typical  ideal 

-1  1/2 

plasmas,  the  average  microfield  ■  [12u(l  +  Z  )nKT] 
is  by  one  to  two  orders  of  magnitude  larger  than  the  char¬ 
acteristic  field  (nearest  neighbor  approximation)  EH 
of  the  Holtsmark  microfield  distribution  P  *  P(E;n).  The 
Holtsmark  theory  and  its  later  extensions  are  shown  to  be 
approximately  valid  for  strongly  nonideal  plasmas  only. 

The  interrelations  between  (average)  kinetic,  interaction, 

collective  microfield,  and  electric  self  energies  is  dis- 

->2  -1 

cussed.  In  particular,  an  equipartition  <E  /8w>  «  3(1+Z  )  x 

nKT/2  among  (average) microfield  and  kinetic  particle  energies 
in  statistical  equilibrium  is  derived  by  means  of  a  thermo¬ 
dynamic  model  of  plasma  formation. 
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INTRODUCTION 


The  probability  distribution  of  the  stochastic  electric  fields  E  produced 
by  the  electrons  and  ions  in  random  thermal  motion,  is  of  basic  interest  for  plasma 
physics  and  for  applications  such  as  the  evaluation  of  the  transport  properties,  the 
Stark  broadening  of  spectral  lines,  and  the  preionization  of  atoms.  Holtsmark  cal“ 
culated  the  probability  ?(E)dE  for  the  dimensionless  electric  microfield  to  be  found 
with  a  magnitude  between  E  and  E  +  dE  for  a  system  of  n/Z  point  charges  Ze  per  unit 
volume  with  the  result  — ^ 

-lr®  -(x/£)3^2 

P(S)  =  (2/ir )E  ]  x  sinx  e  v  '  ’  dx 

0 

where 

E  *  1 E |  / 1 Eh  ,  EH  =  2ir(4/15)2/3Z1/3en2/3 

-5 /2 

This  integral  functional  of  E  shows  the  asymptotic  behavior  P( E)  =»  (3/2 )E  for 

E  -*■  °>  so  that  already  the  second  order  moment  (average  field  energy  density)  does 

not  exist,  <E2>  =■  /  aV{E)EL6E  =  «!  Consecutively,  Gans  — ^  believed  to  have  derived 
0 

the  correct  distribution  in  form  of  an  integral  functional  P (E)  with  converging 

-1/3 

moments  by  considering  non  zero  raddii  rQ  <  (4irn/3Z)  of  the  charged  particles, 
which  make  zero  distance  approaches  and  infinite  fields  impossible.  Both  the  Holts¬ 
mark  and  Gans  distributions  are  independent  of  the  temperature  of  the  field  particles, 
and  describe  essentially  the  fields  of  the  nearest  neighbor  particles  at  an  average 
distance  r  *  (4un/3Z)  since  =  e/r  .  These  theories  are,  therefore,  approxim¬ 
ately  valid  for  strongly  nonideal  plasmas,  in  which  the  thermal  energy  KT  is  negli- 
gible  compared  with  the  average  Coulomb  interaction  energy  Ze  /r  (electron-ion  inter¬ 
actions).  This  conclusion  has  beer,  confirmed  experimentally  by  Vidal— \  who  shoved 
that  observed  Stark  broadening  by  "cold"  microwave  discharge  plasmas  (Ze2n^^  >  KT) 
is  in  good  agreement  with  the  Holtsmark  or  Gans  distributions. 
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The  later  developments  of  microfield  theory  have  been  concerned  with  correcting 
and  extending  the  Holtsmark  theory  under  consideration  of  thermal  effects  and  particle 
correlations.  The  formulation  of  a  rigorous  theory  of  the  microfield  in  thermal 
plasmas  of  arbitrary  nonideality  is  aquivalent  to  the  solution  of  the  general  many- 
particle  Coulomb  interaction  problem,  and  has  not  been  achieved  yet.  However,  sig¬ 
nificant  contributions  to  this  problem  were  made  with  the  help  of  collective  coordi- 

4  5) 

nate  and  discrete  particle  methods  by  distinguished  researchers,  e.g.,  Broyles-’—1 '  , 

Baranger  and  Mozer^’— ^  ,  Hettner  and  Wagner^-’— ^ ,  and  Hunger  and  Larenz^’-^  . 

Most  laboratory  plasmas,  e.g.,  glow  discharges,  arc  discharges  at  not  more 

than  atmospheric  pressure,  and  thermonuclear  fusion  discharges,  are  ideal  systems 

in  which  the  Coulomb  interaction  energy  is  small  compared  to  the  thermal  energy. 

An  interesting  counter  example  for  a  nonideal  plasma  is  the  ball  lightning  phenome- 

18  —3 

non,  which  consists  of  a  highly  ionized  air  plasma  (n  'v  10  cm  )  of  low  tempera- 

ture  (T  'v  10  °K) .  In  ball  lightning,  the  plasma  appears  to  be  in  a  highly  viscous, 

quasi-liquid  state  due  to  the  balance  of  thermal  and  (negative)  Coulomb  interaction 
2  1/3 

energies,  Ze  n  'v  KT,  the  spherical  shape  and  long  life-time  (At  'v  1  sec)  being 

12) 

explainable  by  minimum  energy  considerations. —  In  the  following,  we  are  concerned 
only  with  ideal  plasmas,  for  which  we  derive  the  probability  distribution  of  the 
collective  electric  microfields  from  first  principles,  i.e.  the  results  are  limited 
to  interaction  parameters 

y  =  Ze2n1/3/KT  =  1. 670  x  lO^Zn^T'1*^  1. 

By  the  fundamental  axiom  of  statistical  mechanics  of  ideal  systems  in  thermal 

equilibrium,  all  equilibrium  distributions  can  be  derived  without  consideration  of 

13) 

the  interactions  which  bring  about  the  equilibrium. —  By  extending  this  principle 

for  many-particle  systems  with  discrete  energies  to  continuous  media  with  random 
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4  2 

energy  densities  u  =  E(r/)  / 8n  ,  we  derive  the  probability  distribution  of  the 

collective  microfields  E(r,t)  in  ideal  plasmas  without  approximations.  Previously, 

14) 

we  have  generalized  the  methods  of  statistical  mechanics  for  hydrodynamic —  and  pla¬ 
sma^— ^  turbulence  based  on  the  generalized  entropy  principle  for  nonequilibrium  systems. 
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PROBABILITY  DISTRIBUTION 


Subject  of  the  considerations  is  a  homogeneous,  fully  ionized  ideal  plasma  of 

volume  0  containing  ne  =  n  electrons  and  n^  *  n/Z  ions  per  unit  volume.  In  thermal 

equilibrium,  the  kinetic  energy  densities  of  the  electron  and  ion  components  are 

13) 

given  by  (mg  =  particle  mass) — 


U-l 


m  v  > 
s  sy 


2  “s00* 


N  =  n  0, 
s  s 


e,i. 


(1) 


During  the  random  thermal  motions  of  the  charged  particles,  a  continuous  transform¬ 
ation  of  kinetic  particle  energy  into  potential  electric  energy  occurs,  and  vice  versa, 

due  to  the  particle  interactions  through  their  longitudinal  Coulomb  fields  (transverse 

2 

electromagnetic  interactions  are  negligible  for  m^c  <<  KT) .  By  means  of  a  thermo¬ 
dynamic  model  for  the  formation  of  a  fully  ionized  plasma,  we  demonstrate  that  an 
equipartltion  of  average  random  electric  and  kinetic  energies  exists  in  statistical 
equilibrium  [Eq.  (42)], 


<  £2/8it  >  =  |(1  +  Z-1)nKT  .  (2) 

The  electric  field  ^(r,t)  produced  collectively  by  the  electrons  and  ions  at 

any  point  r  e  n  of  the  plasma  and  the  field  energy  density  u  ■  E(r,t)  / 8 tt  fluctuate 

with  time  t  about  the  average  values  <  >  =  ft  and  <  J2/8tt>  ^  0,  respectively.  The 

proposed  problem  is  to  derive  the  probability  W(£)d^1i  for  finding  the  collective  field 

fluctuation  ?  in  the  volume  element  d^E  =  dE  dE  dE  about  the  point  ?  =  (E  .E  .E  ) 

xyz  r  xJ  y  z 

of  the  field  space  subject  to  the  thermal  equilibrium  conditions  (1)  and  (2). 

In  order  to  determine  experimentally  the  collective  microfield  distribution 
“►  ■"►2 

W(E)  =  W(E  /8n)  in  a  homogeneous  and  isotropic  plasma,  one  would  have  to  measure  the 
fluctuating  field  E  or  the  fluctuating  energy  density  f2/8ir  in  the  vicinity  A^r  of  a 
fixed  field  point  reft  at  consecutive  times  tv  *  v6v,  v  -  1,2,3,...^  within  experi¬ 
mental  errors  At  <<  9  ,  where  9  is  a  time  interval  which  is  large  compared 

v  v  v 


with  the  correlation  time  of  the  stochastic  field  so  that  <lz(r,t)  1j(r,t  +  6  )>  “  0 

v 

(within  these  limitations,  the  magnitude  of  0^  may  be  changed  from  one  measurement 

v  to  the  next  v  +1).  In  a  large  number  of  such  measurements,  N  •*■  »  ,  the  energy 

density  £2/8tt  would  be  observed  times . the  energy  density  ^/8ir  would  be 

observed  N  times,  etc.,  where  £2/8ir  means  an  experimental  value  measured  with  an 
u  a 

error  A(?2/8ir).  The  resulting  step-shaped  energy  distribution  N  «  N  (^/8ir)  is 


a  a  a 


represented  by  the  partition 


12/8tt  I2/87t  ^3 /8ir  ...  $*/8ir  ...  ^/8ir 


where 


N.  +  N0  +  N,  +  . . .  N  +  . . 
12  3  o 


N^/Stt  +  N^/Sx  +  N3t2/8Tr  +  ...N^2/8x  +  ...  N^/Sn  -  N< E2/8n>  .  (5) 

N  is  the  total  number  of  measurements  (N  -*■  ®)  and  N  <  f2/8TT  >  Is  the  total  field  en¬ 
ergy  density  measured  in  the  N  independent  observations.  The  entire  energy  density 
2 

N  <  li  /8t r  >  can  be  distributed  in  a  large  number  H  of  ways  over  sets  {N  }„  of  numbers 

a  N 

13) 

N  .  By  elementary  combinatorics, 
a 


n  -  .V!  /  Nx!  N2!  N3! 


N  !  ...  N  ! 

a  M. 


The  energy  distribution  N  (S  /8ir)  observed  in  statistical  equilibrium  is  the  most 

a  a 

probable  one.  Thus,  N  (?2/87r)is  determined  by  the  condition  for  a  maximum  of  I)  the 

a  a 

number  n  of  realizations  or  ii)  the  entropy  s  •v  Inn  ,  subject  to  the  constraints 
(4)  and  (5). 
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Accordingly,  we  determine  the  probability  distribution  NQ(Ea/8ir)  from  the 
maximum  of  the  function  lnll  =  f(Na)  , 


lnH  =  /V(lrW  -  1) 


-  U 


(InN  -  1) 
a  a  ' 


a=l 


M  1  1 

l  N  E^/8ir  =  tfy(l  +  Z  )nKT 
a=l 


»  #  ->■  00  » 


as  constraints.  Eq.  (8)  holds  by  definition  of  N,  whereas  Eq.  (9)  holds  for  a 
large  number  N  of  measurements  and  the  average  energy  density  <lj2/8Tf*  of  Eq.  (2). 
Addition  of  the  constraints  (8)  and  (9)  multiplied  by  the  Lagrangian  multipliers 
-  X  and  -pto  Eq.  (7)  leads  to  the  compact  maximum  conditions  for  lnH, 


3F(N  )/3N  =  0,  32F(N  )/3N2  <0,  a  *  1,2,...M, 


where 


M  MM, 

F(N  )  =  N(lnN  -  1)  -  y  N  (InN  -  1)  -  A  J  N  -  p  I  N  E^/ 8tt 
a  I.  a  a  a  a  a 

a=l  a=l  a-1 


The  solution  of  Eq.  (10)  gives  the  distribution  of  the  "discrete"  energy  densities 

1)  /8ir  in  the  form 
a 


M  -  A  -  + 

N  =  Ae  a  ,  A=e 
a 


Henceforth,  the  subscript  a  is  dropped  since  1)^  can  be  any  point  1  in  the  field 
space.  The  dimensional  constants  A(X)  and  p  are  then  given  by  the  normalization 
conditions  (8)  and  (9) , 


A fe  uE  /8,T4irE2dE  =  N 
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(14) 


as 


A/(E2/8u)e“^  /8n4irE2dE  =  i4(l+z'1)n  KT 


A  =  [8TT2(l+Z_1)nKT]~3/2^,  p  =>  1/(1  +  Z-1)nKT 


(15) 


For  this  normalization,  which  still  contains  the  number  N  of  measurements,  the 
probability  distribution  (12)  for  the  microfield  energy  density  is 


W^(E2/8tt)  *  — - 


N 


_I2/8tt(1  +  Z_1)nKT 


[8n‘(l  +  Z_1)nKT]3^2 


In  theoretical  applications,  one  is  interested  in  the  probability  dP  =  W(^)d3? 

^  3^. 

for  finding  a  microfield  E  in  the  volume  element  d  E  about  the  point  E  of  the  field 
space,  with  the  normalization  fd?  =  1.  The  corresponding  distribution  function  W(?) 
of  the  collective  microfield  1i  is  obtained  by  renormalization  Qf  ■*  1) : 


2U  +  z'baKT ,  -3/2e-?2/8,  (X+Z-^hKT 


W(E)  =  [  8ir  ( 1  +  Z  )nKT] 


(17) 


This  fundamental  distribution  has  the  form  of  a  Gaussian,  i.e.  all  its  moments  ex¬ 
ist,  e.g., 


<E°>  =  ///  Z°U(t)d*Z  =  t 


(18) 


<E^>  =  ///  t\(h d3E  =  5 


(19) 


<E2>  =  ///  E^W (E)d'*E  =  12tt  (l+Z_i)nKT 


+00 


,3+ 


-U 


(20) 


The  most  probable  (E^)  and  the  r.m.s.  (E^)  collective  microfields  are  by  Eqs.  (17) 
and  (20) 
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(21) 


E  =  [8u  (1  +  Z-1)nKT]1/2 
M 


Ey  =  [  12tt  <1  +  Z-1)nKT]1/2 


For  considerations  concerning  the  fluctuation  of  the  collective  microfield 
2(t)  at  a  point  refi  with  time  t,  temporal  averages  can  be  defined  by 


1*1 


t2  ,  Urn  af  f(t)2dt 

t-h»  2t_t 


The  fluctuation  of  ft(t)  is  defined  by  A?(t)  =  ft(t)  -  ?(t)  with  ft(t)  ■  ft.  In  stat¬ 
ionary  equilibrium,  the  time  averages  are  identical  with  the  ensemble  averages. 

By  Eq.  (17),  the  mean  square  (temporal)  fluctuation  of  ft(t)  is 


AE2  =  f2  -  |f|  =  (3  -  — )4ir(l  +  Z-1)nKT 


2  1/2 

TABLE  1  compares  the  r.m.s.  field  Ew  and  the  r.m.s.  fluctuation  (AE  )  of 

2/3 

the  collective  microfield  with  the  nearest  neighbor  Holtsmark  field  E  *  2*(4/15)  x 

H 

1/3  2/3 

eZ  n  in  dependence  of  the  electron  density  n  for  typical  ideal  plasma  coridi- 

4  2  1/2 

tions  (y  «  1,  T  =  10  °K,  Z  *  1) .  It  is  seen  that  E^  and  (AE  )  are  one  to  two 

18  —3 

orders  of  magnitude  larger  than  E  in  the  range  of  ideal  plasma  densities  n  <  10  cm 
For  these  reasons,  the  Holtsmark  field  E^  represents  a  small  contribution  to  the 
microfield  in  ideal  plasmas.  The  result  Ey  >>  E^  is  readily  understandable  since 
for  ideal  plasmas 


E^/E^  =  (tt/3)(4/15)4/3Z2/3(1  +  Z)"1.  Ze~r  <  y  «  1  . 
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The  probability  for  observing  a  collective  microfield  with  intensity 

E  =  1  E|  in  the  range  between  E  and  E  +  dE  is  P(E)dE  -  W(E)4irE2dE,  where  W(l) 

is  given  by  Eq.  (17).  The  maximum  of  the  probability  density  P(E)  is  P(E^)  ■ 
-1  2  -1  1/2 

4e  [ 8tt  (1  +  Z  )nKT]  by  Eq.  (21).  Accordingly,  the  normalized  probability 
density  is  P(E)/P(EM)  =  (e/4)[8^2(l  +  z'1)nKT]‘1exp[-E2/8ir(l  +  Z_1)nKT]  -  1. 
Fig.  1  presents  P(E)/P(EM)  versus  0  -  E  -  108[Vcm_1]  with  nT  -  1012  -  1022 

_3 

[cm  °K]  as  a  parameter.  This  distribution  is  a  displaced  Gaussian  with  a 

maximum  P(E)/P(Ej>J)  =  1  for  E  =  E^,  which  shifts  to  higher  abscissas  *  [8r  * 
-1  1/2 

(1  +  Z  )nKTJ  with  increasing  nT-values  (pressures).  The  scattering  width 
of  the  distributions  AE  *  ^{8ir(l  +  z  ^)nKT]^2  increases  a(nT)^2  with  increas¬ 
ing  nT-values  (note  logarithmic  scale  of  abscissa) .  The  increasing  quantitat¬ 
ive  importance  of  the  collective  microfield  in  ideal  plasmas  with  higher 
pressures  p  =  (1  +  Z  ^)nKT  is  obvious. 
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ENERGY  RELATIONS 


A  fully  Ionized  plasma  consisting  of  electrons  and  Z  times  ionized  ions  In 

thermal  equilibrium  at  a  temperature  T  exhibits  various  macroscopic  energies,  the 

3  -1 

average  kinetic  energy  <K>  *  — (1  +  L  )nKTO,  the  average  electric  field  energy 

+2 

<U>  =  <E  /8u>  ft,  the  average  interaction  energy  «J»  and  selfenergy  <y>  of  the 

electrons  and  ions.  In  order  to  derive  the  interrelation  between  these  energies, 
the  formation  of  the  plasma  by  an  electric  charging  process  is  considered.  For 
this  purpose,  we  assume  that  the  electrons  and  ions  are  initially  dispersed  at 
infinity  where  they  have  only  selfenergies .  The  plasma  is  then  built  up  by  mo¬ 
ving  one  charge  after  the  other  from  infinity  into  the  volume  fi,  which  requires 
work  against  the  resulting  Coulomb  field  of  the  charges  already  present  in  ft. 

The  thermodynamics  of  the  charging  process  is  illustrated  by  i)  a  reversible 
isothermal  and  ii)  an  adiabatic  or  isotropic  model. 

The  electric  charging  work  expanded  in  moving  N  electrons  of  charge  e^  *  -e 
and  N/Z  ions  of  charge  Ze  against  their  collective  Coulomb  fielf  from  infinity 
into  the  (finite)  volume  0  is  (*  designates  exclusion  of  terms  with  y  «=  v) 


N  .  N  .  „  ^ 

*  «  2|-*e  ->  e  i  x 

2  L.  L .  6  '  ry  "  rv 
y~l  v«l 


i  w  *  « 

A  -  7  I  I 


.  .  N/Z  H/Z 

r1  i  j  *  j  * 


=  1  u-l 


f  _  v  2 1  ■+■!  ->i  |  “  X 

(Ze)  r  -  r 

1  V  U 1 


N  N/Z  . 

-  i  i  ze2\r  -tf1 

y=l  v“l  11  v 


(27) 


->e  i  -f£  X 

where  r^ ’  (r^*  )  are  the  position  vectors  of  the  y-th  (v-th)  electron  (e)  and  ion 
(i)  in  the  volume  ft,  respectively.  The  collective  microfield  of  the  N  electrons 
and  N/Z  ions  at  a  field  point  (r,t)  is  the  superposition 


(28) 
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where  E^r.t)  and  fXr.t)  are  the  individial  Coulomb  fields  produced  at  the  field 
p  v 

point  (r, t)  by  the  p-th  electron  and  the  v-th  ion,  respectively.  By  Eq.  (28),  the 
electric  field  energy  of  the  plasma  f!  is 


u  =  yu  ?(?,t)2d3?  =$  +  '!' 


where 


♦  -  h  frz/w* + fe 

p=l  V>=1  ft  p=l  V“1  ft 


,  ,  N/Z.N/Z.  .  .  , 


,  N  N/Z  .  q 

p=l  v=l  ft 


i  N  0  o  .  N/Z  .  „  , 

'-hyiw*+h\<f?& 

p=l  ft  K  v=l 


(29) 


(30) 


(31) 


are  the  (e-e,  i-i,  and  e-i)  interaction  energy  and  the  (e  and  i)  self energy  of 
the  plasma,  respectively.  Comparison  of  Eq.  (27)  with  Eq.  (29)  reveals  the  inter¬ 
relation 


U  -  V  ~  $>  =  A.  (32) 

Thus,  we  see  that  the  field  energy  U  is  the  sum  of  the  interaction  energy  $  and 
the  selfenergy  T  [Eq.  (29)].  The  charging  work  A  leads  to  an  increase  of  the  in¬ 
teraction  part  $  of  the  field  energy  U  [Eq.  (32)].  The  selfenergy  y  of  the  charges 
is  independent  of  the  spatial  locations  of  the  charges,  i.e.  ¥  is  the  same  before 
and  after  the  charging  process. 

Another  independent  energy  relation  is  obtained  by  multiplication  of  the 

2-+s  2 

coupled  Newtonian  equations  for  the  accelerations  d  r  (t)/dt  of  the  p-th  elec- 

li  >v 

tron  (s  =  e)  and  the  v-th  ion  (s  *  i)  by  their  respective  velocities  vS  (t)  » 

U  »v 
->s 

dr^  ^(t)/dt  and  subsequent  summation  over  all  particles  p  and  v.  The  resulting 
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expression  can  be  brought  into  the  fora  d (K  +  4>)/dt  *  0,  which  demonstrates  that 
the  sum  of  kinetic  (K)  and  interaction  ($)  energies  is  an  invariant  H^, 


where 


K  +  $  -  H 

o 


N  .  ,  N/Z 

J  —  m  (v  )  +  y  —  m  (v  )' 

il  2  e  2  i  v' 


(33) 

(34) 


and  is  defined  in  Eq.  (27) .  Eq .  (33)  expresses  the  conservation  of  kinetic  K 
and  interaction  $  energies  in  a  plasma,  in  which  the  electrons  and  ions  interact 
by  longitudinal  Coulomb  fields. 

The  thermodynamic  functions  of  the  plasma  depend  in  general  on  the  volume  R, 

the  number  Ng  of  particles  in  fl,  and  the  particle  averages  of  the  random  kinetic 

1  -»2  -*2 
energies  —  i“svs  and  the  random  field  energy  densities  E  /8ir.  Accordingly,  we 

assume  Uth  =  Uth(T,e,N  )  for  the  thermal  energy  and  S  *  S(T,e,N  )  for  the  entropy, 

s  s 

where 


3KT/2 


1 

<2 


-*■2 

m  v  > 
s  s 


e  =  <E  / 8tt>  . 


(35) 


For  plasma  formation  by  isothermal  reversible  charging,  the  volume  R  is 

embedded  into  a  heat  bath  of  temperature  T.  The  transfer  of  dN&  charges  eg  from 

infinity  into  the  cavity  R  requires  on  the  average  the  charging  work  dA  * 

8 

d  <U  -  T>  =  d  <U>  by  Eq.  (32),  and  their  thermalization  at  a  temperature  T  con- 
s  s 

th  3 

sumes  on  the  average  the  energy  dU  =  —  KTdN  (s  -  e,i).  The  difference  of  these 

s  t  s 

energies,  dQg,  is  supplied  by  the  heat  bath.  Summation  over  "s"  yields,  in  acc¬ 
ordance  with  the  first  law  of  thermodynamics 


dQ  -  duth 


d<U> 


(36) 
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since  no  other  than  electric  charging  work  is  performed  on  the  system  (dft  *  0) . 
The  associated  entropy  dS  »  dQ/T  is  a  complete  differential. 


ds  =  ?  3T  (ljth  ‘  <U>)dT  +  T  a7  (ljth  -  <U>)de  +  7  l  3|  <uth  -  <U>)dNg-  (37) 

sMe,i  s 

Application  of  the  condition  3£3,j,S  =  to  Eq.  (37)  yields  the  partial  differ¬ 

ential  for  constant  N  and  T, 

s 

3Uth/3e  =  3<U>/3 e  .  (38) 

Since  Uth  =  0  for  e  =  0  (no  thermal  energy  in  ft  before  charging) ,  the  integral 
of  Eq.  (38)  is 

Uth  =  fte  .  (39) 

Eq.  (39)  could  have  been  derived  by  other  thermodynamic  plasma  formation 
processes,  e.g.  by  adiabatic  charging  of  the  cavity  ft.  In  this  case  dQ  »  0, 
and  by  Eq.  (36) 


dQ  =  dUth  -  d<U>  *  0  : 


<U> 


,th 


(40) 


Finally,  <U>  can  also  by  determined  as  that  equilibrium  value  which  maximizes  the 
entropy, 

dS  =  T_1[dUth  -  d<U>]  -  0  :  <U>  -  Uth  '  (41) 

Eqs. (39)-(41)  indicate  that  an  equlpartition  between  thermal  energy  and 
average  microfield  energy  exists  in  statistical  equilibrium.  This  fundamental 
result  is  explicitly 


|  (1  +  Z_1)NKT  =  ft<f2/87T> 


ill 


(42) 


The  virial  equation  for  the  fully  ionired  plasma-^^,  the  averages  of  the 
field  energy  equation  (29),  the  energy  conservation  equation  (33),  and  the  ener¬ 
gy  balance  equation  (42)  represent  four  independent  equations  for  the  average 
kinetic  energy  <A>,  the  average  field  energy  <U>,  the  average  interaction  ener- 


and  the 

energy 

invariant  H  : 

o 

<K>  + 

— <  <t>  >  = 
2 

f  pQ, 

(43) 

<U>  - 

<4>>  = 

y  , 

(44) 

<K>  + 

<$>  = 

H 

(45) 

o 

<U>  = 

<K> 

(46) 

The  pressure  p  of  the  plasma  is  assumed  to  be  known  (measurable).  The  self 
energy  ¥  =  <¥>  is  independent  of  the  random  motion  and  spatial  distribution  of 
the  particles  in  the  case  of  ideal  plasmas,  and  can  be  calculated  from  the  charge 
distribution  in  the  electrons  and  ions. 

As  an  illustration,  the  plasma  energies  are  calculated  by  means  of  Eqs. 

(43)  -  (46)  for  the  case  that  p  and  ¥  are  known:  <U>  * <K>  m  pft  +  ¥/3  >  0,  <4» 

=  pft  -  2¥/3<  0,  Hq  =  2pfl  -  ¥/3  >  0.  By  definition,  ideal  plasmas  are  systems 
with  weak  Coulomb  interactions,  which  preserve  their  electrical  neutrality  and 
collective  behavior.  In  applications  of  Eqs.  (43)  -  (46)  to  ideal  plasmas,  it 
should  be  kept  in  mind  that  the  interaction  energy  <$>  is  always  small  but  never¬ 
theless  nonvanishing, 

0  <  Y  ^  |<$>|</i:>  «  1.  (47) 


112 


CONCLUSION 


In  ideal  plasmas,  the  distribution  of  the  collective  microfields  is  strongly 
temperature  and  density  dependent.  For  typical  temperatures  and  densities  of  ideal 
plasmas,  the  r.m.s.  collective  microfield  is  by  orders  of  magnitude  larger  than  the 
characteristic  Holtsmark  field.  The  temperature  independent  Holtsmark  theory  is 
approximately  valid  for  strongly  nonideal  plasmas  only,  for  which  thermal  effects 
are  negligible.  In  statistical  equilibrium,  a  balance  among  (average)  kinetic  par¬ 
ticle  and  collective  microfield  energies  exists,  which  is  independent  of  the  thermo¬ 
dynamic  process  of  plasma  formation. 
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VII.  FREE  ENERGY  OF  NONIDEAL  CLASSICAL  AND  DEGENERATE  PLASMAS 

By 

H.  G.  Wilhelm  and  A.  H.  Khalfaoui 
Abstract 

A  quantum  statistical  theory  of  the  free  energy  of  a  nonideal  electron-ion 
plasma  is  developed  for  arbitrary  interaction  parameters  0  <  y  <  y 

c 

(y  »  Ze2n  ^3/KT  is  the  ratio  of  mean  Coulomb  interaction  and  thermal 
energies) ,  which  takes  into  account  the  energy  eigenvalues  of  (i) 
the  thermal  translational  particle  motions,  (ii)  the  random  collective 
electron  and  ion  motions,  and  (iii)  the  static  Coulomb  interaction  energy 
of  the  electrons  and  ions  in  their  oscillatory  equilibrium  positions. 

From  this  physical  model,  the  interaction  part  of  the  free  energy  is 
derived,  which  consists  of  a  quasi-lattice  energy  depending  on  the 
interaction  parameter  y,  and  the  free  energies  of  the  quanitized  electron 
and  ion  oscillations  (long  range  interactions).  Depending  on  the  degree 
of  ordering,  the  Madelung  "constant"  of  the  plasma  is  ct(y)  “  a  for  y  >>  1, 
a(y)  5  a  for  y  >  1,  and  a(y)  «y  ^  for  y  <<  1,  where  a  -  1  is  a  constant. 

The  free  energy  of  the  high-frequency  plasmons  (electron  oscillations)  is 
shown  to  be  very  small  for  y  >  1,  whereas  the  free  energy  of  the  low- 
frequency  plasmons  (ion  oscillations)  is  shown  to  be  significant  for 
y  >  1,  i.e.  for  proper  nonideal  conditions.  For  weakly  nonideal  plasmas, 
y  <<  1,  both  the  electron  and  ion  oscillations  contribute  to  the  free  energy. 
Thus,  novel  results  are  obtained  not  only  for  proper  nonideal  (y  >  1) 
but  also  for  weakly  nonideal  (y  <<  1)  plasmas.  From  the  general  formula 
for  the  free  interaction  energy  AF  of  the  plasma  for  0  <  y  <  y  simple 
analytical  expressions  are  derived  for  AF  in  the  limiting  cases,  y  >>  1, 

Y  ?  1,  and  y  <<  1.  Applications  to  astrophysical  problems  are  discussed. 
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INTRODUCTION 


In  the  classical  work  of  Debye  and  Hueckel  on  electrolytes,  the  total 

Coulomb  interaction  energy  is  calculated  from  the  continuum  theoretical 

picture  of  every  ion  interacting  with  its  surrounding  space  charge  cloud. 

Using  more  sophisticated  methods,  similar  results  were  obtained  for 

weakly  nonideal  plasmas  (y<<l)  by  Mayer1  (cluster  expansion),  Ichikawa^" 

3  4 

(collective  variable  approach  ),  Vedenov  and  Larkin  (graphical  density 
expansion),  and  Jackson"*  (hydrodynamic  continuum  interaction  model). 

Based  on  different  methods  and  approximations,  investigations  of  moderately 

(y  >  1)  and  strongly  (y>>l)  nonideal  plasmas  were  given  by  Berlin 

6  7  8 

and  Montroll  ,  The imer  and  Gentry  ,  Ecker  and  Kroell  ,  Ebeling,  Hoffman 

9  10 

and  Kelbg  ,  and  Varobev,  Norman  and  Vilinov  ,  respectively. 

In  spite  of  differences  in  the  theoretical  approaches,  the  leading 

terms  of  the  analytical  results  for  proper  nonideal  plasmas  (y>l) 

give  essentially  the  same  formula  for  the  free  plasma  energy,  AF/NKT  -  -ay  + 

2  1/3 

blny  +  c,  due  to  Coulomb  interaction,  where  y  =  Ze  n  /KT  is  the 

ratio  of  electron  -  ion  interaction  energy  and  thermal  energy,  and 

a,  b,  c  are  constants  depending  on  the  respective  approximations 

and  assumptions.  The  thermodynamic  functions  of  strongly  nonideal 

plasmas  (y>>l)were  also  determined  with  the  help  of  Monte  Carlo  and 

computer  methods  by  Brush,  Sahlin,  and  Teller11,  Hansen1^,  Votobev  et  al15 
.  14 

and  Theimer  ,  respectively.  Although  computer  methods  provide  limited 
physical  insight,  they  are  useful  for  checking  the  quantitative  validity 
of  analytical  theories  . 

At  sufficiently  high  electron  densities,  for  which  Yj  1,  classical 
statistical  theories  fail  due  to  thermodynamic  instability15,  which  is 


inhibited  by  quantum  mechanics.  The  classical  plasma  pressure  would  collapse  for 
y  >  1  dua  to  the  negative  electron-ion  interaction  energy,  whereas  in  reality 
the  pressure  remains  positive  in  a  plasma  due  to  the  Fermi  pressure  (exclusion 
principle)  of  the  electrons.  For  these  reasons,  we  present  herein  a  quantum- 
statistical  theory  for  nonideal  plasmas  based  on  concepts  similar  to  those 
used  by  Debye  for  solids^^ .  The  application  of  this  model  to  proper 
nonideal  plasmas  (y  >  1)  is  justified  since  a  plasma  exhibits  a  quasi¬ 
crystalline  structure  for  y  >  0  before  it  undergoes  a  diffuse  transition 

into  a  solid,  metallic  state  at  a  critical  value  y  .  The  roll  of  the 

c 

longitudinal  phonons  of  the  Debye  theory  is  assumed  by  the  quanta  of  the 
plasma  oscillations  (plasmons)  in  the  case  of  the  quasi-crystalline  plasma. 

The  theory  is  also  applicable  to  weakly  nonideal  conditions,  since  the  quasi¬ 
lattice  energy  reduces  for  weak  ordering,  y  <<  1,  to  the  free  interaction 
energy  of  weaklv  nonideal  plasmas. 

The  theory  to  be  presented  takes  into  consideration  (i)  the  energy 
eigenvalues  of  the  random,  collective  electron  and  ion  oscillations  and 
(ii)  the  static  Coulomb  interaction  energy  (quasi-lattice  energy)  of  the 
electrons  and  ions  in  their  oscillatory  equilibrium  positions.  Thus, 
all  significant  long  and  short  range  Coulomb  interactions  are  considered. 

The  results  are  applicable  to  arbitrary  nonideal  plasmas,  0  <  y  <  y  ,  where  y 

c  c 

is  the  critical  ordering  parameter  at  which  a  phase  transition  into  a 
solid  metallic  state  occurs. 
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PHYSICAL  FOUNDATIONS 


Subject  of  the  theoretical  considerations  are  quasi-homogeneous  high- 
pressure  plasmas  consisting  of  electrons  of  charge  -e  and  density  n  *  N/V  and 
ions  of  charge  +Ze  and  density  n/Z  =  N/ZV,  with  typical  densities  in  the  range 

20  _3  24  -3  3  4  o 

10  cm  <  n  <  10  cm  and  temperatures  of  the  order-of-magnitude  T~10-  10  K. 

2  -h 

For  these  conditions,  the  Debye  radius  D  =  [Anne  (1+Z)/KT]  is  D  ■  6.901  x 

[T/n(l+Z)2  c  10  cm, i . e.,  u  is  smaller  than  the  atomic  dimension  and  the  number 

3  —8 

of  particles  in  the  Debye  sphere  would  be  N^=4irnD  / 3  <  1  for  D<10  cm  and 
24  -3 

n<10  cm  .  It  is  seen  that  the  concept  of  Debye  shielding  completely  breaks 

down,  and  statistical  theories  containing  the  Debye  length  as  a  characteristic 

parameter  would  be  physically  meaningless  for  high  density  plasmas. 

The  nonideal  behavior  of  plasmas  is  determined  by  the  interaction  parameter 

2  1/3 

Y,  which  is  the  ratio  of  the  Coulomb  interaction  energy  ~Ze  n  and  thermal 


energy  -KT, 


Y  =  Ze2n1/3/KT=  1.671  *  10_3Zn1/3/T. 


(1) 


It  follows  that  0.5Z  <  y  S  15Z  for  l02<3cm  3  <  n  <  102^cm  3  and  T  -  10^  °K.  For 
Y  =  1,  the  nature  of  the  plasma  changes  from  a  "thermally  expanding"  (y<1)  to 
an  "electrostatically  contracting"  (y>1)  plasma.  For  Y>1>  the  collapse  of  the 
plasma  due  to  Coulomb  attraction  between  electrons  and  ions  is  inhibited  by  the 
Fermi  pressure  of  the  electrons,  i.e.  by  the  quantum  mechanical  exclusion  prin¬ 
ciple.  Thus,  in  the  region  0  <  y  <  Y  the  plasma  undergoes  a  diffuse  transi- 

c 

tion  from  a  nonideal  classical  plasma  (y  <  1)  to  a  quasi_crystalline  plasma 

(1  <  Y  <  Y  )>  with  an  incomplete  ordering  comparable  to  that  of  a  liquid, 
c 


An  understanding  of  strongly  nonideal  plasmas  has  been  attempted  via  the 
model  of  discrete  interacting  particles  in  a  dense  gas .  For  the  above 
reasons,  however,  it  appears  to  be  more  adequate  to  calculate 


the  thermodynamic  functions  of  proper  nonideal  plasmas  from  the  picture  of 
collective  electron  and  ion  oscillations.  In  this  approach,  the  free  inter¬ 
action  energy  is  due  to  the  static  Coulomb  interaction  of  the  electrone  and 
ions  in  their  "equilibrium  positions"  (Madelung  energy)  and  their  oscillation 
energies  about  the  equilibrium  positions  (plasmon  energies). 

Since  the  plasma  volume  V  contains  N  electrons  and  N/Z  ions,  there  exist 
3N'  (high-frequency  branch)  and  3N/Z  (low  frequency  branch)  characteristic  fre¬ 
quencies  of  longitudinal  oscillations.  Each  plasma  oscillator  of  frequency 
can  only  have  the  energy  (n^  +  *2)110^  ,n^=0, 1 , 2  , .  . .  ,  so  that  the  energy  E{i} 
of  a  plasma  state  with  n.  plasmons  of  frequency  u)^  is 

E{i}  =  £n.-ftu).  (2) 

{i}  1  1 

where  { x }  designates  the  entire  set  of  given  eigenf requencies  Ac¬ 

cordingly,  the  partition  function  Q  of  the  longitudinal  plasma  oscillations  is 


Q 


-  Ie-E{i>/KT  =  11  i/(1-e^VKT). 

{ i}  Ui 


(3) 


From  0,  the  thermodynamic  functions  such  as  the  pressure,  internal  energy, 
entropy,  etc.,  are  derived  in  the  usual  way,  e.g.,  the  free  energy  of  the 
plasmons  is 

F  =  -KT  In  Q  =  KT  £  ln(l-e'fiuJi/KT)  .  (4) 

(i} 


In  the  limit  V-*»  ,  the  discrete  eigenfrequencies  uk  are  replaced  by  con¬ 
tinuous  ones,  w=u)(k),  in  accordance  with  the  dispersion  law  for  space  charge 
waves  of  wave  length  A  =  27r/k,  0  <  k  <  k. 

1.  Electron  Oscillations.  The  high-frequency  branch  of  the  space  charge 
waves  is  due  to  longitudinal  electron  oscillations  .  Their  frequency  u  is  for 
classical  (n  <<  n)  and  completely  degenerate  (n  >>  n)  electrons  given  by^-^ 

[1  +(ice/4ff)ZY~1  (k?e)2J  ,  n<<n  ,  (5) 
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(6) 


w2  =  Up  [1  +  2(J“"  ("g")  ^  (Z/y)  (k  F  )  2  ]  ,  n  >>  ii  , 


where 


n  =  2(2TrmKT/h2) 


2  \  % 


u)  =  (4ime2/m) 

P  -v3 


v2 


r  =  n 
e 


are  the  critical  electron  density,  the  plasma  frequency,  and 

the  mean  electron  distance  (k.S  c  /c  of  the  electrons,  and  m  is 

e  p  v 

their  mass).  Since  k  -  2ir/r  (oscillations  with  X  <  r  are 

max  e  e 

physically  inconceivable),  the  electron  oscillations  propagate, 


to  =  <o(k)  >  Up  ,  in  nonideal  plasmas. 


(7) 

(8) 

(9) 


2.  Ion  Oscillations.  The  low-frequency  branch  of  the  space 
charge  waves  is  essentially  due  to  ion  sound  waves.  Since  the  ions 
are  presumed  to  be  nondegenerate,  the  frequency  of  the  ion  oscillations 
is  given  by^^ 

(a)  =  fiMO^KT/M)1^  (10) 

where 

6(k)  =  [l  +  - — - - 1  '2  ,  n  «  n  ,  (11) 

L  1 +  (iCgMiOZy  (krg)2J 


<$(k)  =1  ,  n  >>  n  ,  (12) 

is  a  correction  factor  of  magnitude-of-order  1,  which  shows  the 

influence  of  the  electrons  on  the  ion  oscillations  (M  =  mass,  ^  * 

c  /c  of  the  ions). 

P  v 

In  weakly  nonideal  plasmas,  y  «  1,  the  electron  sound  waves  are 
strongly  damped  for  wave  lengths  X  <  D,  due  to  trapping  of  the  resonance 
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electrons  with  thermal  speeds  comparable  to  the  wave  speed.  For  proper 

nonideal  plasmas,  y  >  1,  the  number  of  particles  in  the  Debye  sphere 

—  8 

4ttD3/3  is  no  longer  large  compared  with  one  and  D  <  10  cm  is  smaller 

than  atomic  size,  so  that  thermal  Landau  damping  is  no  longer  feasable. 

For  this  reason,  electron  oscillations  should  exist  for  wave  length 

X  >  r  if  r  >  D. 
e  e 

The  ions  are  nondegenerate  since  n,  <<  (2-irMKT/h2 )  3/2  for  the  n  -  T 
region  under  consideration.  The  electrons  are  considerable  degenerate 
for  n  >  n  by  Eq.  (7),  i.e.  their  kinetic  energy  is  essentially  given  by 
the  Fermi  energy  E  =  di2  (3ir2n)2 /3/2m  for  n  >  n.  For  this  reason,  the  non- 
ideality  of  the  electrons  increases  with  increasing  n  as  long  as  n  <  n, 
but  then  decreases  with  increasing  n  as  soon  as  n  i  n.  From  the  condi¬ 
tion  Ze2n1/3  =  E  follows  that  the  electrons  form  again  an  ideal  gas  for 

r 

n  >>  1023  Z3.  This  anomalous  behavior  is  expalined  by  the  stronger  in¬ 
crease  of  E  <*n2/3  with  n  compared  with  the  Coulomb  energy  E  *  n1/? 

r  C 

It  is  recognized  that  the  effects  of  degeneracy  and  nonideality 
on  the  dispersion  of  the  ion  sound  waves,  Eq.  (10),  are  negligible. 
Similarly,  the  effect  of  nonideality  on  the  dispersion  of  the  sound  waves 
of  the  degenerate  electrons,  Eq.  (6),  is  negligible,  but  in  the  disper¬ 
sion  equation  of  the  classical  electrons,  Eq .  (5),  has  to  be  inter¬ 
preted  as  a  polytropic  coefficient,  where  <e(y)  -  Cp/Cv  as  to  order-of- 
magnitude . 
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STATISTICAL  THERMODYNAMICS 

In  the  plasma  under  consideration,  the  electrons  and  ions  interact 

through  their  longitudinal  Coulomb  fields  (transverse  electromagnetic 

2 

interactions  are  negligible  for  KT<<:mc  ).  The  electrons  (s  =  e)  and 

ions  (s  =  i)  have  thermal  velocities  cg  and  random  collective  mean 

mass  velocities  v  due  to  their  oscillatory  wave  motions  about  the 
s 

equilibrium  positions,  so  that  fheir  local  velocity  is  vs  =  +  c8,  with 

>  3 

•-'cs>=  5  and  <vs>  =  ,  where  <Cfg>  =  ///3 s  fs  d  is  the  average  of  ug  with 

respect  to  the  normalized  velocity  distribution  f  of  the  species  s.  The 

resulting  Hamilton  function  with  Coulomb  interaction  leads  to  a  free  energy 

of  the  plasma  of  the  form: 

F  =  I  Fs(0)  +  Em  +  l  Fg  (13) 

s=e,i  s=e,i 


Fg°^  is  the  ideal  free  energy  of  the  noninteracting  plasma  components 
s.  is  the  Coulomb  interaction  energy  of  the  electrons  and  ions  in  their 

equilibrium  positions.  F  .  is  the  free  energy  of  the  electron  and  ion 

e,  l 

oscillations,  i.e.  of  the  high  and  low  frequency  plasmons,  Eq.  (A). 

It  should  be  noted  that  Eq.  (9)  takes  into  consideration  all 

significant  short-range  and  long-range  Coulomb  interactions  by  means  of 

the  Madelung  energy  E^  and  the  plasmon  energies  Fg .  As  is  evident  from 

the  derivation  of  Eqs.  (5}-(6)  and  (10),  in  which  terms  of  order  m/M  are 

neglected  compared  to  1,  Eq.  (9)  contains  the  e-e,  e-1,  and  i-i  Coulomb 

-1/3 

interactions  at  distances  \  >  n 

1.  Free  Energy  Fg°^ .  In  high  pressure  plasmas,  the  electrons  are 

15  3/2  -3 

partially  degenerate  for  densities  n>n  where  n  =  4.828  x  10  T  [cm  ], 

whereas  the  ions  behave  in  general  classically.  Fermi  statistics  gives 

18) 

for  the  free  energy  of  the  ideal  electron  gas 
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< 


e 

where 

u  (u/KT) 
P 

and 


NKT  U3/2(p/KT)/U1/2<u/kT) 

1  f  xpdx 

r(p+D  J  fix- p/kt  +  1 


,  p-l/2,  3/2 


(14) 

(15) 


9  '\l'> 

n  =  2(2irmKT/h  )  ,  U1/2(u/KT) 

19) 

defines  the  Sommer f eld  integrals  ,  and  determines  the  chemical 
potential  y  =  y(n,  T)  of  the  electrons,  respectively.  The  free 
energy  of  the  translational  degrees  of  freedom  of  the  classical 

i  .  ,  18) 

ideal  ion  gas  is 

Fi0)  =  "  (N/Z)  KT  ln[  (27rMKT/h2)3/2  (^j  ] . 


(16) 


(17) 


2.  Quasi-Lattice  Energy  E^.  The  equilibrium  positions  of  the 
electrons  and  ions,  about  which  the  electrostatic  oscillations 
occur,  form  an  electron  "lattice"  and  an  ion  "lattice",  with 
an  incomplete  ordering.  In  the  Wigner-Seitz  approximation, 
the  Coulomb  interaction  energy  of  the  electron-ion  lattices  is, 
independent  of  the  lattice  type, 

Ej,  -  -  ayNKT  ,  «  =  a  =  (4x/3Z)  ,  y  >  1.  (18) 

As  the  ordering  of  the  plasma  increases  with  y,  a(y)  is  a  weak 

function  of  y  such  that  asymptotically  a  =  a  for  y  >>  1.  Eq.  (18) 

2  - 

indicates  that  -Ew/N  -  Ze  / r.  is  of  the  order  of  the  average 


e-i  interaction  energy.  For  weak  ordering,  y  <<  1,  it  will  be  shown 


3.  High-Frequency  Contribution  F  .  Since  the  number  of  longitudinal 
modes  with  wave  numbers  between  k  and  k  +  dk  in  volume  V  is  V4wk2dk/(2n)3 , 
Eq.  (A)  gives  for  the  free  energy  Fg  of  the  high-frequency  electron 
oscillations  of  energy  'fiw(k) 

A 

k 

Fe/KT(V/2ir2)  =  je  in{  1-exp [-ffw(k)/KT] }  k2dk  (19) 


whi_re 


w(k)  =  w  (l+a2k2)  ^2  , 


a2  =  c2/w2  =  (<e/Air)  (Z/y)r2  ,  n  <<  n  , 

*2  s  <3/5>  V''“p  ■  25r<T)1/3(f>2/3  n»»  , 

by  Eqs.  (5)  —  (6) .  The  speed  of  sound  c  and  the  Fermi  speed  v_, 

m  F 


of  the  electrons  are 

cm  =  (<eKT/m)V2  , 


Vp  *  t»(3ir2n)  ^/m. 


The  number  of  modes  in  (0,  kg)  and  V  equals  the  number  3N  of 
degrees  of  freedom  of  the  electron  gas,  i.e.. 


(2tt)-3V 


Airk2dk  =  3N,  k  =  (18ir2n)  '3 
e 


Integration  of  Eq.  (19)  by  parts  yields,  under  consideration  of 

k3  KTV/6ir2  =  3NKT,  for  the  free  energy  of  the  high-frequency  plasmons: 

/  "1/  "too  -v  \ 

F  =  3NKT  £n{l  -  exp[ - 2-  (l+a2k2)  '2]}  -  F( - ak  )  j  C25) 

'  KT  KT  6  / 


where 


F( - s~,  ak  ) 

KT 


J(ak  )“3  ,  x4(l+x2)3dx.  - 

KT  8  0  e  (Ifop/KT)  (1+x2)  2-  1 


/KT  -  (Air)  /2(X Jr  )  (y/Z)  '2  ,  X  -  -R/(mKT)  '2  , 
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ak  =  k  ^2  (9x  ^2/4)  /3(Z/y)  ,  n  <<  n  , 

e  e 

akg  =  2  ^6ir  ^3(  9  / 21^5)  (n/n)  ^(Z/y)  ^  ,  n>>n. 

By  means  of  the  successive  substitutions,  (i)  x=sinh£,  dx  *  cosh  £  d£ 
and  (ii)  e  =  (fiiWKT)  cosh£,  de  =  (-Kmp/KT)  sinh  £  d£,the  integral  (26) 
is  transformed  to 

F(c  ,ak  )  =  (ak  e  )  3  r&  (e2-e2  )  ^2(eE-l)  1  de 
p  e  e  p  J  p 


where 


e  =-Tuo  /KT,  £  =  e  [1  +  (ak  )  ]"*  . 

"  ~  e  p  e 


P  P 


Since  the  leading  expression  in  Eq.  (25)  is  the  logarithmic  term,  it 
is  sufficient  to  give  for  F(ep,  akg)  series  approximation  (Appendix), 

F(ep,ake)/23/2(ake)3ep~3/2  = 


I  e  m£P  I  (n)(2cJ  n  m  (2  +  n\d  +  n,  (e  -e  )m),  e  <  3e  ,  (32) 


n  p 


e  p  e  p 


where 


5  e  -e  3 

,  -  r  +  nep^  +  n 

,5  .  .  v  2  f  2  -mu, 

Y  (-  +  n,  (e  -e  )m)  =  m  J  u  e  du 


is  the  incomplete  gamma  function—^ ,  Since  in  general  y/Z  *  1  for 


ep  <  eg  <  3ep,  the  expansion  (32)  is  useful  where  simple  approximate 
relations  do  not  exist. 
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4.  Low  Frequency  Contribution  F^.  With  the  number  of  modes  in  the 
interval  dk  at  k  and  volume  V  given  by  V4irk2dk/ (2r) 3 ,  Eq.  (4)  yields  for  the  free 
energy  F  of  the  low-frequency  ion  oscillations  of  energy  fim(k) 


F  /KT(V/2x2) 


ln{  1-exp  [-1S(o(k)/KT]}  k2dk 


(34) 


where 


u(k)  =  6 (k)cMk, 


cM  =  (k.KT/M)  12  . 


by  Eqs.  (10)  and  (12).  The  number  of  modes  in  (0,  k^)  and  V  equals 
the  number  ^N/Z  of  degrees  of  freedom  of  the  ion  gas,  i.e., 

h 

(2ir)"3V  l  4Trk2dk  =  3N/Z,  k±  =  (18ir2n/Z) 1/3 


(35) 

(36) 


(37) 


Partial  integration  of  Eq.  (34)  gives,  under  consideration  of 
k^  KTV/6ir2  =  3 (N/Z)KT,  for  the  free  energy  of  the  low-frequency  plasmons: 


-tic. 


Fi  =  3 (N/Z)KT  (  ^n{ 1-exp [- 


M 


KT 


6(k.)k1]} 


G<v) 


(38) 


where 


-fic, 


G(k±)  = 


M  ^  -3  f  1  [  6  (k)+k<5 '  (k)  ]k3dk 

KT  1  1  e(ficiyKT)6(k)k  _  1 


(39) 


Since  the  dispersion  factor  6 (k)  is  a  bounded  function  varying  very 
little  with  k  such  that  l<6(k)  ^  (1+Z)  ^2  for  k  e  (0,  k^) ,  6(k)  can 

be  approximated  by  an  average  value  5  , 

6  (k)  =  6  -  1,  n  <  n  .  (40) 

Since  in  addition  the  logarithmic  expression  is  the  dominant  term  in  Eq.  (38) > 


the  integral  (39)  is  approximated  by 

e 


G(£a) 


T 


e3(eE-l)-1  dt 


(41) 
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where 

e  =  "RCj^Sk/KT,  «=  -fic^fik^/KT. 

G(ei)  has  the  semi-convergent  series  expansions, 
G(ei)  =  +  joei  £i  <<  1* 

G(Ci)  =  iy  ei  3  +  °le  Eil»  >:>  !• 


20) 


(42) 


(43) 

(44) 


This  completes  the  formal  mathematical  aspects  of  the  theory, 
the  physical  implications  of  which  require  further  elaboration. 
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APPLICATIONS 


For  applications  of  the  theory  to  strongly,  intermediate,  and  weakly 

nonideal  plasmas,  it  should  be  noted  that  the  dimensionless  parameters 

Y/Zj'fiw  /KT,  ak  ,  and  n/n  occuring  in  Eq.  (25)  for  the  free  energy  F 
p  e  e 

of  the  high-frequency  plasmons  can  not  be  varied  independently.  Since 
both  y/Z  and  Xe/rg  increase  with  increasing  n  and  decreasing  T, 

Tuitp/KT  ~  (Xg/re)(Y/Z)  ^varies  over  a  large  n-T  region  similar  to  (y/Z)  1 , 

Eq.  (7).  Numerically, 

Y/Z  =  1.670x  KrV^/T,  1wp/KT  =  4. 328>clO"7n1^  /T,  n/n  =  2 . 071xl(f  16nT“^  , 
ake  =  1.586  Ke^(Y/Z)"^  ,  n  <<  n 

-  V,  _U 

akg  =  3. 308 (n/n)  3  (y/Z)  12  ,  n  »  n.  (45) 

E.  g.,  for  T=104°K,  y/Z  £  1  if  n£  1021cm~3  and  1&o  /KT  £  1  if 

n  >  5xl020cm“3.  For  T=103°K,  y/Z  <  1  if  n  <  1018cnf3,  etc.  Thus, 

for  typical  conditions  of  nonideal  plasmas  y/Z  and  T5u>  /KT  are  of  the 

P 

same  order  of  magnitude.  It  is  also  recognized  that  in  general 
n/n  >>  1  if  y/Z  »  1,  and  n/n  <<  1  if  y/Z<<  1. 


In  Eq.  (38)  for  the  free  energy  F^  of  the  low  frequency  plasmons, 
only  one  characteristic  parameter  occurs  since  6(k)~  6-1.  By  Eq. 
(42),  this  parameter  is 


fiVki 


X/2'  ji 


(187T2)  3IC 


i  «  — 
r , 


=  2.158x10'V1/3(^)1/2^/'6«1 


where 


= -ti/ (MKT) 1/2  ,  r±  =  (n/Z)^3  .  (47) 

Accordingly,  for  typical  nonideal  plasma  conditions,  it  is  <<  1 

since  X^/ri<<l  (classical  ions)  although  in  general  Xg/re  >  1  (degenerate 

electrons)  for  y/Z>l  or  liuj  /KT  >  1. 
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The  deviation  AF  of  the  free  energy  of  a  nonideal  plasma  from 


ideality  is  by  Eq.  (13)  due  to  the  quasi-lattice  energy  E^  and  the  plasmon 
energies  ^ , 

AF  =  Em  +  I  Fs-  (48) 

s=e,i 

Since  the  theory  of  electron  oscillations-^-^ has  not  yet  been 
developed  for  arbitrary  degrees  of  degeneracy  (n  >  n),  the  contributions 
of  the  electron  oscillations  to  AF  in  the  cases  n~  n  and  n~  n  have  to 
be  estimated  from  the  dispersion  equations  for  n  <<  n  [Eq.  (5)]  and 
n  >>  n  [Eq.  (6)],  respectively.  Fortunately,  it  turns  out  that 
|  F^  1  <  <  |  AF  |  for  -y/Z>l,  so  that  quantitatively  relyable  approximations 
for  AF  can  be  derived. 


1.  Strongly  Nonideal  Plasmas.  By  Eq.  (6)  the  spectrum  w(k)  of 
electron  oscillations  extends  over  a  band  Am  -  u  above  the  plasma 

-  <  '  -  7k  -1 

frequency  for  y/Z  »  1  since  kre~  kgr  -  1  and  (n/n)  3Zy  ~  1.  Application 
of  the  mean  value  theorem  for  integrals  to  Eq.  (25)  shows  that  the  free 
energy  Fg  of  the  high-frequency  plasmons  vanishes  exponentially  for 


i.e.  y/Z 


F  /3NKT 
e 


Zn{  1  -  exp[-e  (l+a2k  2)  ^2]} 
l  p  e 


e  (ak  ) 

_E _ 


-3 


ak 


exp[Ep(l+x2) V2 ]-l 


x^(l+x2)-  ^2dx 


)*°’ 


0  -  x  -  ak 


(49) 


Accordingly,  |Fe|/3NKT<<l  for  ep  »  1,  i.e.,  y/Z  »  1.  On  the  other 
hand,  the  free  energy  of  the  low  frequency  plasmons  is  by  Eq.  (38)  for 
nondegenerate  ions 
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F.  =  3 (N/Z)KT[£rte  -(1/3)]  « 

11  1/ 

3(N/Z)KT{£rtY  +  ln{  (18tt2/Z4)  1/3  <*1KT/M>  2]  _  (1/3)},  e  «  l.  (50) 

e2/tf  1 

It  is  noted  that  y  /Z  >>  1  is  compatible  with  =  fic^^iS /KT  <<  1  as 
explained  above. 

Equations  (49)  and  (50)  demonstrate  that  the  contribution  of  the 
electron  oscillations  to  the  free  energy  is  negligible  in  strongly 
nonideal  plasmas,  y/Z  >>  1.  In  this  limit,  the  nonideal  part  of  the 
free  energy  is  due  to  the  quasi-lattice  energy  and  the  ion 
oscillations, 

AF/NKT  =  -  ay  +  (3/Z)lny  +  (3/Z)ln (Bc^/Vg)  -  (1/Z) ,  y/Z  »  1,  (51) 

where 

vB  =  e2/fi,  6  =  (187r2Z-4)1/3  .  (52) 

Note  that  Iny  depends  on  both  n  and  T  whereas  ^nSc^/Vg  depends  only 

8  1/ 

on  T,  where  the  Bohr  speed  is  vfi  =  2.118  x  10  cm/sec  >>  -  (k^KT/M)  * . 
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It  is  remarkable  that  the  electron  oscillations  contribute  little 

to  the  free  energy  compared  to  the  ion  oscillations  for  y/Z  >>  1.  This 

result  holds  even  for  moderately  nonideal  conditions,  y/Z  >  1.  Thus, 

we  disagree  with  the  formula  "F  =  ne  +  3NKT&1  (Hw  /KT)"  stated  without 

o  o 

21 

derivation  for  nonideal  plasmas  by  Norman  and  Starostin— ,  according 
to  whom  "all  the  vibrations  have  exactly  the  same  frequency  w 
near  the  plasma  frequency  wp".  The  derivation  of  this  formula  requires 
ffw(k)/KT  <<  1  for  the  electron  oscillations,  which  implies  y/Z  <<  1, 
but  the  latter  inequality  contradicts  their  assumption  w(k)  =  wq  =  wp, 
since  the  frequency  spectrum  extends  over  a  large  Land  Aw  >  wp  above 
wp  for  y/Z  <<  1.  For  these  reasons,  the  free  energy  proposed  by  them 
is  not  applicable  to  proper  nonideal  plasmas,  y/Z  >  1,  nor  is  it  correct 
for  less  nonideal  conditions,  y/Z  <  1. 

2.  Intermediate  Nonideal  Plasmas.  For  intermediate  nonideal 

conditions,  1  f  y/Z  <  10,  the  spectrum  w(k)  of  electron  oscillations 

'  2/  -1 

extends  over  a  region  Aw  <  0(w  ]  above  w  by  Eq.  (6)  since  (n/n)  3Zy  <1  and 

P  P 

kr^l  kfire '  1.  Also  in  this  case,  a  relatively  simple  formula  can  be 
devised  for  the  free  energy.  The  logarithmic  term  in  F  ,  Eq.  (25) 
is  negligible  compared  to  that  in  F^,  Eq.  (38) »  for  y/z  >  1  since 
ep  >>  for  y/Z  >  1  by  Eos.  (45)  and  (46),  respectively. 

Accordingly,  the  nonideal  part  (48)  of  the  free  energy  is  for  inter¬ 
mediate  nonideal  plasmas: 

AF/NKT  =  -ay  +  (3/Z)£ny  +  (3/Z) In (6c„/vn)  -  (3/Z )G(P  ) 

M  D  X 

-  3F(ep,ake),  y/z  >  i.  (53) 
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For  yjz  >  1,  the  ions  can  be  assumed  to  be  non-degenerate  ,  =  "Tic^Ak^/KT  <<  1 

by  Eq.  (46),  so  that  the  ion  integral  (41)  reduces  to 

GU±)  =  1/3,  «  1  .  (54) 

Since  e  >  e  ?  1  and  ak  e  ?  ]  (Eq.  (45)1  for  y/Z  *  1,  the  electron 

e  p  g  p 

integral  (30)  is  significantly  smaller  than  G(e^)  =  1/3, 


0  <  F(e  ,ak  )  <  (e2-e2)  ^  (ak  e  )  3£n(l-e  £e/l-eCF)  <<  1,  y/Z>  1. 
peep  e  p  ' 

The  lower  and  upper  bounds  of  F(e^,ake)  have  been  obtained  by  means  of 


the  mean  value  theorem  for  the  integral  (30), 

-3  . 


F(e  ,ak  )  =  (ak  e  ) 
p  e  e  p 


3/  % 

h  r,e 


'(t2- E2)  /2  J( 


e  -l)de. 


e  <  e  <  rA 

p  -  e 


(55) 


(56) 


While  for  strongly  nonideal  conditions,  the  contribution  of  the 
electron  oscillations  to  the  free  energy  is  completely  negligible, 
this  contribution  is  still  insignificant  for  intermediate  nonideal 
conditions,  y/Z  >  1,  by  Eq.  (55).  For  more  exact  evaluations,  the 
small  term  F(Ep,ake)  in  Eq.  (53)  can  be  computed  from  Eq.  (30)  or  (32). 

3.  Weakly  Nonideal  Plasmas.  Although  the  theory  of  weakly 
nonideal  systems  is  well  understood,1  5^it  is  interesting  to  investigate 
whether  the  present  model  for  proper  nonideal  plasmas  gives  reasonable 
results  in  the  limit  y/Z  <<  1.  For  y/Z  <<  1  it  is  akg  >>  1  by  Eq.  (45), 
and  the  spectrum  w(k)  of  electron  oscillations  extends  over  a  large 
region  Am  >>  above  by  Eq.  (5).  The  electron  integral  becomes  for 


ak  >>  1, 
e 


F (e,ak)  =  e  (ak  ) 
p  e  p  e 


-3 


ak 


/ 


(eGPX-l)  ‘x3dx,  y/Z  <<  1, 


i.e.  , 


F(En’akp)= 1'T^’;r,alco)  +  )  -  ••.),  e  ak  <<  1. 

pe  Jope  20  pe  ’pe 


(57) 


(58) 
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Although  e^akg  is  independent  of  y/Z  by  Eqs.  (27)  and  (28),  the  expansion 
(58)  is  valid  since  the  electrons  are  certainly  nondegenerate,  Xg/r  <<  1 


for  y/Z  <<  1 ,  and 


,  V?  /n_  '4>  /  /  \  V;) 


e  ak  =  (4itk  )  '2  (9x  2/4)'3  X  /r  <<1,  X  /r  <<1. 
pe  e  ee  ’ee 


For  nondegenerate  ions,  the  integral  (41)  is  G(e^)  =  1/3  by  Eq.  (43) 
since  <<  1.  Thus,  one  obtains  from  Eqs.  (18),  (25)  and  (38)  for 
the  interaction  part  of  the  free  energy  of  weakly  nonideal  plasmas: 


AF/NKT  =  -a(y)y  +  (3/Z)£iiy  +  (3/Z)£n(6cM/vB) 

+  3£n(epake)  -  (1+Z-1),  y/Z«l,  (60) 

where  the  logarithmic  term  in  Eq.  (25)  has  been  expanded  for 

c  ak  <  <  1 . 

P  e 

In  Eq.  (60),  a(y)  is  the  Madelung  constant  of  the  weakly  nonideal 
plasma  with  weak  electron  and  ion  ordering,  a(y)->-0  for  y-*-0.  Comparison 
of  the  term  -a(y)y(NKT)  in  Eq.  (60)  with  AF  = 

-(NKT)  (2 / 3) tt  ^2(1+Z)  ^2e3n'2(KT)  ^2  of  the  Debye-Hueckel  theory—^ 

(weakly  nonideal  plasmas)  yields  the  result 

a(y)  =  (2/3)tt1/2  (1+Z_1)3/2  y^2  ,  y/Z  <<  1.  (61) 

The  previous  theories  of  weakly  nonideal  plasmas  do  not  lead  to  the 
logarithmic  terms  in  Eq.  (60)  since  they  do  not  take  into  account  the 
effects  of  electron  and  ion  oscillations. 

4.  Numerical  Illustrations.  Fig.  1  shows  the  (negative)  free  energy  Fq 
of  an  ideal  Z  =  1  plasma  versus  n  and  T  based  on  Eqs.  (14)  -  (17).  F  Serves 
as  a  reference  quantity,  relative  to  which  the  quantitative  significance  of  the 
nonideal  contributions  are  measured.  If  I  increases  with  increasing  n  and  T. 


133 


Fig.  2  shows  the  deviation  AF<0  of  the  free  energy  of  a  Z  *  1 
plasma  from  its  ideal  value  Fo<0  versus  n  and  T  based  on  Eqs. 

(48),  (25)  and  (38).  In  the  n-T  region  under  consideration,  )AF|  is 

of  the  same  magnitude-of-order  as  | FQ j ,  i.e.  is  considerably  larger  than 

the  thermal  energy  -NKT.  AF/F0  exhibits  only  at  large  densities 
19  -3 

n  >10  cm  a  significant  T  -  dependence. 

Fig.  3  shows  the  free  energies  Fe  and  F^  of  the  high  (e)  and  low  (i) 

frequency  plasmons  of  a  Z  =  1  plasma  based  on  Eqs.  (25)  and  (38).  |F  |  is 

considerable  larger  than  |Fe|,  in  particular  at  higher  densities.  The 

T-dependence  of  F  ./  F  increases  with  increasing  density  n.  Comparison 

of  Figs.  2  and  3  indicates  that  AF  -  F  +  F  ,  i.e.  the  quasi-lattice  energy 

[Eqs.  (18)  and  (61)]  is  not  the  dominant  nonideal  effect. 

The  Figs.  2  and  3  demonstrate  the  quantitative  importance  of  the  nonideal 

effects  AF  =  Ej^  +  Fe  +  ,  in  particular  of  the  low  (i)  and  high  (e)  frequency 

plasmon  contributions  F.  and  Fa  (F,  >  F  ) ,  for  the  evaluation  of  the  free 

1  e  1  e 

energy  F  =  FQ  +  AF  of  high  density  plasmas. 

For  quantitative  calculations,  it  is  noted  that  the  free  energy  AF  is 
hardly  affected  by  inaccuracies  in  the  large  maximum  wave  numbers  ke  and  k^, 
which  have  been  determined  in  accordance  with  the  Debye  theory  which  implies 
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Free  energies  F  <  0  of  high  (e)  and  low  (i)  frequency  plasmons 

e  »  i  Q 

versus  n[cm  ]  with  T[  K]  as  parameter  (Z=l) • 


APPENDIX:  Expansion  of  F(e  ,ak  ) 

P  e 

The  integral  (23)  is  conveniently  rewritten  in  the  form 

F(e  ak  )  *  (ak  e  )-3  I(e  ,e) 
p  e  e  p  p 


where 


I(e 


/e)  -  j  (c2-e2) 


3/  _  1 

'z(ee-l)  de,  0<e  <  e  <  °°. 

P 


(Al) 


(A2) 


Since  e  >  0,  i.e.  e  <1,  there  exists  the  series  expansion, 

-  J  t  >0  . 

m=l 


(A3) 


The  substitution,  u  -  e  .  du  =  de,  and  Eq.  (A3)  transform  Eq.  (A2)  to 


I(e  ,e)  *  l  e"mep  7  "V'2  (u+2e  )%  e'mudu. 

P  rl  J  P 


(A4) 


u=o 


For  u<2ep,  i.e.,  e  <  3ep,  the  binomial  expansion, 

X  T  u  ■ %n 


(u+2ep)  /2  =  (2ep)  '2  £  (£)(£-)*,  u/2ep  <  1, 

n=0  P 


(A5) 


is  used,  which  reduces  Eq.  (A4)  to  the  double  series: 


I(£  ,e)  “  (2e  )  /z  £  e~mEp  J  (  n  >  (2 O'"  ®  2  +nV(-|-+  n,  (e-e  )m) , 


m*l  n~0 


£  <  3e  , 
P 


(A6) 


where 


5  ,  e-e  3  . 

c  -  -=-  +  n  ..  p  -*-  +  n 

.5  .  ,  ,  ,  2  C  2  -mu 

Y(-x-+n  ,  (e  -E„)m)  *  m 
1  P 


e  du. 


(A7) 


0 


20) 

is  the  incomplete  gamma  function,  which  is  tabulated.  In  an 


amologous  way,  the  integral  (A2)  can  be  solved  for  u  >  2e  ,  i.e., 

P 

3e  <  e  <  =>. 

P 
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VIII.  FREE  ENERGY  OF  RANDOM  SOUND  OSCILLATIONS 


By 

A.  H.  Khalfaoui  and  H.  E.  Wilhelm 


ABSTRACT 


Thermal  equilibrium  properties  of  a  monatomic  gas  are 
investigated  by  taking  into  account  the  energies  of  the  random 
sound  wave  oscillations.  The  free  energy  is  derived  by  a  quan¬ 
tum  statistical  mechanics  due  to  Bose.  The  system  is  considered 
as  a  macroscopic  continuum  in  which  random  acoustic  oscillations 
are  thermally  excited.  It  is  shown  that  the  contribution  AF  to 
the  free  energy  due  to  the  sound  w-’ves  is  significant  for  high 
density  gases,  in  particular  at  moderately  high  temperatures. 
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INTRODUCTION 


In  solids  and  liquids,  the  effect  of  sound  waves  on  the  theraodyntalc 

quantities  was  studied  by  Landau  and  others^.  In  ionized  gases,  the 

electrical  oscillations  (plasma  oscillations)  affect  the  thermal  equilibrium 
(2-3) 

of  the  system  .  Similarly,  we  are  considering  the  contribution  of  the 
sound  wave  oscillations  to  the  free  energy  of  noncondensed  gases.  In 
thermal  equilibrium  of  gases,  the  acoustic  oscillations  share  the  partition 
of  energy  and  .thus,  change  the  thermodynamic  Functions  of  the  system.  The 

distribution  of  the  sound  wave  quanta  is  determined  by  Bose  statistics, 
which  is  used  herein. 

The  problem  under  consideration  is  concerned  with  gases  as  a  macroscopic 
continuum,  which  exhibits  a  set  of  separate  elementary  exitations,  the  sound 
wave  oscillations.  These  exitations  behave  like  “’quasi-particles"  moving 
in  the  volume  occupied  by  the  gas,  and  have  definite  energies.  The  free 
energy  of  the  gas  evaluated  by  the  theory  to  be  presented  will  take  into 
consideration,  in  addition  to  the  random  thermal  energies  of  the  gas 
particles,  the  energy  of  the  random  sound  wave  oscillations.  It  will  be  shown, 
that  the  effect  of  the  sound  waves  is  important  only  at  high  temperatures 
and  high  gas  densities.  The  results  of  this  theory  are  applicable  at 
temperatures  and  densities  for  which  the  gas  is  not  in  a  condensed  state  t 

(liquid  or  solid) . 

Although  nonideal  effects  due  to  finite  particle  size  are  not  taken  into 
account  explicitly,  it  should  be  noted  that  the  gas  under  consideration  is  not 
an  ideal  one.  The  existence  of  sound  waves  in  the  gas  implies  that  there 
are  particle  interactions,  since  a  gas  can  not  perform  the  ordered,  collective 
mean  mass  motions  of  random  sound  waves  without  such  interactions. 
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THEORY 


Consider  a  gas  as  a  continuum  of  volume  V  containing  N  atoms.  Because 
the  velocity  of  the  gas  in  a  sound  wave  is  in  the  direction  of  propagation, 
the  sound  waves  are  longitudinal.  Each  oscillator  of  frequency  of  the 
longitudinal  sound  waves  can  only  have  the  energies  (fi  »  h/2ir  - 

reduced  Planck  constant) 


E  -  "fi  u)  (n  +-r),  n  *  0,1,2, ...»  . 
n  ao2  o 

o 

The  frequency  of  the  sound  waves  with  wave  number 

a)  =  k  C  ,  C  =  (yKT/M)*4 
o  o  s  s 


where  y  *  cp/cv  is  the  poly tropic  coefficient, 
function  of  the  gas  oscillations  is: 


(1) 

k^  is  (M  *  mass  of  atoms) 

(2) 

Accordingly,  the  partition 


ri  r  -ghC  k  (n  +  -=r)  n 

V  L  e  so  o  2 
n  =0 


-Bfic  k  n 

e  so 

.  -&KC  k 
1-e  s  o 


(3) 


o 

where  3*1 /KT  (K  *  Boltzmann  constant,  T  *  Temperature  of  the  system).  From 
the  partition  function  Z,  the  thermodynamic  quantitites,  such  as  pressure, 
internal  energy  etc.  »  are  derived  in  the  usual  way.  The  free  energy  of  the 
random  sound  oscillations  is  given  by: 

AF  *  -KT  in.  Z  •  (4) 

In  the  limit  V  -*■  ®,  the  discrete  eigenfrequencies  are  replaced  by  a  continuous 
spectrum,  u>  -  w(k),  in  accordance  with  the  dispersion  law  for  sound  waves,  of 
wave  length  X  *  2ir/k, 

u  *  k  C  ,  0<k<k  (5) 

s 

The  theory  to  be  presented  is  sensitive  towards  the  cut-off  wave  number 

k,  which  is  large  in  all  cases  of  interest.  Since  acoustic  waves  with  wave- 

-  -1/3 

lengths  A<max(r,L)  and  mean  free  paths  L<r*n  are  not  possible  in  gases. 
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k  =  2” / X  is  determined  by  the  mean  free  path  L, 


k  =  —  ,  L  =  L(n,T) ,  (6) 

where  n  =  N/V  is  the  density  of  the  atoms.  The  number  of  (longitudinal)  wave 

2  3 

modes  with  wave  numbers  between  k  and  k  +  dk  in  volume  V  is  g(k)dk  ■  V  Air  k  dk/(2n)  . 
Accordingly,  Eqs.  (3)  and  (4)  give 


AF  =  -  f  k 2lnh - 

2tt  J  \l-e 


-pt(Csk/2 

-fiffC  k 
-e  s 


The  integral  in  Eq.  (7)  is  decomposed  into  the  contributions  i)  from  the  ground 

state  (n  =  0)  and  ii)  higher  states  (n  >0).  By  Eqs.  (3)  and  (7) 

a  o 

AF  =  Fl  +  F2  (8) 


where 


VfrC  F  "i 

F  =  - 1  1  k  dk  =  VfiC 

1  4n2  J 

o 


©' 


F  ,  Vki_  t  _  -9KC  k  _  r  k^dk _ 

F2  2  6  ’  ,  2  I  etfc  k  U°' 

66tt  6  tt  J  e  s  -  1 

o 

by  partial  integration.  The  integral  in  Eq.  (10)  can  be  solved  for  "high"  and 

(4  5) 

"low"  temperatures  by  series  expansions)  which  give  ’  : 


F  =Xk 1 
2  2 


f  -x.  1  x.  ^2vX 

en(l  -  e  )  -  3  +  8  +  x  I  (2v  + 


6 3tt2 V  '  '  3  8  v=l  (2v  +  3)2v 


,  (  »  x  <  2  IT 


F  =Vk  fH(1.e-x  1. 
2  68tt2  x3 


*■  n=l  \  nx  nx  nx/Jf 


x  >  1 


where 


x  =  3trc  k 
s 


!  i-  >  ; 


and  are  the  Bernoulli  numbers,  and  c(4)  *  it  /90  is  the  Riemann  ^-function. 

For  comparison  purposes,  AF  and  the  classical  free  energy  F  of  the  ideal 

o 

monatonic  gas  (M  =  atomic  mass)  are  stated: 


F  =  NKT 
o 


H(^  )32  n(  ■  *] 


A  F  =  VffC  1-^-1  + 


}  ) 

\  nx  nx  n  x  /  I  I 


(2v  +  3)  2v! 


x<2ti.  (15) 


6c  (A)  -  x4  l  e  nx 


x  >  1 


for  "high"  and  "low"  temperatures,  respectively. 

It  is  interesting  to  compare  the  specific  heat  of  the  ideal  gas  (Cq)  and 
the  sound  oscillations  (AC)  in  the  high  temperature  limit,  x  =  Sh^k  <<  1. 

By  Eqs .  (14)  and  (15)  , 

C  =-T82F  /3T2  =  3NK/2  (17) 

o  o 

AC  =  -T32AF/3T2-  J^(V/L3)K  <<  Cq  for  r  «  n_1/3  <<  L  .  (18) 

~3  2 

In  the  derivation  of  Eq.  (18),  it  should  be  noted  that  AFs (Vk  /6Bn  )ln  x 

-1/2 

for  x  <<  1  where  x  «  T  .  It  is  seen  that  AC  <  Cq  at  high  temperatures  since 
r  <  L. 

This  completes  the  mathematical  aspects  of  the  problem,  the  physical 


implications  of  which  will  be  discussed  next. 
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DISCUSSION 


In  a  hypothetical  ideal  equilibrium  gas  without  particle  interactions, 
the  average  particle  velocity  is  <c>  =  j  cf(c)  d^c  =  15 ,  i.e.  the  particles  have 
pure  thermal  velocities  c  with  a  Maxwell  distribution  f(c).  No  random  mean  mass 
motions  or  collective  particle  motions,  such  as  sound  oscillations,  exist  due 
to  the  absence  of  parLicle  interactions.  The  free  energy  of  the  ideal  or 
non-interacting  monatomic  gas  is,  therefore,  F  ,  Eq.(14).  In  any  real  gas 
with  particle  interactions,  stochastic  mean  mass  motions  <v(r,t)> 

=  j  v  f(v,r,t)d  v  ^  0  exist  due  to  the  presence  of  thermally  excited  sound  waves 

*>  “►  _ 
(f(v,r,t)  is  the  local  distribution  of  actual  particle  velocities  v  =  <v>  +  c]. 

Since  the  total  energy  of  the  gas  is  distributed  both  over  the  thermal  particle 

-►  “►  -y 

motion  c  and  the  stochastic,  acoustic  mean  mass  motions  <v(r,t)>  ,  a  free 

energy  contribution  from  the  random  sound  waves  exists.  Thus,  the  free  energy 
IF  of  the  random  sound  oscillations  represents  a  nonideal  effect  which  is 
ultimately  due  to  particle  interactions,  which  make  a  hydrodynamic  or  continuum 
description  of  a  gas  possible. 

In  the  contribution  AF  of  the  sound  waves  to  the  free  energy  F  of  an 

ideal  gas  as  given  by  Eq.Q5)  or  (16),  we  identify  two  parts.  i)F^  which  is 

the  contribution  of  the  "zero  oscillation"  mode  which  corresponds  to  n^=0  in 

Eq.(l),  and  ii)  F  the  higher  mode  oscillation  contributions  n  >1.  The 
l  o 

explanation  for  the  increase  of  the  free  energy  of  the  sound  quanta  'flu)  with 
temperature  is  given  by  statistics.  In  the  high  temperature  limit,  the  number 
N  of  sound  quanta  of  frequency  o>  is  N  ?  KT/tfu;,  i.e.  increases  proportional  with  T. 

In  Fig.  1  we  have  drawn  AF/Fq  for  monatomic  helium  gases  over  a  range  of 
temperatures  and  densities  to  show  the  variation  of  AF.  The  overall  contribution 
of  AF  is  larger  at  higher  densities  but  decreases  rapidly  for  lower  densities, 
especially  at  high  temperatures.  Quantitatively,  AF  represents  a  noticeable 
effect  only  at  extremely  high  densities  n  of  the  gas  [n>1021cm3).  For  this 
reason,  the  free  energy  AF  of  the  sound  oscillations  has  to  be  considered  in’the 
evaluation  of  the  thermodynamic  functions  of  high-temperature  gases  only  at 

h  i  eli  dens  i  r  i  os  . 


Cases  with  a  considerable  acoustic  noise  background  are  encountered 
in  various  high  temperature  engineering  systems,  such  as  gas  turbines, 
jet  engines,  rocket  exhausts,  etc.  The  theory  presented  permits  calcula¬ 
tion  of  the  free  energy  AF  of  the  acoustic  degrees  of  freedom  in  such 
systems,  provided  that  the  acoustic  noise  is  in  thermal  equilibrium. 
Considerably  larger  free  energy  contributions  are  to  be  expected  under 
nonequilibrium  conditions,  particularly  if  the  acoustic  fluctuations 
exhibit  intensity  levels  corresponding  to  turbulence. 
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